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Magnetic  circular  dichroism  (MCD)  and  absorption  spectroscopies  have  been  used 

to  study  Ag,  Au,  Mo,  Ti  and  Co  atoms  trapped  in  rare  gas  matrices.    Band  moment 

analyses  have  been  performed  for  the  2P  *-  2S  transition  of  Ag  and  Au  atoms  and  matrix 

induced  changes  in  their  spin-orbit  (SO)  splitting  have  been  assessed.  The  614  cm"1  SO 

splitting  of  gas  phase  Ag  atoms  is  found  to  increase  to  783  cm-1  and  638  cm"1  in  Ar  and 

Kr  matrices,  respectively.  A  reduction  to  524  cm"1  is  observed  in  Xe  matrices.  For  Au 

atoms,  the  observed  matrix  splitting  (Ar:  3354  cm"1,  Kr:  3110  cm"1,  Xe:  2655  cm"1)  is 

greater  than  the  gas  phase  in  all  cases.  These  experimental  values  are  compared  to  values 

predicted  from  a  model  in  which  the  isolated  metal  atom  is  considered  to  occupy  a 

tetredecahedral  hole  formed  by  twelve  nearest  neighbor  rare  gas  atoms.     A  similar 

7  "7 

analysis  for  the  z    P  «-  S  transition  of  Mo  atoms  revealed  an  increase  in  Ar  matrices 
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(126  cm-1)  from  the  gas  phase  value  of  101  cm"1.  Reductions  to  84.5  cm-1  and  57.5 
cm"  were  noted  in  Kr  and  Xe  matrices,  respectively.  The  model  is  extended  to  the  many 
electron  system  of  Mo  atoms  with  good  results.  The  analysis  for  all  three  atoms  also 
revealed  considerable  Jahn-Teller  activity  in  the  excited  state  which  for  the  Ag  and  Au 
atoms  results  in  a  further  splitting  of  the  2P3/2  level.  This  results  in  the  observation  of 
three  absorptions  in  matrices  where  only  two  are  seen  in  the  gas  phase. 

The  observed  transitions  for  matrix  isolated  Ti  and  Co  atoms  have  been  assigned 
by  comparing  the  experimental  spectra  to  simulations  produced  from  gas  phase  data.  For 
Ti  atoms  in  an  Ar  matrix,  this  was  found  to  be  straightforward  with  the  predicted 
spectrum  matching  the  experimental  one.  For  Co  atoms  in  Xe  matrices,  several  observed 
MCD  bands  were  found  to  have  opposite  sign  from  the  predicted  bands.  Further 
investigation  revealed  that  three  bands  had  been  misassigned  in  the  original  gas  phase 
work.  The  same  MCD  patterns  were  observed  in  Kr  and  Ar  matrices,  but  multiple 
trapping  sites  complicated  the  assignment  of  transitions. 
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CHAPTER  I 
INTRODUCTION 


The  ability  of  magnetic  circular  dichroism  (MCD)  spectroscopy  to  provide 
information  about  excited  states  has  produced  much  interest  in  the  technique.  The  excited 
states  of  highly  reactive  species  may  be  probed  using  MCD  if  they  can  be  trapped  in  a 
frozen  rare  gas  matrix.  Magnetic  circular  dichroism  studies  of  the  electronic  and 
structural  configurations  of  small  transition  metal  clusters,  which  are  becoming  increasing 
important  in  catalysis  and  microelectronics,1  have  been  hindered  by  an  incomplete 
understanding  of  matrix  induced  changes  in  the  spectra  of  trapped  species.  The  study  of 
matrix  isolated  atoms  to  characterize  and  quantify  these  matrix  induced  spectral  changes 
is  necessary  before  MCD  characterization  of  clusters  may  begin. 

The  introduction  of  matrix  isolation  by  Whittle,  Dows  and  Pimentel  in  1954  has 
allowed  a  large  number  of  unstable  or  highly  reactive  species  to  be  isolated  in  a  frozen 
gas  lattice  and  studied  spectroscopically.  Many  standard  spectroscopic  techniques, 
including  infrared,  UV- Visible  absorption  and  emission,  Raman,  Mossbauer  and  electron 
spin  resonance,  have  been  used  to  probe  these  molecules.  As  a  first  approximation, 
interactions  between  the  matrix  and  sample  are  assumed  to  be  nonexistent  and  the 
observed  spectra  are  expected  to  be  equivalent  to  gas  phase.  By  comparing  matrix  spectra 
for  species  where  gas  phase  data  are  available  a  number  of  matrix  induced  spectral 
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2 
changes  have  been  noted.   Typically  there  are  slight  shifts  in  transition  energies  when  a 

species  is  trapped  in  a  matrix,  but  these  generally  do  not  prevent  correlation  with  gas 

phase  data.  Multiple  trapping  sites  (different  micro-environments  within  the  matrix)  may 

result  in  overlapping  bands.  These  sites  can  often  be  removed  by  careful  warming  of  the 

matrix  (annealing)  to  allow  the  crystal  lattice  to  rearrange  to  a  more  stable  configuration. 

Jahn-Teller  distortions  and  lattice  mode  interactions  may  also  lead  to  changes  in  the 

matrix  spectra. 

The  ""P  *-  "S  transition  of  the  noble  metals,  Cu,  Ag  and  Au,  has  been  extensively 
studied  in  matrices.  The  transition  is  particularly  sensitive  to  the  matrix  environment. 
In  the  gas  phase  the  2P  state  is  spin-orbit  (SO)  split  into  2P3/2  and  2P1/2  multiplets20 
resulting  in  two  observed  bands.  Trapping  these  molecules  in  the  various  rare  gases  may 
produce  profound  changes  in  the  magnitude  of  the  SO  splitting.  For  Cu  in  matrices  the 
SO  constant  has  been  shown  to  drop21,  from  the  gas  phase  value  of  166  cm-1  to  124 
cm-1  in  At,  95  cm-1  in  Kr  and  to  -23  cm-1  in  Xe.  Additionally,  the  matrix  spectra 
invariably  exhibit  three  absorption  maxima  attributable  to  a  splitting  of  the  ^P3/2  state  by 
a  dynamic  Jahn-Teller  effect.  Absorption  spectra  for  matrix  isolated  Au  atoms 

indicate  that  the  SO  coupling  is  always  greater  than  in  the  gas  phase.  The  spectra  of 
matrix  isolated  alkali  metals,  which  also  have  a  dominant    P  *-  ~S  transition,  show 

similar  matrix  induced  changes  from  the  gas  phase  spectra. 

For  a  species  trapped  in  a  frozen  noble  gas  (the  most  commonly  used  matrix 
material),  the  temperature  is  normally  maintained  near  10K  to  prevent  softening  of  the 
matrix  which  could  lead  to  diffusion  of  the  sample.    At  this  temperature  the  sample 
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molecule  is  vibrationally  and  rotationally  cold  resulting  in  a  great  simplification  of  the 

observed  spectra.  Also,  since  the  trapped  species  are  generally  in  their  electronic  ground 

state,  conventional  spectroscopic  probes  yield  no  information  about  the  excited  states  of 

the  species.      Magnetic  circular  dichroism   (MCD)   provides  a  means   of  extracting 

information  about  excited  states  based  on  polarization  dependent  selection  rules.    MCD 

is  the  measurement  of  the  differential  absorption  of  left-  and  right-circularly  polarized 

(LCP  and  RCP,  respectively)  light  when  a  sample  is  place  in  a  magnetic  field  parallel  to 

the  propagation  direction  of  the  light.  The  theoretical  development  of  this  technique  was 

given  by  Buckingham  and  Stephens27  in  1966.   The  first  application  of  MCD  to  matrix 

^o 

isolated  species  was  reported  by  Douglas,  Grinter  and  Thomson     at  the  University  of 
East  Anglia  in  1973. 

Matrix  Isolation 

Matrix  isolation  involves  the  trapping  of  a  sample  species  within  a  large  excess 
of  frozen  host  gas  at  a  dilution  ratio  high  enough  to  prevent  (or  minimize)  reactions 
between  the  guest  molecules.  The  matrix  (host)  gas  is  usually  chosen  to  be  non-reactive 
(typically  neon,  argon,  krypton  or  xenon)  but  in  some  instances  where  a  reaction  is 
required  to  produce  the  desired  sample  or  wider  temperature  ranges  for  annealing  studies 
are  necessary  other  gases  (CH4,  C3H8,  SF6)  may  be  used.  Several  recent  texts29-32 
provide  more  detailed  descriptions  of  the  varied  techniques  and  procedures  employed  in 
matrix  isolation  than  the  brief  description  given  below. 
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The  choice  of  matrix  gas  determines  the  method  of  cooling  necessary  to  freeze  the 

gas  on  the  sample  substrate.   To  prevent  diffusion  of  the  sample  through  the  matrix  the 

temperature  is  generally  maintained  lower  than  40%  of  the  melting  point  of  the  matrix 

gas.   Liquid  N2  (bp  77K)  and  liquid  He  (bp  4.2K)  are  commonly  used  cryogenic  liquids, 

but  only  liquid  He  provides  enough  cooling  to  form  rigid  noble  gas  matrices.  Liquid  H2 

will  freeze  all  but  neon  but  has  been  abandoned  due  to  the  fire  risk  involved.    The  use 

of   liquid    cryogens    requires   specially   constructed   dewars    and    appropriate   safety 

precautions.     The  introduction  of  closed  cycle  refrigeration  units  employing  He  has 

spurred  many  applications  of  matrix  isolation.  These  units  are  capable  of  reaching  ~10K 

which  allows  them  to  be  used  with  Ar,  Kr  and  Xe  matrices.   The  Heliplex,  marketed  by 

APD  Cryogenics,  uses  a  standard  refrigeration  unit  coupled  with  a  Joule-Thompson 

expansion  to  allow  cooling  to  4K,  enabling  work  with  neon  matrices.  These  systems  are 

typically  equipped  with  thermocouples  for  monitoring  temperature  of  the  substrate  and 

a  resistance  heater  for  annealing  studies  or  for  stabilizing  the  matrix  at  a  temperature 

higher  than  the  lowest  attainable  by  the  refrigerator.   Requiring  little  maintanence  and  no 

handling  of  cryogenic  liquids,  these  units  have  become  the  mainstay  of  many  matrix 

isolation  experiments. 

The  choice  of  cold  substrate  also  plays  an  important  role  in  the  matrix  isolation 

technique.  Two  primary  considerations  are  thermal  conductivity  and  transparency  of  the 

substrate  to  the  probing  radiation.    For  transmission  techniques,  a  window  is  usually 

mounted  in  a  copper  holder  which  is  attached  to  the  cold  finger  of  the  refrigeration  unit. 

Alkali  metal  halides  are  used  in  the  infrared  region  while  quartz,  CaF2  or  sapphire  are 
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employed  in  the  UV-visible  regions.    The  windows  are  mounted  in  the  holder  using 

indium  seals  to  insure  good  thermal  transfer  from  the  copper  to  the  window.    In  Raman 

spectroscopy,  the  cold  surface  is  often  a  polished  Cu  or  Al  block  attached  directly  to  the 

cold   finger.      Single   crystal  sapphire  or  Cu  rods  are  routinely  employed   in  ESR 

investigations. 

The  deposition  system  and  cold  surface  must  be  maintained  in  a  vacuum  (typically 

<  10~5  torr)  to  prevent  condensation  of  atmospheric  gases  onto  the  cold  surface  and  to 

prevent  possible  reactions  of  the  sample  with  the  gases.    For  gaseous  samples  or  those 

with  a  high  vapor  pressure  at  room  temperature  premixing  of  the  sample  and  matrix  gas 

may  be  done  prior  to  deposition.   Host:guest  ratios  of  103:1  to  10  :1  are  typically  used. 

For  samples  with  lower  vapor  pressures,  heating  may  be  employed  to  produce  a  molecular 

or  atomic  beam  which  may  then  mix  with  a  continuous  flow  of  matrix  gas  before 

condensing  on  a  sample  window.  Both  resistive  and  inductive  heating  methods  have  been 

employed.  The  availability  of  high  power,  low  cost  lasers  has  led  to  a  new  technique  for 

vaporization  of  refractory  materials.  The  output  of  a  laser  (often  a  Q-switched  Nd:YAG 

operating  at  532  nm)  is  focused  onto  a  solid  target  within  the  vacuum  chamber.    The 

ablated  sample  is  then  mixed  with  either  a  continuous  flow  of  matrix  gas  or  with  a  pulse 

of  gas  timed  to  coincide  with  the  firing  of  the  laser.    Deposition  times  may  range  from 

minutes  to  hours  depending  on  the  volatility  of  the  sample  and  the  concentration  of 

sample  necessary  for  spectroscopic  study. 


Magnetic  Circular  Dichroism 

Circular  dichroism  results  when  the  absorption  coefficients  for  RCP  and  LCP  light 
differ  for  a  medium.  This  inequality  may  occur  naturally  as  a  result  of  low  symmetry  for 
a  molecule  or  unit  cell  in  the  case  of  crystals.  Placing  the  medium  within  a  magnetic 
field  parallel  to  the  direction  of  light  propagation  will  also  produce  a  difference  in  the 
absorption  coefficients.  The  differential  absorption  of  RCP  and  LCP  light  in  the  presence 
of  the  magnetic  field  is  termed  magnetic  circular  dichroism.  This  results  in  a  signed 
spectrum  (both  positive  and  negative  bands  are  observed)  that  is  governed  by  polarization 
dependent  selection  rules.  The  observed  transitions  can  provide  ground  state  information 
about  atoms  and  molecules  but  are  particularly  useful  in  extracting  excited  state 
information  that  is  otherwise  unavailable. 

In  practice,  the  sample  is  placed  in  a  magnetic  field  colinear  with  the  propagation 
direction  of  alternating  RCP  and  LCP  light.  The  differential  absorption  is  detected  by  a 
photomultiplier  tube  whose  output  is  connected  to  a  phase  sensitive  detector.  The  MCD 
spectrum  is  then  displayed  as  a  wavelength  function  of  differential  absorption.  Additional 
requirements  for  successful  measurement  of  MCD  spectra  are  that  all  optical  components 
are  non-birefringent  and  that  the  sample  does  not  cause  depolarization  of  the  circularly 
polarized  light  by  scattering. 

The  MCD  spectrum  is  actually  the  difference  between  the  RCP  and  LCP 
longitudinal  Zeeman  spectra  and  the  same  information  is  available  from  the  Zeeman 
spectra  if  all  Zeeman  components  are  resolved.  Resolution  of  these  components  is  often 
difficult  and  requires  high  magnetic  fields.    MCD  requires  only  modest  fields  (-5000 
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gauss  in  the  present  study)  and  has  increased  sensitivity  due  to  phase  sensitive  detection 

techniques.   The  signed  nature  of  the  MCD  signal  may  allow  unresolved  absorptions  to 

be  discovered  if  the  MCD  selection  rules  dictate  oppositely  signed  signals  for  the 

transitions.   The  temperature  dependence  of  some  MCD  signals  also  provides  additional 

information  for  molecules  with  degenerate  ground  states. 

Aims  of  this  Investigation 

The  matrix-induced  changes  in  the  SO  coupling  of  matrix  isolated  Ag  and  Au 
atoms  and  the  origin  of  the  tri-peaked  structure  of  the  spectra  will  be  investigated  by 
obtaining  the  MCD  and  absorption  spectra  and  subjecting  the  data  to  a  band  moment 
analysis.  The  analysis  will  be  used  to  quantitatively  extract  SO  coupling  constants  and 
the  relative  contributions  of  totally  symmetric  (cubic)  and  Jahn-Teller  active  nontotally 
symmetric  (noncubic)  modes  to  the  spectral  bandwidth.  The  derived  SO  coupling 
constants  will  be  used  to  confirm  a  recent  model  proposed  by  Pellow  and  Vala33  (PV)  to 
explain  the  wide  variation  in  the  SO  splitting  of  the  2P  *-  2S  transition  of  the  alkali  and 
noble  metals  in  matrices.   The  origin  of  the  three  maxima  will  be  shown  to  result  from 

7  "7 

a  dynamic  Jahn-Teller  effect.  A  similar  analysis  will  be  performed  for  the  P  *-  S 
transition  in  matrix  isolated  Mo  atoms.  The  results  will  be  compared  with  a  theoretical 
prediction  derived  from  an  extension  of  the  PV  theory. 

The  MCD  and  absorption  spectra  of  Co  atoms  and  Ti  atoms  in  several  matrices 
will  also  be  presented.  Electronic  transitions  will  be  assigned  via  a  comparison  to 
simulations  based  on  gas  phase  data.   Analysis  of  the  MCD  bandshapes  will  be  used  to 
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diagnose  several  previously  misassigned  transitions.     The  study  of  transition  metal 

diatomics  using  the  MCD  technique  will  be  suggested  as  future  work  and  modifications 

to  the  apparatus  to  produce  these  species  will  be  outlined. 


CHAPTER  II 
MAGNETIC  CIRCULAR  DICHROISM  THEORY 


Introduction 

The  theoretical  basis  for  the  MCD  phenomenon  and  the  equations  necessary  to 
extract  useful  information  from  experimental  data  will  be  developed  in  this  chapter.  An 
outline  of  Stephens'  34,35  derivation  of  the  parametric  MCD  equations  will  be  presented 
to  detail  the  approximations  and  limitations  of  the  method.  The  extraction  of  excited  state 
information  will  employ  the  method  of  moments  originally  pioneered  by  Henry, 
Schnatterly  and  Slichter.36  Osborne  and  Stephens37  have  recast  the  theory  in  more 
general  terms  and  applied  it  to  MCD.  The  formalism  of  their  derivation  will  be  used  in 
this  work. 

Circularly  Polarized  Light 

A  description  of  circularly  polarized  light  will  be  given  to  aid  in  understanding  its 
interactions  with  a  sample.  A  more  thorough  description  may  be  found  in  reference  38. 
Unpolarized  light  consists  of  a  electric  vector,  E  ,  and  a  magnetic  vector,  B  ,  which  vary 
sinusoidally  in  the  direction  of  propagation  but  with  random  direction.  In  plane  polarized 
light  the  vectors  still  vary  sinusoidally  in  the  direction  of  propagation  but  their  direction 
is  constant.  If  this  plane  polarized  light  passes  through  a  transparent  birefringent  material 
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where  the  speed  of  light  depends  on  the  direction  of  polarization  then  one  part  of  the 

beam  will  be  retarded  with  respect  to  the  other.    By  passing  the  plane  polarized  beam 

through  the  medium  at  45°  to  the  perpendicular  axes  of  the  differing  refractive  indices 

projections  of  the  E   and  B   vectors  will  show  components  of  equal  magnitude  phase 

shifted  relative  to  one  another.    The  resultant  electromagnetic  vectors  will  trace  helical 

paths  as  they  emerge  from  the  medium.  A  clockwise  tracing  of  the  helix  as  an  observer 

looks  toward  the  light  source  is  defined  as  right  elliptical  polarization.    For  elliptical 

polarization  the  E    and  B    vectors  of  the  emergent  light  vary  both  in  direction  and 

magnitude.  Two  special  cases  can  occur  in  this  process.   If  the  phase  shift  is  equal  to  "k 

the  resultant  beam  is  returned  to  plane  polarization.   If  the  phase  shift  is  X/4  the  light  is 

termed  circularly  polarized  because  the  cross  section  of  the  helical  trace  is  a  circle  (Fig 

2.1).   In  this  case  the  E  and  B  vectors  vary  in  direction  but  are  constant  in  magnitude. 

In  MCD  the  phase  shift  must  be  carefully  controlled  to  provide  circular  polarization 

throughout  the  spectral  region  of  interest. 

If  left  (LCP)  and  right  circularly  polarized  (RCP)  light  are  absorbed  in  differing 

amounts  by  a  substance  the  material  exhibits  circular  dichroism.    If  a  magnetic  field  is 

present  along  the  axis  of  light  propagation  all  materials  within  the  magnetic  field  will 

exhibit  a  circular  dichroism  termed  magnetic  circular  dichorism  (MCD). 

Basic  MCD  Equations 

Consider  a  sample  with  ground  state,  A,  and  excited  states,  J  and  K,  placed  in  a 
magnetic  field  that  is  irradiated  by  circularly  polarized  light  collinear  with  the  magnetic 


Figure  2.1  Resolution  of  circularly  polarized  light  into  electric  and  magnetic  field 
vectors  which  vary  sinusoidally  but  with  one  retarded  by  X/4.  A  line 
connecting  the  resultant  vectors  will  trace  a  helical  path  of  circular  cross- 
section,  (ref.  38) 
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field.  This  situation  is  represented  schematically  in  Figure  2.2.  The  intensity  of  the  light 

at  any  point  z  is  given  by  the  Poynting  vector: 


I±(z)  =  -^n±(E0u)2  exp 


-2£k^z 


4ji   *'   u  fcc 

(1) 

0  -2£k±z 

=  I     exp  

±     P       *c 


where  k±  represents  the  absorption  coefficients  of  right(+)  and  left(-)  circularly  polarized 
light,  £=hv  is  the  photon  energy,  t>  is  Planck's  constant  divided  by  2jt,  n±  are  the  complex 
refractive  indices  for  LCP  and  RCP  light,  Eg  is  the  electric  field  at  z=0  and  I±  is  the 
incident  light  intensity. 

The  amount  of  energy  absorbed  per  unit  time  at  z  can  be  expressed  as  -dl(z)/dz 
and  is  dependent  on  the  photon  energy  £,  the  number  of  sample  species,  Na,  and  the 
transition  probability,  Pa_;. 


^=£NaPa-^)£  (2) 

dz  a,j 


The  above  probabilities  can  be  related  to  the  electric  dipole  transition  moments 
using  time  dependent  perturbation  theory  if  the  effective  Hamiltonian  ^0  is  assumed  to 
be  a  sum  of  independent  components.   Considering  only  electric  dipole  transitions 
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%  =  K +  *b'  (3) 


where  W^  is  the  Hamiltonian  in  the  absence  of  light  and  H^  is  the  Hamiltonian  for  the 
perturbation  induced  by  the  light.  At  this  point  the  external  magnetic  field  is  not  yet 
considered.   This  is  denoted  by  the  subscript  0.   This  approach  yields 


Pa_j(z)  =  ^  |E±(z)|2  |<a  |  m±  |  j>|26(£ja-£)  (4) 


where  <a  |  m±  |  j>  are  the  transition  moments,  m±  is  the  electric  dipole  transition  moment 
operator  for  RCP  and  LCP  light,  T.-   =  %■  -  £a,  6  is  the  Kronecker  delta  operator  arising 

from  the  neglect  of  transition  lifetimes  and  a  is  the  Lorentz  effective  field  correction 
factor  which  relates  the  electric  field  due  to  the  light  at  the  absorption  center  to  the 
macroscopic  electric  field. 

By  differentiating  (1)  and  employing  (2)  and  (4)  the  absorption  coefficients  can 
now  be  expressed  as 


k        2^V£  |<a|m    |jH28(E      £)  (5) 

V±       a,j 


Transforming  (5)  to  absorbance  units  gives 
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—  =  —  cz  =  Y±E  ^  l<a  I  m±  I  j>|26(£ja-'£)cz  (6) 


where 


_  N0Ji^azlog10e 
y±  "       250^ctT 


(7) 


Na/N  is  the  relative  population  of  the  ground  state  A,  N0  is  Avogadro's  number,  c  is  the 
speed  of  light,  E±  is  the  molar  absorption  coefficient  and  c  is  the  concentration  of 
absorbers.  This  result  assumes  that  Beer's  law  applies  and  is  not  valid  for  non-linear 
effects  produced  by  intense  light.   The  MCD  is  expressed  as 

AA  =  AL-AR   =A_-A+  (8) 

which  (using  (6))  gives 


AA  N   r 

-  =  V±E  -^   l<a  I  m.  |  j>|2  -  |<a  |  m+  |  j>|2   6(£ia-£)cz        (9) 


£  73   N 


This  equation  (9)  is  general  and  applies  to  natural  CD  of  optically  active  materials.  MCD 
is  not  a  perturbation  of  natural  optical  activity  but  an  additive  effect  resulting  from 
differences  in  polarizabilities  in  the  presence  of  a  magnetic  field. 

To  reflect  the  actual  bandwidths  the  rigid  shift  model  is  used.   This  assumes  that 
the  zero  field  absorption  (ZFA)  band  shifts  along  the  energy  axis  upon  introduction  of  the 
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magnetic  field   but  does  not  change  shape.      The   model  requires  that  the   Born- 

Oppenheimer  approximation  hold  for  both  the  ground  and  excited  states.   Separating  the 

electronic  wavefunctions  into  nuclear  and  electronic  terms  yields 

|Aaa>  =vpA(r)R)Xa(R)  a=ltodA  (10) 


|J)j>  =WJx(r,R)Xj(R)  X  =  ltodj  (11) 

where  r  is  the  electronic  coordinate,  R  is  the  nuclear  coordinate,  Xa  and  X:  are  the 
vibrational  wavefunctions  and  da  and  i  are  the  degeneracies  of  the  ground  and  excited 
states.  To  further  simplify  the  matrix  elements  the  Franck-Condon  approximation  that 
most  electronic  transitions  occur  at  R  near  the  equilibrium  nuclear  separation  R0  is 
embraced. 

Applying  the  rigid  shift  model  and  summing  over  all  vibronic  bands  leads  to  an 
equation  for  the  ZFA 


A0 

_^  =Y±2>0f(£)Cz  (12) 


where  z  is  the  pathlength,  c  is  the  concentration  and 
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=  ~i-  E  I  <A«  I  m.  |  J,>°  |2  +  |  <Aa  |  m+  |  J,>°  |2  (13) 


2dA  a,x 


Vibrational  overlaps  are  explicit  in  the  lineshape  function  f(£). 


f(£)  ■  E 

a.j 


<a|j>|26(2;ja-£)  (I4) 


Jo"f(E)dE-l  (15) 


Equations  (14)  and  (15)  show  that  the  lineshape  is  temperature  dependent  due  to  the 
Boltzmann  term  but  that  the  integrated  intensity  is  temperature  independent. 

The  rigid  shift  model  is  now  applied  to  molecules  in  a  magnetic  field. 
Considering  only  electronic  terms  the  perturbation  of  the  system  Hamiltonian  may  be 
expressed  as 

H'  S  -u^B  =  nB(Lz  +  2SZ)B  (16) 

where  \iz  is  the  electronic  magnetic  moment  along  z  and  \iB  is  the  electronic  Bohr 
magneton.  The  orbital  and  spin  angular  momenta  along  z  are  represented  by  L^  and  Sz, 
respectively.  Equation  (16)  describes  the  Zeeman  splitting  induced  by  the  magnetic  field. 
If  the  electronic  wavefunctions  are  chosen  to  diagonalize  \iz,  then  9i  '  will  be  diagonal 
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in  the  Franck-Condon  approximation.  Stephens35  now  explicitly  includes  intermixing  of 

additional  states  K  to  yield  a  better  approximation  of  the  transition  moments. 

The  ground  state  Zeeman  splitting  leads  to  population  changes.  In  the  case  where 

the  Zeeman  splitting  is  small  compared  to  kT,  the  Boltzmann  term  can  be  expressed  as 

a  Taylor  expansion  which  yields 


N„        1 


£  -^{l  -<AjLz+2Sz|Aa>nBB/kT} 


(17) 


Combining  the  assumptions  of  the  rigid  shift  model  and  applying  them  to  (9) 
produces  the  parametric  equation  for  the  MCD  of  the  A  -*  J  transition 
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The  terms  Av  CB0  and  C0  are  given  by 
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(21) 


where  Wa,  Wj  and  Wk  are  the  energies  of  the  states  and  the  other  terms  are  as  previously 
defined. 

Experimental  forms  of  the  Av  $0  and  C0  terms  are  shown  in  Fig.  2.3.  The  Ax 
term  arises  when  either  the  ground  or  excited  states  are  degenerate.  The  two  lobes,  which 
arise  from  the  small  Zeeman  splitting  of  the  degenerate  level,  are  always  oppositely 
signed.  The  Al  term  intensity  is  influenced  by  the  sharpness  of  the  absorption  band, 
following  its  dependence  on  the  derivative  of  the  lineshape  function,  and  by  the  difference 
in  magnetic  moments  of  the  ground  and  excited  states.   It  the  moments  are  equal  no  Ax 
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term  will  appear.  A  positive  Aj  term  is  defined  when  the  high  energy  lobe  is  of  positive 

sign. 

Degeneracy  of  the  ground  state  is  necessary  for  the  appearance  of  CQ  terms.  These 
depend  on  the  magnetic  moment  of  the  ground  state  and  on  the  circularly  polarized 
transition  moments.   A  C0  term  may  be  identified  by  its  1/T  intensity  dependence. 

A  'Bq  term  is  present  for  all  materials  as  a  result  of  magnetic  field  induced  mixing 
of  electronic  wavef unctions.  These  terms  are  particularly  difficult  to  assess  because  the 
energies  and  multiplicities  of  all  states  near  A  and  J  must  be  known.  Fortunately,  these 
are  also  the  least  intense  MCD  bands. 

The  MCD  spectrum  of  a  species  can  be  deconvoluted  into  contributions  from  each 
of  these  three  terms.  Also,  if  Sip  ®0  and  C0  can  be  calculated  the  MCD  spectrum  can  be 
predicted.  It  is  useful  to  determine  the  relative  contribution  of  each  term  to  the  overall 
MCD.   The  maximum  contributions  from  each  of  the  three  terms  is 

*!  :20:C0«^  :  _?_  :  JL  (22) 

1        u        u       T      AW      kT 

where  Z  is  the  Zeeman  energy,  T  is  the  bandwidth  of  the  electronic  transition,  AW  is  the 

order  of  magnitude  of  an  electronic  energy  gap,  k  is  the  Boltzmann  constant  and  T  is  the 

temperature.   For  a  matrix  isolated  sample,  the  assumptions  T  «  200cm"1,  kT  =  10cm"1 

and  AW  =  2000cm"1  may  be  used,  leading  to 

Aj  :  S0  :  C0  -  10  :  1  :  200  (23) 

which  indicates  that  CQ  terms  should  dominate  the  spectrum  if  they  are  present. 
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Figure  2.3        Appearance  of  MCD  terms. 

(a)  Zero  field  absorption  and  MCD  for  a  positive  !A^  term.  The 
positive  lobe  lies  at  higher  energy. 

(b)  #q  terms  may  be  positive  or  negative  with  the  maximum  coincident 
with  the  absorption  maximum. 

(c)  CQ  terms  may  be  positive  or  negative  with  the  maximum  coincident 
with  the  absorption  maximum  and  will  show  a  temperature 
dependence. 
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Moments 


(24) 


The  general  absorption  and  MCD  moments  are  defined  as 

<A>n  =J"(A/v)(v-v0)ndv 

<AA>n  =  [(AA/vXv^o)11  dv  (25) 

where  the  integration  is  over  the  entire  band  (i.e.,  ^P^  and  P3/2  components);  A  and  AA 
are  the  optical  absorbance  (A)  and  differential  (MCD)  absorbance  per  Tesla.  The  average 
energy  of  the  zero  field  absorption  is  given  by 

v0  =  fAdv/f(A/v)dv  (26) 

where  v  is  the  photon  frequency. 

The  "reduced"  moment  form  represents  the  division  of  the  nth  MCD  moment  by 
the  zeroth  absorption  moment  and  has  the  general  form 


<AA>n        f  Cn' 

<A>n         [  kTj 


HBB 


B D  (27) 


*>o 


^7 

where  An,  #n,  CQ  and  ©0  are  as  defined  in  Osborne  and  Stephens     and  the  other  terms 

are  as  previously  defined.  The  moment  equations  relevant  to  the  present  analysis, 
assuming  harmonic  vibrational  frequencies,  lowest  order  in  SO  coupling  and  first  order 
electron-lattice  interactions  are 
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<AA>1/<A>0  =  2^^  -  a£/kT)  (28) 

<A>2/<A>0  =  <E2>C  +  <E2>NC  +  as^2  (29) 

<AA>3/<A>0  =  6fxBB(gorb  -  a£/kT)(<E2>c  +  <E2>NC/2  +  a£2/2)      (30) 

<AA>0  =  0  (31) 


Here  as  is  defined39  as 


as  =  -(kT/nBB)(l/(2S+l))Emsexp(-msgnBB/kT)  (32) 


ms 


which  in  the  limit  kT  »  ms  g  |iBB  for  S  =  1/2  yields  as  =  1/2.  <E2>C  and  <E2>NC 
represent  the  cubic  and  noncubic  lattice  mode  contributions  to  the  observed  bandwidths. 
X,  is  the  SO  coupling  constant  of  the  excited  2P  state  (the  actual  splitting  between  2P1/2 
and    P3/2  multiplets  is  A  =  2>tJ2),  and  g^  is  the  excited  "P  state  orbital  g  factor. 

From  (28),  the  slope  of  the  plot  of  <AA>1/<A>0  vs.  1/T  yields  the  SO  coupling 
constant,  £,  which  upon  insertion  into  (29)  and  (30)  permits  the  determination  of  <E  >c 
and  <E  >NC.  It  is  expected  that  for  a  system  exhibiting  the  Jahn-Teller  (JT)  effect  in  the 
P  state,  the  value  of  <E  >NC  will  be  nonzero.  Acting  on  a  suggestion  by  Stephens, 
Schatz  and  coworkers24  have  shown  that  a  crystal  field-induced  splitting  (derived  from 
a  distorted  site  symmetry)  may  be  differentiated  from  JT  activity  by  the  temperature 
independence  of  its  <A>2/<A>0  plot.  For  an  active  JT  effect,  the  ratio  <A>2/<A>0  (and 
thus  the  <E  >NC  component)  should  show  a  temperature  dependence     according  to 


24 


<E2>NC  =  <E2>J^C  coth(Vj/2kT)  (33) 

where  <E  >NC    is  the  contribution  of  the  noncubic  lattice  vibration  modes  to  the 
bandwidth  with  T  approaching  OK.   Here 


<E2>NC  s  <eg>  +  <t2g>  (34) 


where  <e  >  or  <t2„>  are  given  by 


<Yg>  =  E(^i2/-^i)vi  cothCvj/ZkT)  (35) 


'g 

i 


in  which  £i  is  the  linear  electron-lattice  coupling  parameter,  k{  the  ground  state  force 
constant  and  v;  the  frequency  of  the  contributing  mode. 


CHAPTER  III 
APPARATUS  AND  EXPERIMENTAL  TECHNIQUES 


Introduction 

The  process  of  studying  matrix  isolated  species  using  MCD  is  conveniently 
divided  into  three  sections:  matrix  isolation  of  the  species  of  interest,  spectroscopic  study 
of  the  sample  and  calibration  of  the  absorbance,  MCD  and  temperature.  This  also 
provides  an  effective  means  of  detailing  the  apparatus  and  techniques  used  in  the  present 
work.  The  first  section  will  deal  with  the  matrix  isolation  assembly  and  matrix 
preparation  procedures,  the  second  will  describe  the  optics  and  electronics  of  the 
spectrometer,  while  the  third  will  discuss  calibrations.  Individual  experimental  details  will 
be  given  for  each  species  in  their  respective  chapters. 

Matrix  Isolation  Equipment 

To  prevent  the  condensation  of  atmospheric  contaminates  on  the  cold  sample 
window,  the  entire  sample  preparation  and  isolation  setup  must  be  maintained  under  high 
vacuum  conditions.  The  high  vacuum  system  employed  consisted  of  a  Alcatel  pumping 
station  which  included  a  roughing  pump  (12  m3/hr),  diffusion  pump  (150  L/sec),  liquid 
nitrogen  trap,  three-way  control  valve  and  internal  butterfly  valve.  The  pumping  station 
was  connected  to  the  furnace-cryostat  assembly  via  a  Edwards  rotary  high-vacuum  valve 
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and  a  length  of  corrogated  stainless-steel  tubing.      All  connections  were  made  with 

Kwik-Flange  connectors.  In  experiments  where  powdered  samples  were  vaporized,  rapid 

initial  pumping  could  result  in  the  powder  being  blown  out  of  the  sample  cell  by  the 

outrush  of  trapped  atmosphere.    The  rotary  valve  allowed  for  fine  control  of  rough 

pumping  speed.    System  pressure  was  monitored  by  thermocouple  (Alcatel  model  CA 

101)  and  Penning  ionization  (Alcatel  model  FA  101)  gauges,  both  of  which  were  mounted 

on  the  pumping  station  control  area. 

Sample  preparation  followed  one  of  two  distinct  routes.  For  metals  with  relatively 
low  melting  points,  vaporization  was  accomplished  by  resistive  heating.  A  small  stainless 
steel  furnace  (4  inch  diameter)  was  constucted  to  allow  the  combined  furnace/cryostat 
assembly  to  remain  inside  the  pole  faces  of  the  electromagnet  and  to  bring  the  source  of 
the  metal  beam  to  within  four  inches  of  the  sample  window.  The  furnace  was  fabricated 
from  a  single  piece  of  stainless-steel  which  was  then  covered  with  another  stainless-steel 
tube.  This  tube  was  heliarc  welded  to  the  outside  of  the  furnace  can  and  fitted  with  inlet 
and  outlet  connections  to  allow  water  cooling  of  the  furnace  body  and  the  internal  heat 
shield.  This  system  is  shown  in  figure  3.1.  This  system  corrected  an  earlier  problem  of 
cooling  water  leaks  into  the  vacuum  space  as  all  welds  were  external. 

Metal  powders  were  placed  into  a  one-inch-long  section  of  tantalum  tubing 
(0.015"  wall  thickness)  with  a  small  hole  drilled  midway  along  the  body.  The  ends  were 
sealed  with  solid  tantalum  endplugs  and  fitted  with  angled  tantalum  brackets.  This 
Knudsen  cell  was  secured  to  water  cooled  copper  electrodes  mounted  on  a  stainless  steel 
flange  that  closed  the  back  end  of  the  furnace  (also  shown  in  Fig  3.1).   The  exit  hole  of 
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the  Knudsen  cell  was  carefully  positioned  to  be  colinear  with  the  center  of  the  sample 

window.    One  of  the  electrodes  was  electrically  isolated  by  using  a  teflon  washer  and 

teflon  screws  for  mounting. 

Resistive  heating  of  the  Knudsen  cell  was  accomplished  by  passing  an  AC  current 
of  up  to  300  Amps,  60  Hz  through  the  electrodes.  The  temperature  of  the  Knudsen  cell 
was  monitored  by  a  IRCON  digital  pyrometer  (model  220)  thru  the  pyrometer  veiwport 
in  the  electrode  flange.  These  samples  generally  required  outgassing  to  drive  off  adsorbed 
H20  and  any  contaminants  more  volatile  than  the  sample.  The  displex  head  was  rotated 
to  protect  the  sample  window  from  any  vaporized  contaminants  and  the  cell  was  heated 
to  approximately  3/4  of  the  expected  deposition  temperature.  This  temperature  was 
maintained  for  20-30  minutes.  For  depositions  the  sample  cell  was  typically  heated  to 
a  temperature  that  resulted  in  a  vapor  pressure  of  10-3  torr  for  the  sample.  The  viewport 
was  equipped  with  a  magnetic  shutter  which  was  opened  only  to  read  the  temperature. 
This  procedure  prevented  metal  from  being  deposited  on  the  viewport  window  which 
would  cause  a  reduced  temperature  reading. 

Highly  refractory  samples  were  more  easily  vaporized  by  direct  laser  ablation 
using  the  second  harmonic  (532  nm)  of  a  pulsed  Nd:YAG  laser  (Spectra  Physics  model 
DCR-11).  A  variety  of  sample  forms  were  available  including  rods,  ingots,  wire  or  discs. 
These  discs  were  made  by  compressing  a  quantity  of  metal  powder  in  a  stainless  steel  die 
using  a  hydraulic  press.  In  most  cases,  the  sample  was  rotated  using  a  small  electric 
motor  to  continuously  provide  a  clean  surface  and  to  prevent  drilling  through  the  sample 
with  the  laser.  Since  resistive  heating  was  not  required  a  simpler  sample  holder  was  built. 
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A  small  brass  block  that  could  be  attached  directly  to  the  cryostat  vacuum  shroud  was 

bored   in   three   directions   to   allow   sample    insertion,   matrix    gas   and    laser   beam 

introduction.   This  method  of  sample  vaporization  is  diagramed  in  Fig.  3.2. 

For  resistively  heated  samples,  initial  system  pump  down  was  controlled  using  the 
rotary  valve  on  the  pump  stack.  Slow  rough  pumping  of  the  system  insured  that  the 
powder  was  not  blown  out  of  the  Knudsen  cell.  Approximately  25-40  minutes  were 
required  to  fully  open  the  roughing  valve  after  filling  (or  refilling)  a  Knudsen  cell.  For 
laser  vaporized  samples,  more  rapid  roughing  of  the  system  was  possible.  After  the 
pressure  in  the  system  reached  ~2xl0~2  mbar,  the  three  way  valve  was  set  to  diffusion 
pump  and  the  internal  butterfly  valve  was  opened.  For  powdered  samples  the  pressure 
was  allowed  to  decrease  to  the  10  mbar  range  before  outgassing  of  the  cell  was 
performed.  Generally,  pumping  of  the  system  was  allowed  to  continue  overnight.  With 
this  arrangement,  ultimate  system  pressures  typically  were  recorded  as  2x10  mbar  with 
the  liquid  N2  trap  filled  and  the  cryostat  off.  With  the  cryostat  cold  the  pressure 
generally  read  5x10     mbar. 

The  cold  surface  for  sample  condensation  was  either  a  CaF2  or  sapphire  window. 
The  sapphire  window  was  used  whenever  possible  due  to  its  greater  thermal  conductivity. 
Only  sapphire  windows  cut  on  the  C  axis  of  a  single  crystal  were  used  to  avoid 
birefringent  signals  from  the  window  overlapping  the  sample  signals.  The  window  was 
held  in  a  holder  made  from  oxygen-free  high  conductive  (OFHC)  copper.  The  windows 
used  were  2.54  cm  in  diameter  and  3  mm  thick.  The  window  holder  was  designed  to 
allow  optical  access  through  a  1.25  cm  diameter  area  at  the  center  of  the  window.   The 
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Figure  3.2        Cutaway  side  view  of  cryostat  with  attached  laser  vaporization  source. 
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remaining  area  of  the  window  was  covered  by  the  copper  holder  to  insure  maximum 

cooling.    Indium  gaskets  were  used  to  insure  good  contact  between  the  window  and  the 

holder.  The  holder  was  screwed  into  the  copper  cold  tip  of  either  an  Air  Products  (now 

APD  Cryogenics)  Displex  or  Heliplex.   Indium  gaskets  were  again  used  to  provide  good 

thermal  contact  between  the  cold  tip  and  the  window  holder.  The  Displex  is  a  two  stage 

closed-cycle  He  refigerator  capable  of  cooling  to  10  K  while  the  Heliplex  is  a  three  stage 

refrigeration  unit  capable  of  cooling  to  4  K  by  use  of  the  Joule-Thompson  expansion  in 

conjunction  with  a  standard  refrigeration  cycle.  An  internal  copper  heat  shield  surrounded 

the  cold  finger  and  window  holder.  This  shield  was  maintained  at  approximately  60-80K 

to  minimize  radiant  heat  transfer  to  the  cold  window.  A  thermocouple  (Au(0.07%Fe)  vs. 

Chromel)  was  mounted  in  a  small  hole  drilled  in  the  copper  window  holder  using  Crycon 

grease.     The  temperature  was  monitored  using  an  APD-E  control  unit  with  digital 

readout.   This  control  unit  also  provided  power  for  the  30  watt  resistive  heater  used  for 

controlling  the  temperature  of  the  window  during  measurement  of  temperature  dependent 

spectra. 

The  head  of  the  refrigeration  unit  was  rotatable  on  O-ring  seals  within  the  vacuum 

shroud.  This  allowed  sample  preparation  at  90°  to  the  optical  path.   The  shroud  had  five 

openings  to  allow  access  to  the  cold  finger.   Three  of  these  were  covered  with  Suprasil 

windows  to  allow  access  to  the  optical  path  and  for  viewing  or  photolysis  of  the  matrix. 

A  Kwik-flange  adapter  was  mounted  over  the  fourth  opening  to  allow  connection  of  the 

furnace.  The  fifth  access  in  the  bottom  of  the  shroud  was  blanked  off  during  experiments 

with  the  furnace  but  was  used  for  pumping  when  the  laser  vaporization  source  was  used. 
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Since  no  gate  valves  were  used  between  the  sample  source  and  the  shroud,  samples  could 

not  be  changed  without  bringing  the  cryostat  to  room  temperature  and  breaking  vacuum. 

A  diagram  of  the  Displex  is  given  in  Fig.  3.3. 

The  high  purity  (Matheson  Research  Grade,  99.9995%)  matrix  gas  was  introduced 

from  a  stainless  steel  manifold  designed  to  have  a  very  small  hold  volume.   A  standard 

two  stage  regulator  was  connected  to  an  analogue  pressure  meter,  a  pumping  port  and  a 

needle  valve  using  1/8"  O.D.  stainless  steel  tubing.  A  small  rotary  vacuum  valve  allowed 

the  pumping  port  to  be  closed  until  a  gas  change  was  necessary.   The  largest  contributor 

to  the  total  volume  was  the  pressure  gauge.  A  coil  of  1/8"  O.D.  stainless  steel  tubing  on 

the  outlet  side  of  the  needle  valve  was  used  for  trapping  any  residual  water  or  other 

contaminants  in  the  matrix  gas.    Typically  this  coil  was  immersed  in  a  Dewar  flask 

containing  a  slush  of  acetone/  liquid  nitrogen.  The  setup  is  shown  in  Fig.  3.4.  Flow  rates 

of  0.25  to  0.33  torr  per  minute  were  used  to  produce  the  matrices  and  were  easily 

controlled  by  adjustment  of  the  needle  valve. 

Spectroscopic  Apparatus 

Previous  attempts  to  study  the  MCD  of  Ag  atoms  in  matrices40  have  shown  that 
the  quantity  of  Ag  atoms  decreases  from  one  scan  of  the  spectrum  to  the  next.  From 
Ozin's  work41-42  this  is  now  known  to  be  the  result  of  cryophotoclustering.  Since  ratios 
of  MCD  and  absorption  moments  are  required  for  the  moment  analysis  of  the  bands  a 
new  instument  was  built  to  allow  the  simultaneous  recording  of  MCD  and  absorption 
spectra.   With  this  approach  the  continual  decrease  in  signal  intensity  did  not  affect  the 
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Figure  3.3        Diagram  of  closed-cycle  helium  refrigerator  (Displex)  expander  section. 
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analysis  since  ratios  of  MCD  and  absorption  moments  were  available  for  each  scan  of  the 

spectrum.    This  section  describes  the  optics,  electronics  and  computer  control  of  this 

psuedo-double  beam  spectrometer. 

The  optical  system  is  represented  schematically  in  Fig.  3.5.  The  output  of  a  300 
watt  xenon  lamp  fitted  with  an  internal  parabolic  reflector  and  sapphire  window  (Eimac) 
was  directed  onto  the  entrance  slit  of  a  3/4  meter  Czerny-Turner  monochromator  (Spex). 
The  wavelength  scan  of  the  monochromator  was  driven  by  a  200  step/revolution  stepper 
motor.  Pulses  to  drive  the  motor  were  produced  on  a  digital  I/O  line  of  a  Metrabyte  Dash 
16  D/A  interface  board  controlled  by  a  IBM  PC  clone  (80286  processor).  A  gaseous  N2 
purge  of  the  monochromator  housing  was  maintained  during  operation  to  prevent 
oxidation  of  atmospheric  contaminants  by  the  UV  output  of  the  lamp  and  subsequent 
degradation  of  the  mirrors  and  grating. 

Circular  polarization  was  provided  by  focusing  the  output  of  the  monochromator 
onto  a  Glan-Thomson  prism  oriented  a  45°  to  the  modulation  axis  of  the  photoelastic 
modulator  (Morvue  PEM-3).  To  maintain  circular  polarization  throughout  the  scanned 
wavelength  region  the  retardation  voltage  applied  to  the  piezoelectric  oscillator  was 
changed.  The  spectrometer  control  program  (see  appendix)  calculated  a  digital  value  for 
each  wavelength  which  was  then  applied  to  the  remote  modulation  control  input  of  the 
PEM  as  an  analog  voltage  from  D/A  channel  0  of  the  Dash- 16  board. 

The  emerging  beam  was  then  chopped  to  provide  reference  and  sample  beams. 
A  second  lens  allowed  focusing  of  the  sample  beam  at  the  center  of  the  gap  between  the 
magnet  pole  faces.   An  axially  bored  magnet  (Alpha  Scientific  model  4600)  allowed  the 
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beam  to  be  collinear  with  the  magnetic  field.  The  maximum  attainable  field  with  the  pole 

pieces  adjusted  to  allow  the  displex  shroud  to  slip  between  them  was  approximately  4750 

gauss  for  initial  experiments.  Later  experiments  utilized  a  newer  magnet  (Alpha  Scientific 

model  4800)  capable  of  reaching  approximately  5500  gauss  for  the  same  pole  gap.  After 

passing  through  the  sample  window  at  the  center  of  the  pole  gap  the  beam  passed  through 

the  second  pole  and  impinged  upon  the  photomultiplier  tube  (PMT)  (EMI  9973QB).  The 

reference  beam  was  directed  around  the  magnet  using  mirrors  and  focused  onto  the  same 

PMT.    All  mirrors  were  front  surfaced  to  allow  measurement  to  near  200  nm. 

A  specially  designed  optical  chopper  wheel,  shown  in  Fig  3.6,  provided  both 
pseudo-double  beam  operation  and  timing  signals  for  spectrometer  control  and  data 
collection.  The  mirrored  sections  directed  light  around  the  magnet  for  the  reference 
beam,  the  blackened  areas  prevented  light  from  the  monochromator  from  reaching  the 
PMT  and  provided  a  zero  signal  while  the  large  open  areas  allowed  light  to  pass  through 
the  wheel  and  strike  the  sample  matrix.  During  the  open  periods  both  sample  absorption 
and  MCD  were  recorded. 

The  square  notches  in  the  outer  edge  of  the  wheel  triggered  opto-interrupters 
which  provided  control  signals  to  switch  the  spectrometer  interface  between  reading 
sample,  reference  and  blank  signals.  A  second  opto-interrupter  was  triggered  by  the 
inner  absorption/MCD  slot  and  caused  the  interface  to  begin  recording  MCD.  The  signal 
timings  are  shown  schematically  in  Fig  3.6. 

The  trigger  signals  were  routed  to  the  gating  logic  board.  Each  signal  passed 
through  two  Schmidt  triggers  to  produce  a  clean  signal  for  a  monostable  multivibrator. 
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The  absorption  control  signals  passed  through  a  JK  flip-flop  network  arranged  as  a  divide 

by  three  counter.  The  three  output  signals  were  connected  to  three  separate  switch  control 
inputs  on  a  ADG-202  quad  switch  on  the  absorption  board.  The  MCD  control  signal 
passed  directly  from  the  multivibrator  to  the  signal  switch  on  the  absorption  board.  A 
clear  pulse  for  the  JK  flip-flop  network  was  also  derived  from  the  opening  of  the  MCD 
gate. 

The  absorption  and  MCD  signals  followed  two  distinct  paths  and  each  will  be 
described.  The  absorption  board  is  shown  in  Fig.  3.7.  The  voltage  input  from  the  PMT 
housing  was  amplified  and  fed  to  three  signal  input  positions  of  the  quad  switch.  The 
timing  signals  from  the  logic  board  closed  the  appropriate  switch  and  sample,  reference 
and  zero  signals  were  sent  to  their  respective  amplifiers.  The  zero  signal  was  then 
electronically  subtracted  from  the  sample  and  reference  signals.  These  zero  corrected 
sample  and  reference  signals  were  then  connected  to  the  inputs  of  a  logarithmic  amplifier 
(Burr-Brown  #4127)  whose  output,  the  absorbance  signal,  was  recorded  on  one  line  of 
a  dual  pen  strip  chart  recorder  (Soltec)  and  digitized  for  computer  storage  by  one  A/D 
channel  (12  bit  resolution)  of  the  Dash-16  board. 

To  prevent  the  MCD  signal  from  being  influenced  by  changes  in  the  single  beam 
absorption  of  the  sample  the  DC  portion  of  the  PMT  output,  which  contained  the 
absorption  information,  was  held  constant  by  a  feedback  circuit  driven  by  the  sample 
minus  zero  signal.  Holding  this  voltage  constant  insured  that  the  AC  portion  (50  kHz 
modulation)  of  the  PMT  signal  was  not  divided  by  the  DC  portion.  The  feedback  circuit 
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controlled  the  high  voltage  input  to  the  PMT  through  a  programmable  high  voltage  power 

supply  (Bertan  PMT-20A-N). 

For  MCD  measurement,  the  PMT  signal  was  not  amplified  but  fed  directly  to  the 
fourth  position  of  the  quad  switch.  During  the  MCD  measurement  window  the  output  of 
the  switch  was  routed  to  a  sharp  bandpass  filter  (50  kHz  ±  2  kHz)  and  then  to  an 
amplifier  circuit  fitted  with  a  variable  set  of  feedback  resistors  to  control  the  gain.  This 
signal  then  passed  directly  to  the  lock-in  amplifier  (Ithaco  model  353)  phase  locked  to 
the  50  kHz  reference  signal  from  the  PEM.  The  DC  output  of  the  lock-in  was  amplified 
and  recorded  on  the  second  pen  of  the  strip  chart  recorder.  The  second  D/A  channel  of 
the  Dash-16  board  was  used  to  digitize  the  MCD  signal  for  disk  storage. 

The  MCD  channel  was  often  found  to  be  extremely  noisy  which  resulted  in 
unusable  data.  This  problem  was  traced  to  ringing  induced  in  the  filter  circuit  when  the 
switch  for  MCD  closed.  The  instantaneous  application  of  a  voltage  caused  the  filter  to 
be  unstable  and  large  oscillations  to  appear  in  the  output  of  the  lock-in.  This  problem 
was  corrected  by  feeding  the  PMT  signal  directly  to  the  filter  circuit  and  from  there  to 
the  quad  switch.  Since  the  filter  was  now  active  at  all  times  no  oscillations  were  induced 
when  the  switch  closed.  After  passing  through  the  switch  the  signal  was  again  sent  to  the 
variable  gain  amplifier  and  then  on  to  the  lock-in.  The  present  setup  is  shown  in  Fig. 
3.8. 
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Calibration  Procedures 

Wavelength  calibrations  were  accomplished  by  passing  the  632.8  nm  line  of  a 
HeNe  laser  through  the  monochromator  and  adjusting  the  wavelength  readout  to  632.8  nm 
at  maximum  throughput.  This  procedure  was  performed  occassionally  throughout  the 
course  of  this  work  and  whenever  the  monochromator  was  moved.  It  was  also  done  when 
new  mirrors  and  the  cleaned  grating  were  installed  in  the  monochromator. 

Calibration  of  the  magnetic  field  was  done  with  a  portable  gaussmeter  (F.W.  Bell 
model  640).  The  probe  was  mounted  securely  in  the  center  of  the  pole  gap.  The  probe 
was  then  covered  with  a  mu  metal  shield  and  the  meter  zeroed.  The  mu  metal  was 
removed  and  current  was  applied  to  the  magnet.  This  was  done  whenever  the  pole  gap 
of  the  magnet  was  changed. 

Aqueous  solutions  of  (+)  -  Co(en)3Cl  •  tartrate  (en  =  ethylenediamine)  and  d10  - 
camphorsulphonic  acid  (CSA)  were  used  to  calibrate  the  absorption  and  MCD  spectra. 
The  absorbance  maximum  of  each  of  these  solutions  was  measured  using  either  a  Cary 
17D  or  Hewlett-Packard  diode  array  spectrometer.  The  spectra  of  these  solutions  was 
then  run  on  the  MCD  apparatus  along  with  a  sample  of  deionized  water  to  serve  as  a  zero 
absorption.  The  digital  value  of  the  water  baseline  at  X^^  was  then  subtracted  from  the 
digitized  value  of  the  sample  at  kmax  to  give  an  absorbance  value  in  digital  units.  The 
known  O.D.  from  the  Cary  or  Hewlett-Packard  runs  was  then  divided  by  this  digital 
value  to  give  a  calibration  factor  in  O.D./byte  which  was  then  used  by  the  moments 
program  (see  appendix)  to  convert  digital  data  units  to  O.D.  units. 
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The  circular  dichroism  of  each  of  these  standard  solutions  was  also  measured  and 

stored  by  the  computer.     Using  the  relation  AA^^/A^^  =  0.225  determined  by 

McCaffery  and  Mason43  for  (+)  -  Co(en)3Cl  ■  tartrate  the  AO.D.  was  determined  and 

divided  by  the  peak  height  (in  digital  units)  to  give  a  calibration  factor  in  AO.D./byte. 

Calibration  of  the  CD  using  CSA  was  more  involved.  The  concentration,  c,  of  the 

solution  was  determined  from 

O.D.  =  eel  (36) 

where  O.D.  is  the  measured  absorbance,  1  is  the  pathlength  and  e  is  the  molar  absorption 
coefficient.  The  molar  ellipticity  of  CSA  at  291  nm  is  7800  degree  cm2/decimole44  and 
is  related  to  the  observed  rotation,  8,  and  the  concentration  by 

[B]291  =  JL  x  100.  (37) 

A  value  for  8  was  determined  and  the  AO.D.  was  calculated  from 


8  =  33  x  AO.D.  (38) 

This  AO.D.  value  was  again  divided  by  the  zero  corrected  digitized  intensity  of  the  CSA 
circular  dichroism  to  give  a  calibration  factor  in  AO.D./byte. 

Uncertainty  in  absolute  temperature  of  the  sample  represents  a  large  source  of 
error  in  values  derived  from  plots  of  moment  ratios  vs.  1/T.  Sample  temperature  was 
monitored  with  a  thermocouple  (chromel/Au  (0.07%  Fe))  mounted  in  the  copper  window 
holder.  These  temperatures  were  corrected  using  a  procedure  reported  previously     which 
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uses    the    matrix-isolated   metal    itself   as   a   probe    of  the   temperature   of   its   own 

environment.  For  a  paramagnetic  species  the  temperature  dependence  of  the  moment  ratio 

may  be  expressed  as 

<^>0/<A>0  =(C0/'D0)(vLBB/kT)  (39) 

For  two  temperatures  (T2  >  Tj), 

T2  =  Tj   +  AT.  (40) 

Denoting  the  moment  ratio  as  %  and  using  the  ratio  XjIXv  li  can  ^e  shown  that  a  plot  of 
AT  vs.  {(Xi/X2)-1/  w^'  nave  a  slope  equal  to  the  inverse  of  the  lowest  temperature  (T^1) 
attainable.  The  interval  AT  is  taken  from  the  thermocouple  readings  and  the  moment  ratio 
from  the  lowest  energy  component  of  the  three  peaks  observed  for  both  Ag  and  Au 
containing  matrices.  This  procedure  was  performed  for  each  matrix  and  the  difference 
between  the  calculated  lowest  temperature  and  the  recorded  lowest  temperature  was  used 
as  a  correction  factor  for  all  temperatures.  The  Displex  thermocouple  reading  was 
generally  about  0.5K  too  low,  while  the  Heliplex  thermocouple  reading  was  consistently 
about  2.0K  too  low.  This  procedure  was  verified  in  the  following  manner.  A 
polyvinylalcohol  (PVA)  film  containing  K3Fe(CN)6  was  deposited  from  solution  directly 
onto  the  sample  window.  Absorption  and  MCD  spectra  were  recorded  over  the  range  10- 
40K.  After  temperature  correction  (using  the  above  procedure)  a  plot  of  AAo/(A0(j.BB) 
vs.  1/kT  for  the  400  nm  Fe(CN)63"  MCD  band  was  found  to  be  linear  with  a  slope  of 
0.84.  A  similar  study  by  Bryce  Williamson  at  the  University  of  Canterbury,  New  Zealand 
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using  a  free-standing  PVA  film  containing  K-,Fe(CN)6  in  a  immersion  cryostat  from  1.5 

to  300K  found  a  slope  of  0.86  for  the  same  plot.    In  this  case  the  film  was  surrounded 

either  by  liquid  He  (4.2K  and  below)  or  by  thermally  equilibrated  He  gas  (above  4.2K) 

ensuring  that  the  film  was  at  the  temperature  indicated  by  the  vapor  pressure  of  the  liquid 

He.    This  suggests  that  after  calibration  our  temperatures  may  still  be  -0.2K  too  low. 

This  uncertainty  will  translate  to  a  3-6%  uncertainty  in  the  SO  coupling  constants 

reported  later.  Additionally,  the  temperature  dependence  of  the  MCD  from  a  ground  state 

spin-doublet  such  as  possessed  by  Ag  or  Au  is  expected  to  vary  as  tanh((i.BB/kT)  which 

for  kT  »  fxBB  becomes  ~[xfiB/kT.   The  temperature  variation  at  a  fixed  magnetic  field 

followed  within  3.5%  what  would  be  expected  from  the  latter  limiting  expression.  These 

observations  provide  confidence  that  the  temperatures  used  in  the  moment  analyses  are 

accurate  to  a  reasonable  degree. 


CHAPTER  IV 

ABSORPTION  AND  MCD  OF  MATRIX  ISOLATED 

SILVER  AND  GOLD  ATOMS 


Introduction 

The  origin  of  the  three  absorption  maxima  observed  when  the  noble  metals  are 
isolated  in  rare  gas  matrices  has  stimulated  much  discussion.  Early  studies  '  favored 
a  distorted  site  model  where  an  axial  distortion  of  the  matrix  cage  atoms  surrounding  the 
metal  induces  a  static  "crystal  field"  splitting  of  the  2P  excited  state.  An  ESR  study  by 
Kasai  and  McLeod8  showed  that  Cu,  Ag  and  Au  atoms  sit  in  single,  substitutional, 
octahedral  (i.e.  nondistorted)  sites  in  the  matrix.  These  contradictory  conclusions  were 
reasonably  reconciled  by  suggestions1  9''  involving  simultaneous  SO  and  dynamic 
Jahn-Teller  (IT)  coupling  in  the  P  state.  These  suggestions  were  proven  for  Cu  atoms 
in  rare  gas  matrices  '  by  subjecting  temperature  dependent  MCD  and  absorption  data 
to  a  moments  analysis.  Noncubic  (Jahn-Teller  active)  lattice  modes  were  found  to  be  the 
major  contributors  to  the  observed  bandwidth  indicating  that  a  dynamic  JT  effect  was 
operating  within  the  matrix  environment.  The  effect  of  matrix  isolation  on  the  SO 
constants  was  also  assessed.  The  SO  constant  for  Cu  was  found  to  decrease  from  166 
cm-  in  the  gas  phase  to  124cm  in  Ar,  95  cm"1  in  Kr  and  reverse  sign  in  Xe,  to  -23 
cm-1.  Early  studies  of  Au  atoms  in  matrices10  indicated  from  studies  of  absorption 
spectra  that  the  SO  constant  increased  over  its  gas  phase  value     of  2543  cm"  . 
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Experimental 

Matrices  containing  Ag  and  Au  were  prepared  using  the  resistively  heated 
Knudsen  cell.  For  Ag,  powdered  metal  (Spex)  was  placed  in  the  Knudsen  cell  and  after 
outgassing  was  allowed  to  pump  overnight.  All  subsequent  deposits  were  made  without 
refilling  the  cell.  The  cell  was  heated  to  900-950°C  and  deposited  with  excess  rare  gas 
for  30-60  minutes.  The  Heliplex  was  used  to  cool  the  sample  window  for  all  Ag 
containing  matrices.  To  prepare  gold  containing  matrices  the  Knudsen  cell  technique  was 
again  used  but  Au  wire  (Spex)  was  used  instead  of  metal  powder.  This  cell  was  heated 
to  1200-1250°C  during  depositions  which  also  lasted  from  30-60  minutes.  A  Displex 
was  used  to  maintain  the  cold  window  during  all  Au  experiments.  These  matrices  were 
annealed  (to  -20K  for  Ar  and  -25K  for  Kr  and  Xe)  prior  to  recording  data  to  eliminate 
possible  interferences  from  multiple  sites. 

For  Ag  matrices,  moments  were  calculated  directly  from  the  digitized  spectra.  The 
Au  matrices  showed  a  large  background  absorption,  which  is  most  likely  attributable  to 
Rayleigh  scattering  by  the  matrix  itself.  Below  220  nm  the  absorption  signal  declined 
rapidly.  This  was  due  to  a  decrease  in  reference  signal  caused  by  reflectivity  losses  in 
the  chopper  wheel  mirrors.  The  background  of  each  Au  absorption  spectrum  was 
corrected  by  fitting  the  baseline  to  a  6th  order  polynomial  and  subtracting  the  calculated 
fit.  Normal  moments  calculations  were  then  performed  on  the  baseline  corrected  spectra. 
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Results 

Spin-Orbit  Interactions  in  the  "P  State 

The  absorption  and  MCD  spectra  of  Ag  atoms  in  Ar,  Kr  and  Xe  matrices  are 
shown  in  Figs.  4.1-4.3  with  the  observed  band  positions  listed  in  Table  4.1.  Temperature 
dependent  C  terms  dominate  the  spectra  as  shown  in  the  figures.    The  Ag/Xe  system 

shows  the  four  band  structure  that  has  been  observed  previously.4,1  '  The  origin  of  this 
structure  is  still  unsettled  although  several  proposals  have  been  introduced  that  will  be 
discussed  later. 

The  SO  coupling  constant,  £,  for  Ag  in  the  three  rare  gases  was  extracted  from 
the  slope  of  a  plot  of  the  reduced  first  MCD  moment  vs.  \IT.  Fig.  4.4  shows  the  plots 
which  are  straight  lines  as  expected  from  (28).  In  Ar  and  Kr  matrices,  the  extracted  £s 
were  found  to  be  greater  (Ar:  783  cm-1,  Kr:  638  cm-1)  than  the  gas  phase  value20  of  614 
cm-1,  while  in  Xe  a  reduction  to  583  cm-1  was  noted.  These  values  are  also  listed  in 
Table  4.1. 

In  deriving  the  equations  for  determining  the  SO  constants  second-order  effects 
were  neglected.  The  importance  of  these  effects  can  be  assessed  in  two  ways.  First, 
using  the  relation 

<AA>V<A>0  =  asfxBB0;')2/kT  (41) 

second-order  effects  can  be  shown  to  be  small  if  the  t,  obtained  is  close  to  that  obtained 
from  (28).  The  £'  represent  the  combined  contributions  from  first  and  second-order 
effects.  The  values  shown  in  Table  4.1  show  that  the  t,  obtained  from  this  expression  are 
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Figure  4.1        Experimental  MCD  (top)  and  absorption  (bottom)  spectra  of  Ag  atoms  in 
a  At  matrix  at  5.9  K. 
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Figure  4.2        MCD  (top)  spectra  of  Ag  atoms  in  a  Kr  matrix  showing  temperature 
dependence  of  the  signal.   Absorption  (bottom)  spectra  at  5.8  K. 
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Figure  4.3        Experimental  MCD  (top)  and  absorption  (bottom)  spectra  of  Ag  atoms  in 
a  Xe  matrix  at  6.5  K. 
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Table  4.1     Observed  Absorption  Band  Positions  (in  cm    )  and  Parameters  Derived  from 

the  Moments  Analysis  for  Ag  Atoms  in  Rare  Gas  Matrices. 


Gasa 

At 

Kr 

Xe 

Excited  State 

2p 

M/2 

29552 

31606 

30902 

29842 

2p 

r3/2 

30473 

32733 

31807 

30481 

33278 

32216 

30611 
30888 

^/cm"1 

614 

783 

638 

583 

^VcnT1 

614 

766 

616 

524 

<E2>N0c/cm"2 

- 

1.6xl05 

1.5xl05 

0.7xl05 

<E2><9  /cm"2 

- 

4.9xl04 

2.0xl04 

l.lxlO4 

v^/cm"1 

- 

32.1 

35.8 

24.2 

aValues  from  Ref.  20. 
b 


First  order,  obtained  from  (28). 
cSecond  order,  obtained  from  (41). 
d Values  from  fitting  (33). 
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Figure  4.4  Reduced  MCD  first  moment  plots  vs.  inverse  temperature  for  the  three 
Ag/RG  systems.  For  clarity,  the  Ag/Kr  and  Ag/Xe  plots  have  been  shifted 
by  -5  and  -10  y  axis  units,  respectively. 
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within  experimental  error  (-10%)  of  those  obtained  from  (28).  Secondly,  while  the  zeroth 

MCD  moment  is  expected  to  be  zero  to  first  order  it  should  become  nonnegligible  and 

vary  linearly  with  \IT  if  second-order  effects  are  present.   Plots  of  <AA>(>/<A>0  vs.  T"1 

with  'Bq/'Dq  (intercept)  values  of  the  order  10"9  -  10~10  and  Cq/2?0  (slope)  values  of  about 

10      indicate  that  second-order  SO  coupling  is  negligible  for  Ag  in  Ar,  Kr  and  Xe 
matrices. 

Baseline  corrected  absorption  and  MCD  spectra  for  Au  in  the  three  gases  are 
shown  in  Figs.  4.5-4.7  and  are  again  dominated  by  temperature  dependent  C  terms.  The 

band  positions,  listed  in  Table  4.2,  agree  well  with  those  reported  by  Gruen  et  al. 
Straight  lines  were  obtained  for  plots  of  <AA>1/<A>0  vs.  1AT,  shown  in  Fig.  4.8,  from 
which  the  SO  coupling  constants  were  extracted.  These  constants  are  listed  in  Table  4.2. 
In  this  case  the  extracted  values  (Ar:  3354  cm-1,  Kr:  3110  cm-1,  Xe:  2665  cm-1)  were 
all  greater  than  the  gas  phase20  value  of  2543  cm-1. 

Second-order  effects  were  again  assessed  using  (41)  which  is  shown  plotted  in 
Fig.  4.9.  The  t,  values  obtained  are  listed  in  Table  4.2  and  are  not  close,  for  Au/Ar  and 
Au/Kr,  to  the  first  order  values  obtained  from  (28).  For  Au/Xe  the  difference  in  the 
calculated  ts  is  within  the  experimental  error  (-4%)  but  the  difference  for  Au/Kr  is  -17% 
and  is  near  30%  for  Au/Ar.  This  is  also  reflected  in  the  zeroth  MCD  moment  plots 
shown  in  Fig.  4.10.  While  considerable  scatter  is  present  in  all  three  plots  the  Au/Xe  and 
Au/Kr  plots  appear  temperature  independent  but  the  Au/Ar  plot  shows  a  distinct 
temperature  dependence  (i.e.,  linear  with  1/T).  The  importance  of  second-order  spin  orbit 
effects  appears  to  decrease  in  going  from  Ar  to  Xe  matrices. 
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Figure  4.5        Experimental  MCD  (top)  and  absorption  (bottom)  spectra  of  Au  atoms  in 
a  At  matrix  at  10.5  K. 
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Figure  4.6        MCD  (top)  spectra  of  Au  atoms  in  a  Kr  matrix  showing  temperature 
dependence  of  the  signal.   Absorption  (bottom)  spectra  at  10.3  K. 
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Figure  4.7        Experimental  MCD  (top)  and  absorption  (bottom)  spectra  of  Au  atoms  in 
a  Xe  matrix  at  12.0  K. 
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Table  4.2    Observed  Absorption  Band  Positions  (in  cm  *)  and  Parameters  Derived  from 
the  Moments  Analysis  for  Au  Atoms  in  Rare  Gas  Matrices. 


Gasa 

Ar 

Kr 

Xe 

Excited  State 

2p 

37359 

38895 

38124 

36644 

2p 
r3/2 

41174 

43365 
43937 

42301 
42790 

40258 
40734 

^/cra"1 

2543 

3354 

3110 

2655 

^/cm"1 

2543 

2357 

2572 

2555 

a:bd 

- 

0.44:0.90 

0.85:0.52 

0.98:0.20 

tf/caT1 

- 

2120 

2577 

2534 

<E2>N°c/cm-2 

- 

1.2xl06(g) 

1.3xl06(g) 

0.3xl06 

<E2><5  /cm"2 

- 

-1.7xl06 

-l.lxlO6 

-O.lxlO6 

v^/cm-1 

- 

- 

- 

23.2 

aValues  from  Ref.  20. 
bFirst  order,  obtained  from  (28). 
cSecond  order,  obtained  from  (41). 
Coefficients  of  mixing:   ax|>'(2P3/2)  +  bvJ>"(4P3/2)- 
elst  SO  parameter,  calculated  from  (44). 
Values  from  fitting  (33). 
8Average  values  -  no  fit  performed. 
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Figure  4.8  Reduced  MCD  first  moment  plots  vs.  inverse  temperature  for  the  three 
Au/Rg  systems.  For  clarity,  the  Au/Kr  and  Au/Xe  plots  have  been  shifted 
-5  and  -10  y  axis  units,  respectively. 
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Figure  4.9        Reduced  MCD  second  moment  plots  vs.  inverse  temperature  for  the  three 
Au/Rg  systems. 
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Figure  4.10      MCD  zeroth  moment  plots  vs.  inverse  temperature  for  the  three  Au/Rg 
systems. 
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The  fact  that  second-order  SO  interaction  in  Au/RG  decreases  in  the  order  RG  = 

Ar  >  Kr  >  Xe  shows  the  dramatic  effect  the  matrix  material  can  have  on  this  quantity. 
Second-order  SO  coupling  may  arise  from  either  out-of-state  mixing  of  the  Au  6p 
orbital  with  the  appropriate  orbitals  on  nearby  rare  gas  atoms  and/or  with  other  Au  atom 
excited  states.  Mixing  with  higher-lying  p  orbitals  on  the  rare  gas  neighbors  can  be 
immediately  discounted  since  the  trend  in  the  energies  of  the  appropriate  mixing  states" 
is  opposite  to  the  observed  trend.  Table  4.3  lists  the  energies  of  the  states  for  Au.  Note 
that  the  energy  of  the  2P3/2  multiplet  (the  average  of  the  two  bands)  increases  for  Au/Ar 
and  Au/Kr  compared  to  the  gas  phase,  whereas  for  Au/Xe  it  decreases.  This  multiplet 
may  second-order  SO  couple  to  the  4P3/2  multiplet.  Coupling  of  the  P1/2  state  with  the 
?l/2  state  was  neglected  here  due  to  the  significantly  larger  energy  separation.  The  4P3/2 
state  was  assumed  to  be  unshifted  in  the  various  matrices  because  no  information  on  its 
location  was  available  from  this  study  (it  lies  at  shorter  wavelengths  than  the  present 
apparatus  can  reach).  It  can  then  be  seen  that  the  energy  spacing  between  the  P3/2  ar>d 
4P3/2  multiplets  decreases  in  the  order  Xe  (-6500  cm"1)  >  Kr  (-4500  cm-1)  >  Ar  (-3300 
cm"  ).  On  this  basis,  a  larger  second-order  SO  coupling  would  be  predicted  in  Au/Ar 
than  in  Au/Xe.  By  way  of  contrast,  the  spacing  of  the  Ag  2P3/2  and  4P3/2  multiplets  is 
in  the  28,000  cm-1  range,  effectively  precluding  any  appreciable  second-order  mixing. 
The  second-order  SO  mixing  between  excited  states  of  a  matrix-isolated  atom  or 
ion  may  be  explicitly  determined  using  a  procedure  that  has  been  previously  reported. 
Rivoal  and  Vala  and  coworkers39,47  showed  that  for  the  mixing  of  an  allowed  state  u/  and 
a  forbidden  state  \\>"  with  mixing  coefficients,  a  and  b,  respectively,  the  zeroth  MCD 
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Table  4.3         Selected  Excited  States,  Configurations  and  Gas  Phase  Energies  of  Au 

Atomsa 


Au 

State 

Configuration 

v/cm 

Gas 

At 

Kr 

Xe 

JSl/2 
^Pl/2 

"P3/2 

4P 

/3/2 

?W 
Avc 

5d106s 

5d106p 

5d106p 

5d96s6p 

5d96s6p 

0 
37358.9 
41174.3 
47004.4 
53197.5 

0 
43651b 

3353 

0 
42546b 

4458 

0 
40496 

6508 

aFrom  Ref.  20 

Taken  as  average  of  two  split  bands,  cf.  Table  4.2. 
cAv  ^  v  GAS  (4P3/2)  -  vRG  (2P3/2),  in  cm"1. 
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moment,  <AA>q  for  the  allowed  band,  is  given  by 

<AA>a  __  (2J+3)(a2-l)m2  (42) 

0  (J+1)(2J+1)         ^ 

where  J  is  the  total  angular  momentum  of  the  ground  state,  m  is  the  transition  dipole 
moment  matrix  element  and  <SZ>  is  the  thermal  average  over  the  Zeeman-split  ground 
state  levels.  The  value  of  m  was  obtained  from  <A>0  =  m  ,  the  zeroth  absorption 
moment  over  the  entire  band.  <AA>q  and  <S^>  were  calculated  at  each  temperature 
and  were  combined  with  m"  in  (42)  to  give  a  value  for  a^.  The  mixing  coefficient  a  was 
taken  from  the  square  root  of  the  average  of  all  a2  values.  From  a  +b  =1,  the  value  of 
b  was  calculated.   Finally,  using  the  perturbed  first  moment  MCD  expression    • 

<AA>*E  =  [2gtu.BB  +  2^<Sz>]m2  <43) 

where 


C  -t 


j  +  b2,     r2J,3)     /AW  _  j) 
2   (J+1X2J+1)(   ^         J 


(44) 


the  first  order  SO  constant  t,  for  the  allowed  2P  state  was  extracted  from  the  usual  MCD 
first  moment  plot.  The  experimental  SO  constants  (t,',  which  include  first  and  second 
order  contributions)  for  the  Au/Ar,  Au/Kr  and  Au/Xe  cases  have  been  given  previously 
(cf.  Table  4.2).  Using  b  and  AW  (the  energy  separation  between  allowed  and  forbidden 
partners)  along  with  (44)  the  first  order  SO  constant  t,  was  calculated.  The  values  found 
for  the  a  and  b  coefficients  and  the  £'  and  t,  constants  are  given  in  Table  4.2.  If  the 
mixing  partners  proposed  here  are  valid,  the  £  values  obtained  should  be  close  to  the  ones 
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determined  from  (41).    These  values  were  found  to  lie  within  3-10%  of  each  other: 

Au/Ar  (2120  cm"1  vs.  2357  cm-1),  Au/Kr  (2577  cm"1  vs.  2572  cm"1),  and  Au/Xe  (2534 

cm"1  vs.  2555  cm"1).   The  a  and  b  coefficients  for  Au/Ar  (a  =  0.44,  b  =  0.90)  indicate 

that  the  "P3/2  and  4P3/2  are  thoroughly  mixed  by  SO  interaction.   The  slight  discrepancy 

between  the  two  t,  values  may  arise  from  other  states  such  as  2P1/2>  4Pi/2  not  included  in 

the  calculation  or  from  the  inaccuracy  of  the  assumed  energy  separations  (AW).   In  any 

event,  the  close  correspondence  found  strongly  indicates  that  the  second-order  SO  effect 

arises  primarily  from  the  mixing  between  the  T^  and  4P3/2  states.  This  also  establishes 

that  the  intra-atomic  SO  coupling  in  gold  is  strongly  influenced  by  the  matrix  material 

with  the  largest  effect  occurring  in  Ar  and  the  smallest  in  Xe. 

Origin  of  the  Observed  Band  Structure 

Two  primary  explanations  have  been  offered  for  the  appearance  of  three  bands  in 
the  spectrum  of  the  P  *-  S  transition  of  matrix  isolated  noble  and  alkali  metals:  a  static 
distorted  site  model  or  a  dynamic  Jahn-Teller  (JT)  effect.  Lund  et  al.24  have  shown  that 
absorption  spectra  and  MCD  spectra  may  be  combined  to  conclusively  establish  the 
dynamic  JT  effect  as  the  source  of  the  three  bands  for  the  matrix  isolated  alkali  metals. 
The  present  spectra  for  matrix  isolated  Ag  and  Au  atoms  differ  from  the  alkali  metals  in 
several  ways.  The  measured  SO  constants  for  Ag  and  Au  are  much  larger  than  the  alkali 
metal  SO  constants  and  do  not  change  sign  from  their  gas  phase  values,  indicating  that 
the  2P3/2  multiplet  remains  at  higher  energy  than  the  2P1/2  multiplet.  Also,  the  2P1/2  band 
is  well -resolved  and  separated  from  the    P3/2  doublet. 
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Previous  authors  have  proposed  explanations  for  the  observed  band  structure  of 

matrix-isolated  Au  atoms.   Gruen  and  coworkers     proposed  that  an  axial  site  distortion 

was  responsible  for  the  splitting.    Forstmann,  Kolb,  Leuthoff  and  Schulze,13  following 

Gruen,  also  proposed  such  an  origin  for  the  appearance  of  the  three  peaks.    Weinert  et 

al.,     used  a  moments  analysis  on  the  absorption  bands  of  the  individual  multiplets  within 

the  extended  configuration  coordinate  model,  but  noted  a  number  of  problems  with  this 

approach.    They  suggested  in  a  subsequent  paper4    that  a  moments  analysis  utilizing 

MCD  and  absorption  would  yield  valuable  information.    In  the  following,  a  distorted  site 

model  will  be  tested  and  shown  to  be  incapable  of  reproducing  the  observed  spectra,  and 

thereafter,  a  moments  analysis  of  the  MCD  and  absorption  data  will  be  presented  to 

confirm  the  presence  of  JT  activity. 

Following  the  procedure  used  by  Lund  et  al.,24  the  metal  atom  was  assumed  to 

reside  in  a  site  of  D4  (or  D3)  symmetry.   A  crystal  field  calculation  on  the  2P  state  then 

yielded  the  energies  of  the  three  split  levels  as  a  function  of  £,  the  one  electron  SO 

coupling  constant,  and  A,  the  crystal  field  (or  tetragonal10)  parameter  as 


E^bE')  =  -l(-£+A+2S) 


E^E")  =    I(-£+A)  (45) 

E3(aE')  =  -I(-£+A-2S) 


where 
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S  =  (^A2+^A£+:^2)2  =E,-Et>0 

4  2  4  3     l 


(46) 


The  state  descriptions  in  parentheses  refer  to  the  D4  double  group  (D4).  For  the  correct 
ordering  of  energy  levels  for  the  Ag/RG  and  Au/RG  cases,  where  the  P3/2  (U')  SO 
multiplet  lies  above  the  2P1/2  (E')  level,  the  expressions  of  Lund  et  al.  have  been 
modified  by  substituting  -£,  for  t  and  -A  for  A.  The  energy  levels  are  shown 
schematically  in  Fig.  4.11.    Solving  the  above  equations  for  t,  and  A  yields 


< 

=  a-(a2 

l 

A  = 

=  a+(a2- 

1 

-q)1 

(47) 


(48) 


where 


a  =[(E2-E3)HZ2-Ei)V3 


(49) 


q  =2[(E2-E1)(E2-E3)]/3.  (50) 

The  relevant  MCD  parameters,  (A  v  "20,  C0,  ©0),  defined  in  Lund  et  al.,     were  calculated 
and  are  given  for  Ag/RG  and  Au/RG  in  Table  4.4. 

Using  these  parameters  the  MCD  and  absorption  spectra  were  simulated  from 


AA'/e  =  152.5  £ 

i=l 


-*i(i) 


3fi(e) 


de 


«oO) 


(kT)j 


m 


c£ 


(51) 


and 
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3 


A/e  =  326.6  £   V0(\)ix(t)c<l  (52) 

i=i 

where  AA'  is  the  MCD  per  tesla,  e  is  the  photon  energy  (in  cm  ),  c  the  sample 
concentration  and  2  its  pathlength.  The  band  shape,  f,(e),  was  taken  as  a  normalized 
Gaussian 

fi(e)  =(l/&lV^)exp(-(e-eI)2/&f).  (53) 

As  seen  in  Table  4.4,  the  MCD  parameters  Av  22  0,  and  CQ  are  completely  determined  by 
the  two  crystal  field  parameters,  t,  and  A,  which  in  turn  are  determined  by  the  absorption 
band  peak  positions.  Most  importantly,  the  scaling  of  the  simulated  MCD  spectrum  is 
completely  fixed  by  the  experimental  absorption  band  intensity.  For  the  present 
simulations,  both  simulated  and  experimental  absorption  spectra  were  normalized  on  the 
lowest  energy  (bE')  band. 

Using  this  procedure  absorption  and  MCD  spectra  were  simulated  for  both  Ag  and 
Au  atoms  in  the  three  rare  gases.  The  best  fit  simulations  are  compared  with 
experimental  spectra  in  Fig.  4.12  (Ag/Ar)  and  Fig.  4.13  (Au/Xe).  The  bandwidths  for 
each  band  (i=l-3)  were  adjusted  to  best  match  the  absorption  spectra  contours  and  then 
retained  when  simulating  the  MCD  spectra.  For  Ag/Ar  the  absorption  spectrum  is 
matched  reasonably  well  except  for  the  highest  energy  band  whose  intensity  is  not  well- 
reproduced.  The  Ag/Ar  MCD  simulation  shows  three  distinct  bands,  as  observed,  all  of 
which  are  predicted  to  be  too  intense.  Similar  observations  hold  for  the  Ag/Kr  and  Ag/Xe 
cases,  with  absorption  being  simulated  well  but  with  even  larger  discrepancies  between 
predicted  and  observed  MCD  intensities. 
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Figure  4.11      Ordering  of  the  spin-orbit  split  levels  of  the  excited  2P  state  in  octahedral 
(Oh)  and  distorted  tetragonal  (D4)  fields. 


Table  4.4         Crystal  Field  Model:   Absorption  and  MCD  Parameters3 
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E'(2S)  — 

bE' 

E" 

aE' 

*,» 

x+y 

1/3 

x-y 

7d  b 

®21  +  «31 

-%2l  +  232 

-*31  "  «32 

r  b 

1/6  +  u 

-1/3 

1/6 -u 

©0b 

1/3 

1/3 

1/3 

^is  table  is  taken  from  Ref.  24. 

bAll  MCD  and  absorption  parameters  are  in  units  of  the  reduced  transition  moment  matrix 
element  squared.   The  symbols  in  the  Table  are  defined  as: 


x  ■  5g  +  5?g  +  1-  y 

•  2  „02  6' 


5Ccf 


'CF 


2S 


12S 


CF 


CF 


4S 


^CF  ....    3  ^CF 

-,    U  —  . 


iCF 


CF 


4S 


CF 


4S 


CF 


9ACF       SAqjCcf       9^ 
S   =  =  +  =  E3-EJ 


2j:'s  given  in  Ref.  24. 
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Figure  4.12  MCD  (top)  and  absorption  (bottom)  spectra  of  Ag  atoms  in  a  Ar  matrix  at 
6K  (solid  curves).  Simulation  using  the  CF  model  with  variable 
bandwidths  (dotted  curves).  Simulation  using  the  CF  model  with 
bandwidths  determined  from  eq.  54  (dashed  curves).  Simulations  are 
normalized  to  lowest  energy  experimental  absorption  band. 
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Figure  4.13  MCD  (top)  and  absorption  (bottom)  spectra  of  Au  atoms  in  a  Xe  matrix 
at  10.5K  (solid  curves).  Simulation  using  the  CF  model  with  variable 
bandwidths  (dotted  curves). 
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In  the  case  of  Au/Xe,  similar  results  were  obtained.  The  intensities  of  both  higher 

energy  bands  (E"  and  aE')  were  not  reproduced  well  by  the  crystal  field  model.    The 

MCD  simulation  again  shows  intensities  which  are  too  strong  for  all  bands.  Furthermore, 

the   bandwidths  necessary  to  produce  the  absorption  simulation  result  in  a  MCD 

simulation  with  only  two  distinct  bands.  The  two  higher  energy  bands  show  no  evidence 

of  separation  in  the  simulated  MCD.    This  model  produced  considerably  less  accurate 

simulations  for  the  other  Au/RG  cases.  Overall,  the  crystal  field  model  does  not  correctly 

simulate  the  amplitudes  of  the  absorption  and  MCD  spectra  of  any  of  the  Ag/RG  or 

Au/RG  cases. 

Another  approach  to  the  crystal  field  model,  introduced  by  Lund  et  al.,     starts 

with  a  specific  D4  (or  D3)  model  and  develops  moment  equations  for   MCD  and 

absorption,  which  are  then  equated  to  the  experimental  moments.    This  approach  was 

applied  to  Ag  in  the  rare  gases.   The  parameters,  £  and  A,  are  then  deduced  from  these 

moments,    the    MCD    and   absorption   are   simulated   and   finally   compared    to    the 

experimental  spectra.  Lund  et  al.24  found  that  (28)  is  still  valid  if  gorb  is  set  equal  to  one 

and  2<E2>NC  is  replaced  by  A2,  where  A  is  the  tetragonal  crystal  field  parameter  defined 

in  (48).  The  standard  first  MCD  moment  plot  was  used  to  deduce  t-  The  value  of  A  was 

taken  as  the  square  root  of  2<E2>°NC  obtained  from  (29)  and  (30),  vide  infra.    These 

values  have  been  listed  in  Table  4.2.     With  £  and  A,  the  band  energies  (Ej)  were 

calculated  using  (45)  and  thence  the  A  j  (i),  #  0  (i),  C  0  (i)  and  i>  0  (i)  obtained  from  the 

definitions  in  Table  4.4.   The  bandwidths  were  not  varied  arbitrarily  to  fit  the  absorption 

widths  (as  in  the  previous  model),  but  instead  the  relation  used  was 
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£6f=6<EV  (54) 

i 

This  form  is  dictated  since  the  only  source  of  line  broadening  in  the  crystal  field  model 
comes  from  cubic  lattice  modes  (i.e.,  <E2>C)  and  because  the  dipole  strengths  T>  0  (i)  for 
transitions  to  each  of  the  levels  are  identical.  The  values  of  <E  >c  deduced  from  (29) 
and  (30)  are  given  in  Table  4.2.  The  simulations  using  bandwidths  determined  from  (54), 
shown  as  a  dashed  curve  in  Fig.  4.12  for  Ag/Ar,  generally  fit  the  spectra  very  poorly. 
(This  procedure  was  not  performed  for  Au/RG  due  to  the  negative  values  of  <E  >c 
obtained  from  the  experimental  data,  vide  infra.)  Indeed,  neither  approach  using  the 
crystal  field  model  is  able  to  fit  both  MCD  and  absorption  well.  The  inability  of  the 
crystal  field  model  to  reproduce  the  observed  spectra  agrees  with  the  conclusions  of  Kasai 
and  McLeod8  that  the  atoms  sit  in  a  substitutional  site  and  requires  that  a  different 
explanation  for  the  tri-peaked  appearance  be  found. 

In  earlier  work  on  Cu  in  the  rare  gases    •     and  the  alkali  metals  in  the  rare 
gases,  which  also  feature  a   P  *-   S  transition,  it  was  conclusively  established  that 

the  Jahn-Teller  effect  is  operative  in  the  excited  2P  state  and  that,  together  with  the  SO 
splitting  of  the  1/2  and  3/2  J  components,  it  is  responsible  for  the  tri-peaked  appearance 
of  the  absorption  band.  The  assumption  that  a  similar  mechanism  holds  true  for  the  Ag 
and  Au  cases  is  warranted  since  a  matrix-dependent  SO  splitting  has  been  shown  and  the 
transition  involved,  "P  *-  "S,  is  identical.  The  proof  of  the  existence  of  JT  activity  in 
matrix-isolated  atoms  rests  on  two  major  experimental  observations,  viz.,  that  the 
contributions  of  noncubic  lattice  vibrations  to  the  bandwidths  (i.e.  <E2>NC)  are 
nonnegligible  and  that  their  values  vary  with  temperature  as  the  coth(Vj/2kT)  function. 
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To  test  this  assumption  (29)  and  (30)  were  used  to  determine  the  cubic,  <E2>C, 

9  9  9 

and  non-cubic,  <E  >NC,  contributions  to  the  bandwidths.  Plots  of  <E~>C  and  <E  >NC 
vs.  T  for  the  Ag/RG  cases  are  shown  in  Figs.  4.14-4.16.  The  <E2>NC  and  <E2>C  values 
found,  together  with  the  effective  lattice  frequencies  (vs),  obtained  by  fitting  (53),  have 
been  listed  in  Table  4.1.  The  <E2>C  terms  are  all  small  compared  to  the  <E2>NC. 
Furthermore,  all  the  <E  >NC  data  fit  a  coth  function  reasonably  well,  although  with 
considerable  scatter.  These  observations  fit  the  criteria  for  JT  activity  in  the  2P  state  and 

9  9 

support  the  assumption  that  the  tri-peaked  appearance  of  the  ""P  *—  S  transition  for 
matrix  isolated  silver  atoms  is  a  result  of  simultaneous  SO  splitting  and  a  dynamic  JT 
effect. 

Confirmation  of  the  JT  effect  in  the  Au/RG  cases  was  similar  but  not  without 

some  complications.    Figure  4.17  shows  plots  of  <E2>NC  and  <E2>C  as  a  function  of 

9     n 
temperature  for  Au  atoms.    For  Au/Xe,  expression  (33)  was  fit  well  with  <E  >NC=  3.0 

x  10    cm      and  Vj  =  23  cm    .    However,  the  <E  >c  values  for  Au/Xe  are  negative  (a 

physically  unrealistic  result)  and,  furthermore,  became  increasingly  negative  at  higher 

9  9 

temperatures.  For  Au/Ar  and  Au/Kr,  the  <E  >NC  and  <E  >c  plots  are  problematic.  They 
show  practically  no  T-dependence  although  their  values  for  <E2>NC  are  greater  than  for 
<E  >q.  These  anomalous  results  may  arise  from  anharmonicity  and/or  higher  order 
electron  phonon  coupling  terms  not  specified  in  the  moments  derivation  and  are  discussed 
further  below.  However,  since  the  crystal  field  model  could  not  correctly  simulate  the 
observed  absorption/MCD  (recall  that  the  Au/Xe  case  was  better  simulated  than  either  the 
Au/Ar  or  Au/Kr  case)  and  the  calculated  <E  >NC  values  are  substantial,  it  appears  that 
the  JT  effect  is  active  in  the  three  Au/RG  cases. 


77 


35 

30 

1 

25 

E 

20 

o 

15 

Ld 

10 

5 

0 

-5 


0 


Ag/Ar 


o    <EZ> 


E™WDQpDDDnD    O 


_i L 


J i I i I 


10 


20 


T/K 


30 


40 


Figure  4.14  Temperature  dependence  of  the  cubic  <E2>C  and  noncubic  <E2>NC 
bandwidth  contributions  of  Ag/Ar.  Solid  line  fit  follows  the  relation: 
<E2>NC  =  1.6X105  coth(32.1/2kT),  in  cm-1. 
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Figure  4.15      Temperature  dependence  of  the  cubic  <E  >c  and  noncubic  <E  >NC 

bandwidth  contributions  of  Ag/Kr.    Solid  line  fit  follows  the  relation: 
<E2>NC  =  UxlO5  coth(35.8/2kT),  in  cm"1. 
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Figure  4.16      Temperature   dependence  of  the  cubic  <E  >c  and  noncubic  <E  >NC 

bandwidth  contributions  of  Ag/Xe.    Solid  line  fit  follows  the  relation: 
<E2>NC  =  6.9xl04  coth(24.2/2kT),  in  cm"1. 
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Figure  4.17      Temperature     dependence    of    the    cubic    and    noncubic    bandwidth 
contributions  for  the  three  Au/RG  systems. 

Legend:   (<E2>NC  values:   O  Au/Ar;  ■  Au/Kr;  a  Au/Xe;  <E2>C  values: 
+  Au/Xe;  T  Au/Kr;  O  Au/Ar). 
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The  expressions  given  in  (28)  -  (34)  were  derived  under  a  number  of  limiting 

assumptions:   harmonic  ground  state  vibrational  wavefunctions,  first  order  SO  coupling, 

and  linear  electron-phonon  coupling.37    In  relaxing  the  latter  condition  and  allowing 

"2*7 

quadratic  electron-phonon  coupling,  expression  (29)  becomes 


<A>2/<A>0  =  £2/2  +  Iw2  -  lo? .  (55) 

2  4 


Here  o)2  involves  only  linear  electron-phonon  coupling  and  is  equivalent  to  (34).  The 
additional  term,  -l/4o)j  ,  incorporates  the  influence  of  quadratic  electron-phonon  mixing, 
which  is  given  by 

(!>!  =J2(Ki/2Ai)v1coth(v1/2kT)  (56) 

i 

where  K4  is  the  quadratic  electron-photon  coupling  parameter  and  other  terms  are  as 
defined  previously.  The  magnitude  of  (Oj  may  be  determined  from  the  variation  with  T 
of  the  barycenter  energy,  v0  (as  defined  in  (26)): 


v0  ■  vp  +  ^1  s  vp  +  Npcoth(v,/2kT)  (57) 


A  plot  of  v0  vs.  T  is  shown  in  Fig.  4.18  for  Au/Xe.  Knowing  the  magnitude  of  0)j  at 
various  temperatures  now  permits  a  qualitative  understanding  of  the  form  of  the  Au/RG 
<A>2/<A>0  vs.  T  curves.  Figure  4.19  shows  that  Au/Ar  decreases  at  higher  T,  while  both 
Au/Kr  and  Au/Xe  increase  at  higher  T.  This  behavior  can  be  interpreted  as  arising  from 
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Figure  4.18      Band  peak  position  vs.  temperature  for  Au/Xe.   The  solid  line  fit  follows 
the  relation:  v0  =  39283  -  148coth(29.5/2kT),  in  cm"1. 
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Figure  4.19      Reduced  absorption  second  moment  plots  vs.  temperature  for  the  three 
Au/RG  systems.    Lines  are  least  squares  fits. 
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a  partial  cancellation  of  the  a>2  and  (o1  terms,  with  the  (Oj  term  larger  in  Au/Ar  and  the 

co2  term  larger  in  the  Au/Kr  and  Au/Xe  samples.    However,  as  expected,  the  SO  term 

dominates  in  all  three  cases.    The  intercept  (T=0)  in  these  plots  should  equal  £2 72,  and 

does.    The  values  observed  are  (in  parenthesis  the  t,  values  from  first  moment  MCD 

slopes):  Au/Ar:  3350  cm-1  (3354  cm"1);  Au/Kr:  3200  cm"1  (3110  cm"1);  Au/Xe:  2650 

cm"1  (2655  cm"1). 

The  question  of  why  the  <E2>C  curves  are  negative,  and  increasingly  so  at  higher 

temperatures,  can  be  assessed  qualitatively  in  the  following  manner.   Labelling  the  LHS 

of  (29)  as  R  and  the  ratio  of  the  temperature  dependent  parts  of  (30)  and  (28)  as  P,  i.e., 

R  s  <A>2/<A>0  =  <E2>C  +  <E2>NC  +  £2/2  (58) 

P  =  <AA>3/<AA>!  =  3<E2>C  +  3<E2>NC/2  +  3£2/4  (59) 

and  rearranging,  gives  expressions  for  <E2>C  and  <E2>NC: 

<E2>NC  =  2R-2P/3-C2/2  (60) 

<E2>C  =2P/3-R  (61) 

Thus  if  R  >  2P/3,  <E  >c  can  become  negative.  The  question  is  now,  under  what 
conditions  can  <E2>C  be  expected  to  be  negative.  Recall  that  R  =  £2/2  +  l/2co2  -  1/4(0^ 
(cf.  (55)).  Knowing  that  (ol  »  -304.6  cm"1  at  10K  (from  the  v0  vs.  T  plot)  and  that  its 
(negative)  contribution  to  R  is  only  -3%  of  the  total  value  (at  50K  it  climbs  to  -16%), 
the  value  of  R  will  increase  with  T,  but  not  sufficiently  if  only  linear  electron-coupling 
terms  are  included  in  the  co2  term.  The  fitting  of  the  <A>2/<A>0  curve  shows  this  is  to 
be  the  case.  Thus,  it  appears  clear  that  higher  order  coupling  effects  and/or  anharmonicity 
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terms  should  be  included  in  w2  to  account  for  the  further  decrease  of  <E^>C  below  zero 

with  increasing  temperature.  This,  of  course,  demonstrates  that  the  values  of  <E">NC  and 

<E2>C  obtained  above  are  only  approximate.   Still,  the  determination  of  non-negligible 

<E  >NC  values  provides  strong  evidence  that  the  Jahn-Teller  effect  is  operative  in  these 

systems. 

The  appearance  of  four  bands  in  the  Ag/Xe  absorption  and  MCD  spectra  deserves 

further  comment.    The  crystal  field  (CF)  model  by  itself  is  incapable  of  explaining  this 

observation.   Previous  authors     have  proposed  that  the  fourth  band  arises  from  the  CF- 

split    P3/2  ±  3/2,  level  which  interacts  with  the  Xe  6s  level  in  a  Ag-Xe  molecule.   Others 

ID  'J 

have  suggested  that  it  is  the  result  of  a  resonant  interaction  of  the  P3/2  multiplet  with 
a  near-lying  ^D5/2  state,  coupled  by  nontotally  symmetric  vibrations  of  the  lattice. 
Having  confirmed  the  importance  of  the  JT  effect  in  this  sample,  it  is  possible  to 
speculate  that  the  additional  structure  arises  directly  from  transitions  to  the  three  JT-SO 
excited  state  potential  manifolds  whose  inherent  sharp  band  structure  is  not  completely 
smoothed  out  by  coupling  to  the  symmetric  modes.  A  computation  of  the  band  contours 
which  may  lead  to  this  situation  would  be  very  interesting  and  would  possibly  lead  to  a 
conclusive  argument  for  the  origin  of  the  fourth  band. 

Discussion 

Spin-Orbit  Coupling 

Recently,  a  theory  of  spin-orbit  coupling  in  the  T  excited  state  of  alkali  metal 
and  noble  metal  atoms  in  rare  gas  matrices  was  proposed.33  The  theory  assumed  a  model 
in  which  the  metal  atom  resides  at  a  tetradecahedral  substitution  site  surrounded  by  a  first 
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shell  cage  of  twelve  rare  gas  atoms.     It  was  found  that  with  this  "super-molecule" 

approach  the  magnitude  and  sign  of  the  SO  splitting  in  the  matrix  was  controlled  by  three 

factors:  the  metal-rare  gas  overlap  integral  (metal  p  orbital  with  rare  gas  p  orbitals),  the 

molecular  orbital  mixing  coefficients,  and  the  ratio  of  the  rare  gas-to-metal  SO  coupling 

constants.  Three  categories  of  experimental  SO  coupling  constants  were  recognized:   (I) 

positive  and  greater  than  the  gas  phase,  (II)  positive  and  less  than  the  gas  phase,  and  (III) 

negative  and  less  than  the  gas  phase  value.     Using  computed  overlap  integrals  and 

empirical  SO  coupling  constants,  the  theory  was  shown  to  predict  these  three  categories. 

Quantitative  comparisons  were  attempted  using  a  semi-empirical  M.O.  approach  to  find 

the  M.O.  mixing  coefficients.    At  the  time,  comparison  for  the  noble  metals  was  only 

possible  for  Cu/RG.  A  match  was  attempted  for  Au/RG  for  which  t,  was  extracted  from 

the  absorption  band  positions,     a  method  approximately  correct  because  of  the  large  SO 

coupling  constant  of  Au.     Ag/RG  values  were  not  available  when  the  theory  was 

proposed.    Table  4.5  shows  the  experimental  SO  "reduction"  factors  (i.e.,  Crg^gas)  anc* 

compares  them  to  two  versions  of  the  theory.    The  subscripts  "np"  and  "ave"  on  the 

reduction  factors  in  Table  4.5  refer  to  the  choice  of  RG  orbital  used  for  the  calculation 

of  the  M.O.  mixing  coefficient:   "np"  denotes  the  (n+l)p  orbital  (unoccupied)  and  "ave" 

means  the  energy-weighted  average  of  the  np  (occupied)  and  (n+l)p  orbitals  on  the  rare 

gas. 

Predictions  for  all  the  Au/RG  cases  for  both  the  np  and  ave  approaches  fall  into 

Class  I  as  do  all  experimental  values.    For  Au/Ar  and  Au/Kr  the  predicted  values  are 

somewhat  low.     For  both  k      and  kave  Ag/Ar  is  predicted  to  be  in  Class  I  and  is 

experimentally  shown  to  be  so.  Ag/Kr  is  predicted  to  be  either  in  Class  I  (kave)  or  Class 
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II  (knp)  and  is  experimentally  in  Class  I.   Ag/Xe  is  similar,  but  now  experiment  shows 

it  to  be  in  Class  II.    Presumably,  a  modification  of  the  weighting  factors  of  the  average 

orbital  approach  would  reproduce  the  observed  sequence.     The  Cu/RG  systems  are 

included  here  for  completeness  sake  although  they  have  been  discussed  previously.33  The 

k      approach  provides  the  best  results  but  fails  to  predict  the  negative  k  observed 

experimentally.  The  supplemental  column  labelled  k^  denotes  a  calculation  for  a  matrix 

cage  arbitrarily  collapsed  25%.   With  shorter  metal-rare  gas  distances  it  can  be  seen  that 

the  negative  sign  is  reproduced.   The  SO  constants  derived  here  provide  a  rigorous  test 

of  the  proposed  theory,  which  accounts  for  many,  but  not  all,  of  the  observed  features  of 

the  noble  metal  SO  splittings  in  the  rare  gases. 

By  the  use  of  various  moments  of  the  MCD  and  absorption  spectra  the  presence 

of  second-order  SO  coupling  between  the    P3/2  and  4P3/2  states  of  Au/RG  has  been 

detected  and  quantified.    The  extracted  intra-atomic  mixing  coefficients  show  significant 

mixing  in  both  Au/Ar  and  Au/Kr  cases.    This  mixing  was  suggested  by  Gruen  and 

coworkers  in  1974  as  an  explanation  for  a  strong  "two-electron"  emission  transition  of 

gaseous  Au  observed  between  the  2P  (5d106p)  and  2D(5d96s2)  states.10  Mixing  of  the  2P 

state  with  others  of  5d96s2  configuration  was  put  forth.   The  former  coupling  is  exactly 

what  has  been  described  here,  except,  of  course,  that  here  the  coupling  is  enhanced  by  the 

matrix-induced  shift  of  the  ^P  level. 

Jahn-Teller  Activity 

An  investigation  of  the  distorted  site  model  as  the  explanation  of  the  tri-peaked 
band  structure  showed  that  this  model  for  the  matrix  isolated  species  could  not  simulate 
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the  observed  MCD  and  absorption  spectra  well.    It  is  also  worth  noting,  as  Lund  et  al. 

have  done  for  the  Li/RG  cases,"4  that  the  £  parameters  extracted  from  the  distorted  site 
model  are  very  different  from  those  obtained  from  the  MCD  first  moment  plots  (e.g.,  for 
Au/Ar  the  distorted  site  model  gives  t,  =  542  cm-1  whereas  the  first  moment  plot  yields 
3354  cm").  The  values  obtained  from  the  moment  plots  cannot  be  seriously  in  error 
since    they    are    practically    model-independent.  Nevertheless,    even    when    the 

experimental  t,  values  are  used  in  the  distorted  site  model,  good  simulations  of  MCD  and 
absorption  band  contours  were  impossible.  This  evidence  confirms  that  the  tri-peaked 
appearance  of  the  observed  spectra  does  not  result  from  a  distorted  matrix  site. 

On  the  other  hand,  it  has  been  demonstrated  that  the  JT  effect  is  active  in  both 
Ag/RG  and  Au/RG  systems.  The  contributions  of  totally  symmetric  modes  in  the  Ag/RG 
systems  were  shown  to  be  three  to  six  times  smaller  than  the  nontotally  symmetric  ones. 
Noncubic  contributions  to  the  mean  square  band  widths  were  also  shown  to  be  substantial 
for  the  Au/RG  systems.  The  observation  of  negative  cubic  contributions  was  ascribed  to 
the  presence  of  higher  order  electron-phonon  coupling  or  anharmonicity  terms  not 
included  in  the  derivation  of  the  moment  equations  used.  Although  this  casts  doubts  on 
the  values  extracted  for  the  noncubic  contributions,  the  conclusion  of  JT  activity  remains 
valid.  Previous  calculations  of  the  band  contours  of  Ag/Ar  by  O'Brien50,  using  a  (e  + 
t^  equal  frequency  JT  model  with  SO  interaction,  strongly  supports  the  present 
explanation.  The  present  and  past  work  on  the  noble  metals  taken  together  with  previous 
work  on  the  alkali  metals  in  matrices  present  conclusive  evidence  for  significant  matrix- 
mediated  effects  in  the  form  of  Jahn-Teller  activity  and  enhanced  or  reduced  spin-orbit 
coupling. 


CHAPTER  V 
MOLYBDENUM  ATOMS  IN  RARE  GAS  MATRICES 


Introduction 

Moments  analysis  of  MCD  and  absorption  spectra  have  successfully  shown  that 
matrix  mediated  SO  couplings  and  the  Jahn-Teller  effect  are  responsible  for  the 
appearance  of  the  bandshapes  of  matrix  isolated  atoms  which  possess  a  2P  «-  2S 
transition.23-26-21'22  These  effects  are  not  limited  to  2P  *-  2S  transitions  as  evidenced  by 
the  MCD  studies  of  the  6P  *-  6S  transition  of  Re  atoms51  and  the  7P  *-  7S  transition  of 
Cr  atoms.    •       Molybdenum  atoms,  belonging  to  the  same  group  as  Cr,  also  possess 
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several    P  *-    S  transitions  that  may  be  influenced  by  these  matrix  mediated  effects. 

The  absorption  spectra  of  matrix  isolated  Mo  atoms  has  been  reported 
previously.  ^  ■       Band  assignments  were  made  by  comparison  with  gas  phase  spectra. 

S7  7  7 

Green  and  Gruen  found  no  unexpected  bands,  as  seen  in  the  P  *-  S  cases,  and 
concluded  that  no  crystal  field  effects  were  necessary  to  explain  the  spectra.   They  also 

S7  7 

concluded  that  "the  effects  of  the  matrices  on  the  z  P  sublevel  frequency  splittings  and 
relative  intensities  are  not  simple  interactions  and  appear  to  involve  a  change  in  the 
configuration  interaction  as  well  as  the  spin-orbit  coupling." 

Using  a  moments  analysis  of  MCD  and  absorption  spectra  this  investigation  will 
show   that   matrix   influenced   SO   coupling  and  the   Jahn-Teller  effect  are   indeed 
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responsible  for  the  changes  in  the  z  band  system  noted  by  Green  and  Gruen.  Second- 
order  SO  interactions  will  also  be  considered  as  a  source  of  intensity  for  forbidden 
transitions  observed  in  some  of  the  present  spectra  and  also  by  Green  and  Gruen. 
Finally,  an  extension  of  the  theory  proposed  by  Pellow  and  Vala  to  explain  the  SO 
coupling  in  matrix  isolated  metal  atoms  with  a  'P  *-  ~S  transition  will  be  used  to  predict 
the  SO  reduction  factors  for  matrix  isolated  Mo  atoms. 

Experimental 

Molybdenum  metal  was  vaporized  using  a  modification  of  the  resistive  heating 
technique  used  for  Ag  and  Au  atoms.  To  produce  a  vapor  pressure  of  -10  torr  Mo 
must  be  heated  to  near  2000°  C.  Heating  a  Knudsen  cell  to  this  temperature  would 
require  a  current  of  approximately  400  amps  with  the  furnace  used  in  this  work.  This  is 
beyond  the  current  handling  capability  of  the  cables  connecting  the  furnace  to  the  variac. 
To  overcome  this  problem  several  strands  of  Mo  wire  (2mm  diameter,  Aesar)  were 
twisted  together  and  fastened  between  the  electrodes  of  the  furnace.  With  the  lower 
resistance  only  a  few  amps  were  required  to  heat  the  wires  to  between  2000  and  2050°C. 
Great  care  was  necessary  to  prevent  passing  too  much  current  through  the  wires  and 
melting  them.  Even  with  care  the  wires  generally  melted  through  after  three  or  four 
depositions. 

Moments  were  calculated  from  the  spectra  as  recorded.  The  MCD  baseline 
generally  sloped  upward  as  the  spectrum  was  scanned  toward  higher  energy  but  the  CD 
baseline  also  followed  this  pattern  and  showed  no  underlying  bands.    Since  the  bands 
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were  not  completely  resolved  and  no  band  was  separated  from  the  other  components,  it 

was  not  possible  to  use  the  metal  as  a  temperature  probe  of  its  own  environment  as  was 

done  in  the  Ag  and  Au  cases.    An  average  temperature  correction  (-0.5  K)  for  the 

Displex  was  calculated  from  the  corrections  recorded  for  Au  atoms.    This  was  used 

throughout  the  moments  analysis.   Calibrations  of  O.D.  (absorption)  and  AO.D.  (MCD) 

were  carried  out  by  the  method  described  for  Ag  and  Au  atoms  using  the  two  standard 

solutions. 

Results 

The  spectra  in  all  three  matrices  show  two  distinct  band  structures:  a  resolved 
triplet  at  lower  energy  and  a  unresolved  broad  band  to  higher  energy.  These  are  shown 
in  Figures  5.1-5.3  with  band  positions  and  assignments  listed  in  Table  5.1.  These  bands 
correlate  well  with  those  reported  by  Green  and  Gruen  although  the  positions  reported 
here  are  consistently  at  slightly  lower  energy  than  those  previously  reported.  Based  on 
gas  phase  data20  the  lower  energy  band  is  assigned  to  the  z  7P  *-  7S  transition  while  the 
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higher  energy  band  is  assigned  to  the  y  P  *-  S  transition.  The  MCD  temperature 
dependence  of  the  z  7P  *-  7S  band  for  Mo/Ar  is  shown  in  Fig.  5.4.  The  other  matrices 
show  similar  temperature  dependence  of  the  MCD.  The  spacing  of  the  sublevel  energies 
of  the  z  band  increases  as  the  matrix  gas  is  changed  as  expected  from  the  previous 
work.  In  the  gas  phase  the  components  appear  at  26321,  25872  and  25614  cm  with 
spacings  of  449  and  258  cm-1.  In  all  three  rare  gases  the  spacings  are  increased  with  the 
largest  increase  seen  in  Ar  (666  and  549  cm")  and  the  least  in  Xe  (456  and  419  cm"  ). 


Figure  5.1        Experimental  MCD  (top)  and  absorption  (bottom)  spectra  of  Mo  atoms  in 
a  Ar  matrix  at  10.3  K. 
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Figure  5.2        Experimental  MCD  (top)  and  absorption  (bottom)  spectra  of  Mo  atoms  in 
a  Kr  matrix  at  10.9  K. 
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Figure  5.3        Experimental  MCD  (top)  and  absorption  (bottom)  spectra  of  Mo  atoms  in 
a  Xe  matrix  at  11.3  K. 
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Table  5.1     Observed  Absorption  Band  Positions  (in  cm  l)  and  Parameters  Derived  from 
the  Moments  Analysis  for  Mo  Atoms  in  Rare  Gas  Matrices. 


Gasa 

Ar 

Kr 

Xe 

Excited  State 

z7P2 

25614 

27885 

27317 

26492 

z7P3 

25872 

28479 

27807 

26911 

z?P4 

26321 

29145 

28375 

27367 

y7P2 

31300 
31533 

34114(sh) 
34479 

33365(sh) 
33625 

32268 
32563 

y7p4 

31913 

33942 

32930 

^/cm"1 

101 

126 

84.5 

57.5 

<E2>N°c/cnT2 

- 

1.9X105 

1.7x10s 

1.8xl05 

<E2><9  /cm"2 

- 

5700 

19900 

5060 

v^/cm"1 

- 

44.7 

47.7 

43.2 

aValues  from  Ref.  20. 

bFirst  order,  obtained  from  (28). 

cValues  from  fitting  (33). 


97 

8- 

a    —    11.3  K 

6- 

Mo/Ar          a!  —   13.5  K 

IaI  —  17-5  k 

4- 

CM 

1 

o 

-     2  - 

i— 
< 

<   o- 

-2- 

-4- 

-6  - 

26000   27000   28000   29000   30000   31000 

FREQUENCY  /cm"1 

Figure  5.4        MCD  temperature  dependence  for  the  z  7P  *-  7S  transition  of  Mo  atoms 
in  a  At  matrix. 
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Green  and  Gruen     also  report  evidence  of  several  forbidden  transitions  which  are 

not  readily  observed  in  the  present  spectra.    Table  5.2  shows  the  reported  gas  phase 

energies  for  the  observed,  allowed  transitions  as  well  as  several  symmetry  or  orbitally 

forbidden  transitions.  A  small  negative  MCD  band  at  33312  cm-1  for  Mo/Ar  may  result 

7  7 

from  forbidden  D  «-  S  transition  which  would  be  predicted  to  have  a  negative  C0  term. 
A  more  concentrated  sample  of  Mo/Kr  was  made  which  shows  (Fig.  5.5)  several  small 
MCD  bands  between  the  two  allowed  transitions  which  correlate  with  the  forbidden  bands 
reported  by  Green  and  Gruen.  A  similar  situation  was  reported  by  Vala  et  al.  in  the 
study  of  Cr  atoms  in  Kr  and  Xe  matrices  where  a  moments  analysis  showed  that  isolation 
in  the  matrix  caused  a  considerable  enhancement  of  second-order  SO  coupling.    A 

7  7  5 

orbitally  forbidden  D  state  was  found  to  mix  with  the  z  P  state  and  a  spin  forbidden  P 
state  was  determined  to  mix  with  the  y  7P  state.  As  Table  5.2  shows  similar  forbidden 
states  lie  near  the  allowed  states  for  Mo  atoms. 

7  7 

Moments  analysis  were  performed  on  the  z  P  *-  S  transition  for  all  three  rare 
gases.  The  SO  coupling  constant,  £,  was  extracted  from  plots  of  <AA>1/<A>0  vs.  T"1 
which  are  shown  in  Figs.  5.6.  In  the  gas  phase20  the  total  splitting  is  707  cm-1  which 
is  related  to  the  SO  coupling  constant  by: 

A  =(2S+1)C-  (62) 

In  At,  the  SO  coupling  constant,  t»  was  found  to  increase  from  the  gas  phase  value  (101 
cm-1)  to  126  cm-1  while  in  Kr  and  Xe  matrices  it  was  found  to  decrease  to  84.5  and 
57.5  cm-1,  respectively.  The  decrease  in  t,  in  going  from  Ar  to  Kr  to  Xe  matrices  is  the 
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Table  5.2     Gas  Phase  Absorption  Band  Positions  and  gf  Values  for  Allowed  and  Spin- 
Forbidden  or  Orbitallv  Forbidden  Transitions  for  Mo  Atoms.21,13 


v  (cm  *) 

X  (nm) 

gfc 

Excited  State 

z7P2 

25614 

390.3 

0.27 

z7P3 

25872 

386.4 

0.77 

*X 

26321 

379.8 

0.94 

z5P2 

28837 

346.7 

0.018 

z% 

28924 

345.6 

0.056 

y7D2 

31155 

320.9 

0.26 

y7P2 

31300 

319.4 

0.68 

y7P3 

31533 

317.0 

0.82 

y?D3 

31655 

315.8 

0.57 

y7P4 

31913 

313.3 

1.4 

y7o4 

32123 

311.2 

0.14 

aValues  from  Ref.  20. 
bValues  from  Ref.  55. 
cgf  (g:g  factor;  f:oscillator  strength) 
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Figure  5.5  Expansion  of  MCD  from  a  concentrated  Mo/Kr  matrix.  Bands  marked 
with  a  *  may  be  due  to  forbidden  transitions.  (+)  band  for  5P  *-  7S  and 
(-)  band  for  7D  *-  7S.   No  absorption  was  visible  for  either  band. 
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Figure  5.6        Reduced  MCD  first  moment  plots  vs.  inverse  temperature  for  the  three 
Mo/RG  systems. 
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same  trend  as  found  for  the  Cu21'22,  Ag,  and  Au  atoms  as  well  as  for  Cr39,47  atoms. 

7  7 

Calculations  for  the  y    P  *-   S  transition  were  not  performed. 

The  presence  of  second-order  SO  coupling  was  again  assessed  by  plotting  the 
reduced  second  MCD  moment  vs.  T"1  (eq.  41)  with  as  =  2S(S+l)/3  =  8.  In  all  cases  the 
plots  were  too  scattered  to  give  a  reliable  fit.  This  unfortunately  also  prevents  a 
quantitative  determination  of  mixing  coefficients.  A  second  method  for  determining  the 
presence  of  second-order  SO  effects  involves  the  zeroth  MCD  moment,  <AA>0.  In 
Au/Ar,  where  second  order  effects  were  quite  large,  <AA>0  was  found  to  show  a 
temperature  dependence  (i.e.  linear  with  1/T).  In  the  Mo/RG  cases  plots  of  <AA>0  vs. 
1/T  had  essentially  zero  slope  although  there  was  considerable  scatter.  This  leads  to  the 
conclusion  that  matrix  enhancement  of  second-order  SO  effects  is  negligible. 

Some  comparisons  between  the  Cr  and  Mo  spectra  also  point  to  little  or  no 
enhancement  of  second-order  SO  splitting  by  the  matrix.  In  the  Cr  atom  studies,39 
forbidden  transitions  were  seen  quite  strongly  even  though  their  gf  (g:  g  factor,  f: 
oscillator  strength)  values  are  two  orders  of  magnitude  lower  than  the  gf  values  of  the 
allowed  transitions.  For  Mo,  these  transitions  are  not  readily  seen  even  though  the 
forbidden  gf  values  (Table  5.2)  are  only  an  order  of  magnitude  less  than  the  allowed 
transitions.  This  indicates  that  isolation  in  the  matrix  is  not  causing  any  significant 
donation  of  oscillator  strength  to  the  forbidden  transitions. 

The  presence  of  substantial  non-cubic  contributions  to  the  observed  bandwidth 
was  once  again  used  as  a  test  for  Jahn-Teller  activity  in  the  excited  state.  For  Mo  atoms 
equations  (58)  and  (59)  become 
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R  =  <A>7/<A>n  =  <E2>r  +  <E2>Nr  +  8£2  <63) 


P  =  <AA>3/<AA>X  =  3<E2>C  +  3<E2>NC/2  +  12£2  (64) 

which  may  be  rearranged  to  give  expressions  for  <E  >c  and  <E  >NC: 

<E2>NC  =  2R-2P/3-8£2  (65) 


<E2>C  =  2P/3-R.  (66) 

The  calculated  values  are  shown  plotted  vs.  T  in  Figs.  5.7-5.9.  The  <E  >NC  curves  were 
fit  using  (33)  and  values  for  <E2>N°C  and  v;  are  listed  in  Table  5.1.  This  confirms  that 
the  Jahn-Teller  effect  is  active  in  these  matrix  isolated  species  since  the  non-cubic 
contributions  are  substantially  greater  than  the  cubic  contributions  and  show  the  expected 
hyperbolic  cotangent  dependence  on  temperature. 

Discussion 

The  SO  splitting  for  Mo  atoms  in  matrices  has  been  shown  to  be  greatly 
influenced  by  the  matrix  environment  with  an  increase  over  the  gas  phases  value 
occurring  in  Ar  and  a  decrease  in  Kr  and  Xe  matrices.  In  none  of  the  matrices  has  the 
SO  coupling  constant  become  negative  which  would  indicate  an  inversion  of  the  sublevels 
of  the  excited  state.  This  effect  was  noted  for  Cr  atoms  in  both  Kr  and  Xe  matrices.39 
Isolation  of  the  atoms  in  the  rare  gas  matrices  also  does  not  appear  to  be  enhancing 
second-order  SO  coupling. 
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Figure  5.7  Temperature  dependence  of  the  cubic  <E  >c  and  noncubic  <E  >NC 
bandwidth  contributions  of  Mo/Ar.  Solid  line  fit  follows  the  relation: 
<E2>NC  =  1.9X105  coth(44.7/2kT),  in  cm-1. 


105 


30 

Mo/Kr 

25 

T      20 

E 

0 

+      15 
o 

>  io 

UJ 

v      5 

o  cP^^ 

Sb,T^^ 

o   *c 

0 

<#«%><&#><>  o 

-5 

1,1,!, 

0          10         20         30         40 

T/K 

Figure  5.8        Temperature  dependence  of  the  cubic  <E  >c  and  noncubic  <E  >NC 
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Figure  5.9  Temperature  dependence  of  the  cubic  <E  >c  and  noncubic  <E  >NC 
bandwidth  contributions  of  Mo/Xe.  Solid  line  fit  follows  the  relation: 
<E2>NC  =  1.8x10s  coth(43.2/2kT),  in  cm-1. 
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The  Jahn-Teller  effect  was  also  shown  to  be  strongly  active  in  all  three  matrices 

as  evidenced  by  large  contributions  to  the  observed  bandwidths  from  non-cubic  lattice 

modes.    This  activity  along  with  the  SO  splitting  is  responsible  for  the  resolved  bands 

found  for  the  z    P  *-    S  transition.    This  parallels  the  case  of  Cr  atoms  in  matrices. 

Although  no  analysis  was  performed  on  the  y  band  system  it  may  be  expected  that  the 

Jahn-Teller  effect  is  less  intense  (or  absent)  in  this  transition.    Similar  results  were 

observed  for  Cr  atoms39  where  the  y  band  system  remained  unresolved  and  little  or  no 

Jahn-Teller  activity  was  found. 

Concurrently  with  the  recording  of  these  spectra  predictions  were  made  for  the  SO 

reduction  factors.    The  application  of  a  theory     derived  for  a  one-electron    P  *-    S 

transition  to  a  six-electron  7P  *-  7S  transition  may  seem  inappropriate  at  first.  However, 

it  can  be  shown     that  since  the  five  unpaired  electrons  in  the  4d  shell  of  Mo  are  coupled 

to  a   S  state  the  SO  coupling  constant  depends  only  on  the  one  electron  5p  SO  constant 

as 

^7P  =  V6 
In  analogy  with  the  Cr  case  large  second-order  effects  were  expected  and  following  this 
assumption  the  energy  of  the  overlapping  metal  orbitals  was  calculated  as  a  weighted  (2J 
+  1  weighting  factor)  average  of  the  allowed   P  and  forbidden   P  transition  energies.  The 
resulting  k     and  kave  values  are  listed  in  Table  5.3. 

After  the  moments  analysis  the  assumed  second-order  effects  were  not  seen  and 
the  energy  of  the  overlapping  metal  orbitals  was  recalculated  without  contributions  from 
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the  forbidden  5P  state.   The  resulting  k's  are  also  given  in  Table  5.3.    It  is  interesting  to 

note  that  for  the  two  cases  the  1^    and  kave  values  have  almost  exactly  switched. 

Comparing  the  experimental  k  values  with  those  predicted  without  5P  influence 

shows  that  kave  closely  follows  the  observed  pattern  of  Ar>Kr>Xe.  Additionally,  the  Ar 

system  is  found  to  be  case  I  (t,  positive  and  larger  than  the  gas  phase)  while  the  Kr  and 

Xe  cases  are  found  to  be  case  II  (t,  positive  but  less  than  the  gas  phase).  For  Au  and  Ag, 

kave   was   also   found   to   best   follow   the  observed  pattern   and   magnitude   of  the 

experimental  k's.    For  Mo/Xe  the  observed  value  is  considerably  lower  than  the  value 

predicted  by  kave.     It  is  actually  closer  to  the  value  predicted  by  k^    when  the    P 

influence  is  included.    This  may  indicate  that  some  enhancement  of  second-order  SO 

coupling  is  occurring  in  the  Xe  matrix  although  the  effect  is  too  weak  to  observe 

experimentally.  Overall,  the  results  indicate  that  it  is  possible  to  apply  this  "one  electron" 

theory  to  a  many  electron  system. 
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Table  5.3         Comparison  of  Experimental  and  Theoretical  SO  Coupling  "Reduction" 
Factors  for  Molybdenum  Atoms  in  Rare  Gas  Matrices 


Rare 

Metal  Atom 

Gas 

Mo 

kno  ('P)3 

kaVe  ('P)3 

K« 

ave 

k     b 

Ar 

1.06 

1.15 

1.14 

1.05 

1.3 

Kr 

0.80 

1.08 

1.07 

0.82 

0.8 

Xe 

0.69 

1.01 

1.01 

0.72 

0.5 

Includes  mixing  with  JP  state. 
b±0.1 


CHAPTER  VI 
SPECTRA  OF  MATRIX  ISOLATED  COBALT  AND  TITANIUM  ATOMS 


Introduction 

The  studies  of  Au,  Ag  and  Mo  atoms  detailed  here  and  many  of  the  references 
cited  in  this  work  have  shown  the  power  of  moments  analysis  when  applied  to  MCD  and 
absorption  spectra  in  extracting  information  about  excited  states  that  would  otherwise  not 
be  available.  The  application  of  this  technique  to  the  study  of  transition  metal  diatomics, 
both  hetero-  and  homo-nuclear,  is  a  logical  step  in  the  extension  of  the  present  studies. 
General  expressions  for  MCD  of  randomly  oriented  linear  molecules  have  been  recently 
published  by  Pellow  and  Vala.  Since  most  matrices  which  contain  a  diatomic  species 
will  also  contain  isolated  component  atoms,  no  MCD  studies  of  matrix  isolated  diatomics 
can  proceed  until  spectra  of  these  atoms,  both  MCD  and  absorption,  are  recorded  and 
analyzed. 

While  both  hetero-  and  homo-nuclear  diatomics  of  a  number  of  transition  metals 
have  been  formed,  two  recent  reviews  by  Weltner  and  Van  Zee58  and  Morse  show  that 
the  first  row  transition  metal  diatomics  have  been  most  extensively  studied.  Of  the  first 
row  transition  metals  only  Cr39'47,  Mn60,  Fe61'62,  Ni63  and  Cu21'22  have  been  studied  in 
matrices  using  MCD.    The  studies  are  extended  here  to  include  Co  and  Ti.    Matrix- 
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Ill 

induced  energy  shifts  and  band  assignments  will  be  made  by  comparing  the  observed 
spectra  with  simulated  spectra  produced  from  gas  phase  data. 

Experimental 

Initial  investigations  were  conducted  by  vaporizing  Ti  metal  from  a  Knudsen  cell 
at  temperatures  of  approximately  1950°C  and  condensing  the  atomic  beam  with  a 
continuous  flow  of  argon.  The  currents  needed  to  raise  the  Knudsen  cell  to  this 
temperature  were  nearly  400  Amps,  well  beyond  the  rated  capacity  of  the  furnace  leads. 
These  depositions  required  several  hours  to  produce  sufficient  quantities  of  trapped  atoms 
for  study  and  required  constant  monitoring  to  prevent  overheating.  The  long  depositions 
often  resulted  in  cloudy  matrices  due  to  the  amount  of  gas  deposited.  The  availability  of 
a  Nd:YAG  laser  prompted  the  construction  of  the  laser  ablation  source  described  in 
Chapter  III.  Using  the  continuous  flow  of  argon  for  the  matrix  the  output  of  the  Q- 
switched  laser,  operating  at  532  nm  and  10  Hz,  was  focused  onto  a  target  rod  of  Ti.  The 
rod  was  continuously  rotated  to  prevent  the  laser  from  drilling  deep  pits  into  the  rod. 
Deposition  times  ranged  from  15  minutes  to  one  hour.  In  these  attempts  multiple  matrix 
sites  were  always  present,  many  of  which  could  not  be  annealed  out  without  seriously 
degrading  the  clarity  of  the  matrix.  These  matrices  provided  little  chance  for  assignment 
of  bands. 

A  more  recent  construction  of  the  source  involves  the  use  of  a  pulsed  valve  in  the 
gas  inlet  line.  Careful  timing  of  the  gas  pulse  and  laser  firing  will  successfully  entrain 
the  metal  in  the  gas  and  carry  it  onto  the  matrix  window.    More  details  of  this  apparatus 
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will  be  presented  in  the  following  chapter.  Depositions  with  this  source,  still  using  argon 

as  the  matrix  gas,  produced  clear  matrices  which  match  those  produced  by  the  furnace 

with  little  or  no  evidence  of  multiple  sites.   These  spectra  will  be  analyzed  here. 

The  investigations  of  Co  atoms  also  began  with  a  Knudsen  cell  technique. 
Unfortunately,  at  high  temperatures  Co  alloys  with  the  Ta  Knudsen  cells  and  causes  them 
to  break  open  which  causes  a  interruption  of  the  current  flow  and  ends  the  deposition. 
A  ceramic  liner  may  be  used  to  prevent  this  problem  but  a  large  amount  of  current  was 
again  required  to  reach  the  temperatures  necessary  to  produce  the  required  vapor  pressure 
of  Co.  The  laser  ablation  source  was  used  to  produce  the  spectra  described  here  but  these 
experiments  were  completed  before  the  pulsed  valve  was  introduced  to  the  system. 
Consequently,  the  Ar  and  Kr  matrices  show  evidence  of  multiple  sites. 

To  aid  in  band  assignments  simulated  absorption  and  MCD  spectra  were  produced 
in  the  following  manner.  A  listing  of  gas  phase  transitions  observed  from  the  electronic 
ground  state  and  their  wavelengths  was  made.  The  absolute  intensities  (i.e.,  gf  values) 
for  each  of  these  transitions  were  also  noted.  The  simulated  absorption  spectra  were 
produced  using  the  gas  phase  wavelengths  as  band  positions,  the  gf  values  as  intensities 
and  Gaussian  bandshapes.  The  bandwidths  of  the  simulation  were  varied  to  produce  a 
spectrum  that  closely  approximates  the  bandwidths  of  the  experimental  spectrum.  In 
some  cases  narrow  bandwidth  simulations  were  used  to  determine  area  where 
experimental  bandwidths  may  be  concealing  transitions.  The  MCD  spectra  were 
simulated  in  a  similar  manner  but  the  relative  intensities  and  signs  of  the  bands  were 
produced  by  multiplying  the  gf  values  by  the  theoretical  Cq/©0  values.    These  latter 
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theoretical  values  were  calculated  from  general  formulas61,64  which  are  shown  in  Table 

6.1.  By  comparing  these  simulations  with  the  observed  spectra  assignments  can  be  made 

and  matrix  induced  energy  shifts  can  be  assessed.    The  absorption  spectra  are  usually 

sufficient  to  assign  the  bands  in  cases  where  few  transitions  are  observed.  In  cases  where 

a  large  number  of  bands  are  observed,  the  sign  and  pattern  of  the  MCD  bands 

corresponding  to  the  absorptions  often  provides  a  more  confident  assignment  than  the 

absorption  spectrum  alone.    In  the  Co  case,  the  MCD  data  will  be  used  to  show  that 

several  gas  phase  energy  levels  have  been  misassigned. 

Titanium 

The  gas  phase  data  used  in  the  simulations  of  absorption  and  MCD  spectra  are 
listed  in  Table  6.2.  A  large  number  of  transitions  from  the  F2  ground  state  have  been 
observed  in  the  gas  phase.  In  producing  the  simulated  spectra  the  bandwidths  of  the 
Gaussian  bands  were  adjusted  to  closely  match  those  observed  in  the  experimental 
spectra.  This  procedure  has  obscured  many  of  the  gas  phase  transitions  due  to  band 
overlap.  The  bands  responsible  for  the  majority  of  intensity  in  the  simulations  have  been 
marked  with  an  asterisk. 

Titanium  Atoms  in  Argon  Matrices 

The  initial  laser  experiments  produced  matrices  with  a  number  of  multiple  sites 
which  made  conclusive  assignments  difficult.  The  appearance  of  these  spectra  is 
compared  with  those   from   the  pulsed  valve/laser  vaporization  setup   in  Figs.   6.1 
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(absorption)  and  6.2  (MCD).  The  simplification  of  the  spectra  by  the  appearance  of  only 

one  site  will  be  shown  to  make  band  assignments  particularly  straightforward  in  this  case. 

The  absorption  spectrum  of  Ti/Ar  in  the  410-250  nm  range  is  shown  in  (a)  of  Fig. 
6.3.  The  lower  trace  (b)  in  the  figure  is  a  simulation  produced  using  gas  phase  data  with 
bandwidths  set  to  approximate  the  experimental  bandwidths.  Except  for  some  shifting  of 
transition  energies  to  the  blue  in  the  experimental  spectrum  the  two  are  essentially 
identical.  The  corresponding  observed  and  simulated  MCD  for  this  region  are  shown  in 
(a)  and  (b),  respectively,  of  Fig.  6.4.  Small  bands  appear  in  the  MCD  on  either  side  of 
several  of  the  major  bands.  These  may  result  from  a  matrix  induced  shift  of  one  of  the 
many  underlying  Ti  atom  bands,  from  a  small  amount  of  secondary  site  or  from  inpurities 
in  the  sample  rod.  The  close  match  between  simulated  and  observed  spectra,  both 
absorption  and  MCD,  allow  conclusive  assignment  of  each  transition  observed  in  the 
matrix  isolated  sample.  These  are  listed  in  Table  6.3  and  parallel  those  published  by 
Gruen  and  Carstens65  based  on  their  study  of  absorption  spectra. 

The  region  to  longer  wavelength  might  be  expected  to  show  several  bands  based 
on  the  gas  phase  observations  but  no  absorptions  were  noted  in  most  matrices.  In  several 
very  concentrated  matrices  two  broad  absorptions  could  be  identified  at  484.2  and  443.4 
nm.  The  positive  and  negative  MCD,  respectively,  associated  with  these  absorptions 
allows  them  to  be  assigned  to  Dj  (gas  phase  501.6  nm)  and  3G3  transitions  (gas  phase 
465.8  nm). 

The  assignment  of  bands  in  this  case  was  simplified  by  the  appearance  of  only  one 
site  and  the  observation  that  the  simulated  and  observed  spectra,  both  absorption  and 


Table  6.1  Theoretical  CqADq  (MCD)  Terms  for  the  Atomic  Transitions  J'  *-  J. 
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J' 

Cq/Vq  (MCD)a 

J  + 1 

J 
J  -  1 

"JgA 

Sa 
(J  +  l)gA 

a  Ref.  64. 


Table  6.2     Gas  Phase  Data  for  Titanium  Atoms3 
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Ti  Atoms 

F2  ground  state 

excited  state 

Xgas  (nm) 

VMS  (cm_1) 

gfb 

gf(Co/2>o) 

*3F2 

398.3 

25107 

0.82 

0.541 

3F3 

396.4 

25227 

0.17 

-0.224 

*3D, 

395.0 

25316 

0.81 

1.60 

3D2 

393.1 

25439 

0.16 

0.106 

3P2 

392.2 

25497 

0.052 

0.034 

*1 

391.6 

25536 

0.0044 

0.0087 

3D3 

390.0 

25641 

0.016 

-0.021 

5D2 

389.1 

25700 

0.014 

0.0092 

5D3 

387.6 

25800 

0.051 

-0.067 

*3F2 

373.1 

26802 

0.65 

0.429 

3F3 

371.8 

26896 

0.10 

-0.132 

^ 

365.6 

27352 

0.072 

0.143 

3D2 

364.7 

27420 

0.045 

0.029 

3D3 

363.9 

27480 

0.031 

-0.040 

*3G3 

363.6 

27503 

1.2 

-1.56 

**. 

337.1 

29665 

0.43 

0.85 

3D2 

335.9 

29771 

0.11 

0.073 

3D3 

334.3 

29913 

0.047 

-0.062 

*3G3 

334.3 

29913 

2.3 

-3.04 

*3G3 

318.7 

31377 

1.2 

-1.58 

agas  phase  data  from  Ref.  20. 
bgf  values  from  Ref.  55. 


Table  6.2  —  continued 
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Ti  Atoms 

F2  ground  state 

excited  state 

Xgls  (nm) 

Vs  (cm'1} 

gfb 

gf(Co/2>o) 

3F2 

297.1 

33659 

0.066 

0.04 

3F3 

296.9 

33681 

0.024 

-0.032 

*3F2 

294.3 

33979 

1.0 

0.66 

3F3 

293.4 

34083 

0.18 

-0.24 

3G3 

3 

266.3 

37552 

0.16 

-0.21 

r3 

265.8 

37622 

0.073 

-0.096 

3F2 

265.6 

37651 

0.029 

0.019 

«* 

264.2 

37850 

1.3 

2.57 

3D2 

263.3 

37979 

0.36 

0.238 

*3F2 

260.1 

38447 

0.69 

0.455 

*3F3 

259.4 

38550 

0.50 

-0.66 

♦3D, 

252.1 

39667 

0.47 

0.93 

3D2 

251.8 

39714 

0.13 

0.086 

agas  phase  data  from  Ref.  20. 
bgf  values  from  Ref.  55. 


Figure  6.1        Absorption  spectra  of  Ti  atoms  in  an  Ar  matrix. 

(a)  Multiple  sites  from  laser  vaporization  with  continuous  gas  flow. 

(b)  Single  site  from  laser  vaporization  with  pulsed  gas. 


Figure  6.2        MCD  spectra  of  Ti  atoms  in  an  Ar  matrix. 

(a)  Multiple  sites  from  laser  vaporization  with  continuous  gas  flow. 

(b)  Single  site  from  laser  vaporization  with  pulsed  gas. 
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Figure  6.3       (a)  Absorption  spectrum  of  Ti  atoms  in  an  Ar  matrix,  410-250  nm. 

(b)  Simulated  absorption  spectrum  produced  from  gas  phase  band  positions 
and  intensities. 
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Figure  6.4        (a)  MCD  spectrum  of  Ti  atoms  in  an  Ar  matrix,  410-250  nm. 

(b)  Simulated  MCD  spectrum  produced  from  gas  phase  band  positions  and 
intensities. 


Table  6.3    Assignments  for  Titanium  Atoms  in  an  Ar  Matrix 
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Ti  Atoms  in  an  At  Matrix       3F,  ground  state 

excited  state* 

Kin™) 

vm  (cm"1) 

Av  (cm-1)  b 

obs.  MCD 

*3F2 

3891 

25700 

593 

+ 

3F3 

*% 

3843 

26021 

705 

+ 

*2 

3P2 

*i 

3D3 

5D2 

% 

*% 

3634 

27533 

731 

+ 

% 

^ 

3D2 

3D3 

*3G3 

3561 

28082 

579 

- 

*3Dl 

3303 

30276 

611 

+ 

3D2 

3D3 

*3G3 

3269 

30590 

677 

- 

*3G3 

3081 

32457 

1080 

+ 

agas  phase  data  from  Ref.  20. 

Av  =  v    -  v 

ra         gas 


Table  6.3  —  continued 
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Ti  Atoms  in  an  Ar  Matrix 

Ft  ground  state 

excited  statea 

^W 

vm  (cm-1} 

Av  (cm-1)  b 

obs.  MCD 

3F2 

3F3 

*% 

2872 

34819 

840 

+ 

% 

% 

% 

3F2 

**! 

2555 

39139 

1289 

- 

3D2 

*3F2 

*% 

JO, 

*. 

agas  phase  data  from  Ref.  20. 

Av  =  v     -  v 

m        gas 
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MCD,  appear  almost  identical  except  for  energy  shifts  in  the  matrix  spectra.   While  the 

MCD  spectra  served  primarily  as  a  confirmation  of  assignments  in  this  case,  the  following 

section  will  show  the  MCD  spectra  as  a  severe  test  of  gas  phase  band  assignments. 

Cobalt 

For  the  initial  investigations  Co  atoms  were  isolated  in  an  Ar  matrix.  The 
recorded  absorption  spectra  agreed  well  with  those  of  Ozin  and  Hanlan.  Absolute 
assignments  were  complicated  by  the  presence  of  several  stable  sites  within  the  matrix. 
The  spectra  of  the  other  group  VIIIB  metals  in  argon  matrices  were  also  complicated  by 
numerous  sites  but  the  spectra  recorded  in  krypton  matrices  were  free  from  influence  of 
other  sites.  The  spectra  of  Co  atoms  in  Kr  matrices  were  recorded  and  once  more 

showed  a  small  amount  of  secondary  site  that  remained  stable  to  prolonged  annealing. 
The  desire  to  obtain  a  single  site  matrix  led  to  a  deposition  using  Xe  as  the  matrix 
material.  These  matrices  appear  to  produce  only  one  site  and  will  be  used  to  establish 
the  band  assignments  before  the  complications  of  overlapping  bands  in  the  Kr  and  Ar 
matrices  are  addressed.  Table  6.4  lists  the  gas  phase  wavelengths,  frequencies  and  gf 
values  used  in  the  simulation  of  Co  spectra. 

Cobalt  Atoms  in  Xenon  Matrices 

The  assignment  of  transitions  for  the  Xe  matrix  will  be  pursued  first  due  to 
apparent  lack  of  matrix  sites.  For  convenience  the  spectrum  has  been  divided  into  two 
sections:  380-280  nm  and  270-220  nm.    No  bands  were  observed  between  the  two 


Table  6.4      Gas  Phase  Data  for  Cobalt  Atoms3 
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Co  Atoms         F9/2  ground  state      [3d    4s  ] 

excited  state 

\as  (nm) 

veas  (cm_1) 

gfb 

configuration 

gK<V2>o) 

%n 

217.4 

45971 

.081 

3d8(3P)4p 

+0.593 

4r 

ull/2 

227.5 

43952 

.042 

3d7(a5P)4p 

-0.252 

2F 

r7/2 

229.5 

43555 

.048 

" 

+0.351 

4F 

r9/2 

230.9 

43295 

.54 

ti 

+0.718 

u9/2 

236.5 

42269 

.27 

3d74s(b3F)4p 

+0.359 

4F 

238.4 

41918 

.47 

if 

+3.44 

4r 

ull/2 

240.7 

41529 

2  2 

M 

-13.18 

4F 

r9/2 

242.5 

41226 

1.9 

M 

+2.53 

'7/2 

243.6 

41041 

.35 

It 

+2.56 

4D7/2 

252.1 

39649 

1.9 

II 

+13.91 

2<37/2 

292.9 

34134 

.0056 

3d8(a3F)4p 

+0.041 

4F 

r7/2 

298.7 

33467 

.11 

11 

+0.805 

2<~ 

u9/2 

298.9 

33440 

.11 

11 

+0.146 

4°7/2 

301.4 

33173 

.049 

» 

+0.359 

4F 

r9/2 

304.4 

32842 

.44 

ll 

+0.585 

uii/2 

308.3 

32431 

.086 

11 

-0.515 

4°7/2 

312.1 

32028 

.040 

M 

+0.293 

%n 

341.3 

29295 

.20 

3d74s(a5F)4p 

+  1.46 

4r 

ull/2 

346.6 

28845 

.42 

M 

-2.52 

4F 
r7/2 

347.4 

28777 

.64 

II 

+4.68 

4F 
r9^ 

352.7 

28346 

.48 

II 

+0.638 

agas  phase  data  from  Ref.  20. 
bgf  values  from  Ref.  55. 
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regions  or  to  longer  wavelength  and  limitations  of  the  optics  prevented  scanning  to  shorter 

wavelengths. 

The  absorption  spectrum  for  the  380-280  ran  region  is  shown  in  (a)  of  Fig.  6.5 
with  the  corresponding  simulation  shown  in  (b).  In  this  region  the  long  wavelength  group 
of  bands  appears  blue  shifted  while  the  shorter  wavelength  group  appears  red  shifted. 
Overall,  the  observed  spectra  agree  well  with  the  simulation.  Only  the  simulated  band 
at  347.4  nm  is  particularly  stronger  than  the  corresponding  experimental  transition.  This 
can  be  easily  resolved  by  noting  that  the  simulated  band  is  the  summation  of  two 
overlapping  bands  (4F7/2:  gf=0.64  and  4G11/2:  gf=0.42)  with  similar  predicted  intensity. 
These  two  transitions  can  be  seen  in  the  observed  spectrum  as  the  peak  at  344.1  nm  and 
the  shoulder  at  342.6  nm.  The  observed  spectra  compare  well  with  those  of  Ozin  and 
Hanlan  which  leads  to  the  conclusion  that  assignment  of  the  bands  will  be 
straightforward. 

Inspection  of  the  observed  and  simulated  MCD,  Fig.  6.6  (a)  and  (b),  respectively, 
shows  that  band  assignments  may  not  be  as  straightforward  as  inspection  of  the 
absorption  spectra  indicated.  While  the  simulated  short  wavelength  region  (330-320  nm) 
matches  reasonably  well  with  the  observed  spectrum,  the  simulation  of  the  longer 
wavelength  region  (380-330  nm)  does  not  correspond  to  the  observed  MCD  pattern.  In 
this  region  four  bands  are  predicted  with  a  strong  positive  band  from  the  4F7/2  state 
dominating.  A  moderately  strong  negative  band  from  the  4G11/2  transition  is  expected  to 
lie  just  to  the  blue  of  the  positive  band.  The  observed  spectra  certainly  do  not  show  the 
pattern  predicted  by  the  simulation.    No  evidence  of  any  positive  band  from  the    F7/2 
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Figure  6.5       (a)  Absorption  spectrum  of  Co  atoms  in  a  Xe  matrix,  380-280  nm. 

(b)  Simulated  absorption  spectrum  produced  from  gas  phase  band  positions 
and  intensities. 
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Figure  6.6        (a)  MCD  spectrum  of  Co  atoms  in  a  Xe  matrix,  380-280  nm. 

(b)  Simulated  MCD  spectrum  produced  from  gas  phase  band  positions  and 
intensities. 
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state  is  seen  although  a  small  shoulder  is  noted  on  the  negative  band  corresponding  to  the 

G11/2  transition.    This  is  probably  indicative  of  a  small  amount  of  secondary  site. 

Another  anomalous  result  occurs   in  the  270-220  nm   spectral   range.      The 

experimental  and  simulated  absorptions,  shown  in  Fig.  6.7  (a)  and  (b),  respectively,  are 

again  in  excellent  agreement  with  all  matrix  shifts  being  to  the  red  except  for  the  band 

at  230.0  nm  which  is  blue  shifted.   Comparison  of  the  observed  and  simulated  MCD  for 

this  region,  Fig.  6.8  (a)  and  (b),  respectively,  again  reveals  a  discrepancy.    A  positive 

band  is  predicted  to  occur  at  243  nm  which  results  from  the  overlap  of  two  gas  phase 

bands.   The  first  of  these  bands  occurs  at  242.5  nm  in  the  gas  phase  and  results  from  a 

transition  to  a  J  =  7/2  level  of  some  unidentified  state."     The  other  band  is  a  transition 

to  a  4F9/2  level  which  occurs  at  243.6  nm.   A  strong  negative  band  is  predicted  to  lie  to 

the  blue  of  this  band  at  240.7  nm.     This  strong  negative  band  is  observed  in  the 

experimental  spectrum  at  246.8  nm  but  no  evidence  for  the  predicted  positive  band  on  its 

red  edge  is  observed.  Also,  the  observed  negative  band  is  quite  broad  which  may  indicate 

the  overlapping  of  two  negative  bands.   This  would  indicate  that  one  or  both  of  the  two 

bands  responsible  for  the  predicted  positive  peak  may  be  misassigned.   A  listing  of  the 

observed  bands  and  their  assignments  based  on  the  absorption  data  is  given  in  Table  6.5. 

In  the  discussion  section  an  explanation  for  the  differences  in  observed  and  simulated 

MCD  band  structure  will  be  pursued.  The  reassignment  of  certain  transitions  will  result 

from  this  discussion  and  these  are  shown  in  the  table  in  parenthesis. 


Figure  6.7        (a)  Absorption  spectrum  of  Co  atoms  in  a  Xe  matrix,  270-220  nm. 

(b)  Simulated  absorption  spectrum  produced  from  gas  phase  band  positions 
and  intensities. 
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Figure  6.8        (a)  MCD  spectrum  of  Co  atoms  in  a  Xe  matrix,  270-220  nm. 

(b)  Simulated  MCD  spectrum  produced  from  gas  phase  band  positions  and 
intensities. 
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Table  6.5      Band  Assignments,  Matrix  Shifts  and  Observed  MCD  Sign  for  Co  Atoms  in  a  Xenon 

Matrix 


Co  Atoms  in  Xe  matrix 

4F9/^  ground  state 

excited  statea'b 

\n  (nm) 

vm  (cm"1) 

Avc 

configuration3 

obs.  MCD 

4p 

r9/2 

230.0 

43478 

183 

3d7(a5P)4p 

+ 

4°9/2 

242.5 

41237 

-1032 

3d74s(b3F)4p 

+ 

4F 

r7/2 

243.7 

41034 

-884 

M 

+ 

4r 
ull/2 

246.8 

40519 

-1010 

Tl 

- 

4F 
r9/2 

248.1 

40306 

-920 

II 

- 

J7/2  Cll/2) 

II 

4D7/2 

257.8 

38790 

-859 

l» 

+ 

2°7/2 

3d8(a3F)4p 

4F 

^7/2 

300.4 

33289 

-189 

II 

+ 

2G9/2 

-167 

11 

4G 

u7/2 

303.4 

32960 

-213 

n 

+ 

4F 

h9/2 

307.0 

32573 

-269 

11 

+ 

4G 

ull/2 

309.8 

32279 

-152 

11 

- 

4D7/2 

314.8 

31766 

-262 

11 

+ 

4D7/2 

338.2 

29568 

273 

3d74s(a5F)4p 

+ 

4G11/2  (4F7/2) 

342.6(sh) 

29188 

343 

11 

- 

4F7/2  (4Gn/2) 

344.1 

29061 

284 

11 

-(w) 

4F 

r9/2 

349.4 

28260 

274 

11 

+ 

agas  phase  data  from  Ref.  20. 

Corrected  assignments  in  parenthesis,  see  Discussion 


cAv  =  v    -  v 
"v         ra         gas 


133 


Cobalt  Atoms  in  Krypton  Matrices 

The  experimental  absorption  spectrum,  (a),  for  the  380-280  nm  region  of  Co 
atoms  in  a  Kr  matrix  and  the  simulated  spectrum,  (b),  are  shown  in  Fig.  6.9.  The 
appearance  is  almost  identical  to  that  of  Co/Xe  for  the  same  region  except  that  the  matrix 
bands  are  shifted  to  the  blue  for  both  band  groups.  No  shoulder  appears  on  the  band  at 
336.5  nm  as  was  seen  on  the  corresponding  band  in  Xe  (338.2  nm)  but  in  Kr  the  band 
is  broadened  and  shows  no  well  defined  maximum.  The  MCD,  observed  and  simulated, 
for  this  region  is  given  in  Fig.  6.10  (a)  and  (b),  respectively.  The  observed  spectrum  also 
parallels  that  observed  in  Xe  matrices  but  with  an  increased  blue  shift.  The  MCD  pattern 
for  the  long  wavelength  band  group  of  this  region  is  again  not  simulated  well  by  the 
published  gas  phase  data. 

The  270-220  nm  region  for  Co  in  a  Kr  matrix  shows  considerable  evidence  for 
multiple  trapping  sites.  The  absorption  spectrum  (a),  again  with  the  simulated  spectrum 
(b),  are  shown  in  Fig.  6.11.  The  major  transitions  are  readily  assigned  by  comparison  of 
their  shifts  and  intensity  with  the  simulation.  Small  bands  on  the  low  energy  side  of  each 
major  transition  result  from  a  second  site.  The  major  site  appears  blue  shifted  from  the 
gas  phase  while  the  secondary  site  appears  to  the  red  of  the  gas  phase  transitions.  The 
corresponding  MCD  spectra  are  shown  in  Fig.  6.12.  In  the  observed  MCD  spectrum  the 
existence  of  a  third  site  is  apparent  from  the  two  shoulder  bands  on  the  low  energy  side 
of  the  band  at  250.9  nm.  The  overlap  from  multiple  sites  obscures  some  information  but 


Figure  6.9        (a)  Absorption  spectrum  of  Co  atoms  in  a  Kr  matrix,  380-280  nm. 

(b)  Simulated  absorption  spectrum  produced  from  gas  phase  band  positions 
and  intensities. 
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Figure  6.10      (a)  MCD  spectrum  of  Co  atoms  in  a  Kr  matrix,  380-280  nm. 

(b)  Simulated  MCD  spectrum  produced  from  gas  phase  band  positions  and 
intensities. 


Figure  6.11      (a)  Absorption  spectrum  of  Co  atoms  in  a  Kr  matrix,  270-220  nm. 

(b)  Simulated  absorption  spectrum  produced  from  gas  phase  band  positions 
and  intensities. 
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Figure  6.12      (a)  MCD  spectrum  of  Co  atoms  in  a  Kr  matrix,  270-220  nm. 

(b)  Simulated  MCD  spectrum  produced  from  gas  phase  band  positions  and 
intensities. 
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the  major  features  of  the  primary  site  are  readily  identified.  The  assignments  of  several 

bands,  based  on  the  combined  absorption  and  MCD  spectra,    are  listed  in  Table  6.6. 

Again  corrected  assignments  from  the  discussion  are  given  in  parenthesis. 

Cobalt  Atoms  in  Argon  Matrices 

Multiple  sites  are  readily  discemable  even  in  the  absorption  spectra  of  the  380- 
280  nm  region  of  the  Ar  matrices  containing  Co  atoms.  Figure  6.13  shows  the 
experimental  spectrum  and  simulation  as  parts  (a)  and  (b),  respectively.  In  the  360-320 
nm  range  only  the  assignment  of  the  352.3  nm  band  to  a  4D7/2  transition  (gas  phase  341.3 
nm)  seems  appropriate.  Several  intense  bands  appear  in  the  320-280  range  and 
assignment  of  these  bands  is  possible.  The  overlapped  bands  from  multiple  sites  prevent 
conclusive  assignment  of  all  bands.  The  MCD  of  this  region,  shown  in  Fig.  6.14  (a), 
provides  little  additional  help  in  assignment  of  bands  due  to  overlap.  However,  it  is 
possible  to  discern  the  same  MCD  pattern  seen  in  Xe  and  Kr  matrices  for  the  lower 
energy  bands.  Two  distinct  occurrences  of  this  pattern  are  observed.  In  this  case  all 
bands,  including  those  from  other  sites,  are  blue  shifted  from  their  gas  phase  counterparts. 

The  270-220  nm  region  is  more  problematic  but  strong  absorptions  from  the 
major  site  allow  some  assignments  to  be  made.  As  the  absorption  spectrum  shows,  Fig. 
6.15  (a),  there  are  three  substantially  populated  sites,  all  blue  shifted  from  the  gas  phase. 
The  appearance  of  three  negative  bands  in  the  MCD  spectrum,  Fig.  6.16  (a),  confirms  this 
assessment  of  the  absorption  spectrum.  Assignments  are  based  on  the  most  intense  site 
and  parallel  those  of  Ozin     and  Mann  and  Broida    .   Attempts  to  reduce  the  number  of 
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Table  6.6      Band  Assignments,  Matrix  Shifts  and  Observed  MCD  Sign  for  Co  Atoms  in  a 

Krypton  Matrix 


Co  Atoms  in  Kr  matrix  (major  site) 

?9/2  ground  state 

excited  state3' 

Xm(nm) 

vm  (cm  2> 

Avc 

configuration3 

obs.  MCD 

4F 

r9/2 

226.6 

44131 

836 

3d7(a5P)4p 

+ 

4G 

u9/2 

235.1 

42535 

266 

3d74s(b3F)4p 

+ 

4F 

^7/2 

236.8 

42230 

312 

II 

+ 

V 

ull/2 

239.2 

41806 

277 

if 

- 

4F 

r9/2 

240.5 

41580 

354 

11 

hn  ^11/2) 

II 

4D7/2 

250.9 

39856 

207 

11 

+ 

u7/2 

3d8(a3F)4p 

4F 

^7/2 

292.2 

34223 

745 

11 

+ 

2G 

u9/2 

767 

it 

4G7/2 

M 

4F 
r9/2 

298.9 

33456 

614 

11 

+ 

4G 

ull/2 

301.5 

33167 

736 

II 

- 

4D 

u7/2 

307.1 

32563 

535 

If 

+ 

4D 

U7/2 

330.6 

30248 

953 

3d74s(a5F)4p 

+ 

Gll/2    (  F7/2) 

336.5 

29718 

866 

H 

- 

F7/2    (  Gll/2) 

933 

11 

-(w) 

4F 

r9/2 

341.7 

29265 

919 

11 

+ 

agas  phase  data  from  Ref.  20. 

corrected  assignments  in  parenthesis,  see  Discussion 

Av  =  v    -  v 

ra        gas 


Figure  6.13      (a)  Absorption  spectrum  of  Co  atoms  in  an  Ar  matrix,  380-280  nm. 

(b)  Simulated  absorption  spectrum  produced  from  gas  phase  band  positions 
and  intensities. 


141 


Figure  6.14      (a)  MCD  spectrum  of  Co  atoms  in  an  Ar  matrix,  380-280  nm. 

(b)  Simulated  MCD  spectrum  produced  from  gas  phase  band  positions  and 
intensities. 
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Figure  6.15      (a)  Absorption  spectrum  of  Co  atoms  in  an  Ar  matrix,  270-220  nm. 

(b)  Simulated  absorption  spectrum  produced  from  gas  phase  band  positions 
and  intensities. 
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Figure  6.16      (a)  MCD  spectrum  of  Co  atoms  in  an  Ar  matrix,  270-220  nm. 

(b)  Simulated  MCD  spectrum  produced  from  gas  phase  band  positions  and 
intensities. 
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sites  by  annealing  were  only  partially  successful.   Maintaining  the  matrix  at  28K  for  one 

hour  resulted  in  only  minor  shifts  in  band  intensity  among  the  sites  and  a  further  hour  at 

34K  brought   about   considerable  matrix   degradation  but  all   three   sites   were  still 

discernable  in  the  MCD.   Table  6.7  lists  the  assignments  possible  from  this  matrix. 

Discussion 

The  good  agreement  between  experimental  and  simulated  absorption  spectra  for 
the  Co/Xe  case  allows  the  assignment  of  some  observed  bands  based  on  gas  phase  data. 
However,  two  regions  are  identified  by  comparison  of  observed  and  simulated  MCD 
spectra  that  show  a  discrepancy  between  predicted  and  observed  band  patterns.  This 
prevents  concise  assignments  of  the  observed  transitions.  A  similar  situation  exists  in  the 
At  and  Kr  matrices  but  the  experimental  spectra  are  complicated  by  multiple  matrix  sites. 
The  following  discussion  will  resolve  the  conflicts  between  predicted  and  observed 
spectra  by  analysis  of  the  MCD  bandshapes.  The  discussion  will  focus  on  the  Xe  matrix 
since  only  one  site  appears  populated  but  the  conclusions  will  extend  to  the  other 
matrices. 

The  360-330  nm  region  will  be  considered  first.  A  weak  positive  band  in  the 
observed  spectrum  at  349.4  nm  is  readily  assigned  to  a  4F9/2  (352.7  nm)  gas  phase 
transition.  Moving  to  higher  energy  and  using  the  published  gas  phase  data"  a  strong 
positive  band  (347.4  nm,  F7/2)  is  predicted.  Immediately  adjacent  to  this  band  a 
moderately  strong  negative  band  (346.6  nm,  4G11/2)  should  be  seen.  Examination  of  the 
observed  spectrum  shows  a  moderately  strong  negative  band  closely  followed  by  a  strong 
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Table  6.7      Band  Assignments,  Matrix  Shifts  and  Observed  MCD  Sign  for  Co  Atoms  in  an 

Argon  Matrix 


Co  Atoms  in  Ar  matrix  (major  site)  4F9^,  ground  state 


excited  statea'b 


K  (nm) 


vm  (cm"1) 


Avc 


configuration3 


obs.  MCD 


r9/2 

°9/2 

4F 

^7/2 

ull/2 

4F 
r9/2 

J7/2    ^11/2) 

4n 

u7/2 

2G, 


'7/2 


7/2 


'9/2 


'7/2 


9/2 


'11/2 


lD 


>D 


'11/2 


7/2 
7/2 
(%2) 


%2    (4Gll/2) 


9/2 


223.0 
228.6 
229.8 
233.0 
235.1 

245.8 


292.4 
295.0 
301.8 
325.3 


44843 
43745 
43516 
42918 
42535 

40683 


34200 
33898 
33135 
30741 


1548 
1476 
1598 
1389 
1309 

1034 


1358 
1467 
1107 
1446 


3d7(a5P)4p 
3d74s(b3F)4p 


3d8(a3F)4p 


3d74s(a5F)4p 


agas  phase  data  from  Ref.  20. 

bcorrected  assignments  in  parenthesis,  see  Discussion 


uAv  =  v     -  v 
"  m         gas 
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positive  band.   This  is  opposite  to  the  predicted  pattern.   Since  only  a  transition  to  a  J  = 

11/2  sublevel  can  produce  a  negative  MCD  band  from  the  4F9/2  ground  state  and  no  other 

allowed  transitions  to  a  J  =  1 1/2  level  lie  near  this  region  the  observed  pattern  can  only 

result  if  the  energy  levels  of  the  4G11/2  and  4F7/2  states  are  exchanged. 

To  test  this  assumption  the  4F7/2  state  was  assigned  to  the  transition  at  346.6  nm 
and  the  Gn/2  state  was  assigned  to  the  347.4  nm  gas  phase  band.  New  values  of 
gf(C()/2)0)  were  calculated  and  another  simulated  spectrum  produced.  Bandwidths  were 
adjusted  to  approximately  match  those  of  the  observed  spectra.  This  new  simulation, 
along  with  the  original  simulation  and  observed  spectrum  are  shown  in  Fig.  6.17. 
Exchanging  the  two  gas  phase  assignments  clearly  produces  a  better  simulation  of  the 
observed  spectrum  although  two  positive  bands  are  shown  on  the  high  energy  side  of  the 
predicted  negative  band.  A  shift  of  the  reassigned  4G11/2  band  by  50  cm-1  to  higher 
energy  causes  the  two  positive  bands  to  become  partially  overlapping.  The  simulation 
produced  with  this  change  is  also  shown  in  Fig.  6.17.  This  simulation  matches  the 
experimental  band  structure  slightly  better.  The  strong  positive  experimental  band  at 
338.2  nm  is  likely  the  result  of  overlap  of  the  4F7/2  and  4D7/2  bands. 

The  reproduction  of  the  observed  MCD  confirms  that  the  gas  phase  band  at  347.4 
nm  should  be  assigned  as  G11/2  while  the  band  at  346.6  nm  should  be  assigned  as  a 
transition  to  a  F7/2  level.  This  reassignment  also  produces  a  better  agreement  with  the 
Lande  interval  rule  within  both  the  z  G  and  z4F  manifolds.  The  Lande  interval  rule 
states68  that  the  splitting  between  two  adjacent  level  within  a  term  is  proportional  to  the 
higher  J  value.  The  previous  and  present  assignments,  energy  levels,  energy  differences 
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Figure  6.17      Experimental  and  simulated  spectra  of  Co  atoms  in  a  Xe  matrix.    The 
simulations  have  been  shifted  to  align  with  the  experimental  spectrum. 

(a)  Simulation  using  gas  phase  data  published  in  Ref.  20. 

(b)  Experimentally  observed  spectrum. 

(c)  Simulation  with   G11/2  and    F7/2  energy  levels  exchanged. 

(d)  Simulation  with  newly  assigned  4G11/2  shifted  50  cm-1  to  higher 
energy. 
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and  proportionality  constants,  X,  are  listed  in  Table  6.8.     Within  both  manifolds  the 

reassignments  have  brought  the  constants  into  closer  agreement  which  supports  the 

conclusion  from  the  MCD  spectra  that  the  two  energy  levels  were  previously  misassigned. 

The  broad  nature  of  the  negative  band  observed  in  the  270-220  nm  region  and  the 
lack  of  a  predicted  small  positive  band  on  its  low  energy  edge  presents  a  slightly  more 
complicated  problem.  Using  the  procedure  outlined  above,  but  in  this  case  switching  the 
energies  assigned  to  the  4G11/2  and  4F9/2  levels,  the  observed  bandshape  can  be  simulated. 
However,  upon  inspection  of  the  sublevel  energy  spacings  produced  by  this  exchange  it 
is  observed  that  the  two  manifolds  that  originally  followed  the  Lande  interval  rule  are 
now  seriously  distorted.  This  observation  prevents  the  switching  of  the  assignments  from 
being  an  acceptable  explanation  for  the  observed  spectra. 

Another  band,  1°  J  =  7/2,  lies  at  243.6  nm  in  the  gas  phase.  The  assignment  to 
J  =  7/2  is  based  on  the  observed  g- factor  of  1 .40.  Of  the  allowed  transitions  from  a  4F9« 
ground  state,  the  transition  to  a  4D7/2  level  has  a  predicted  g-factor  of  1.43.  The  inability 
of  the  observers  to  assign  a  state  designation  and  the  listing  of  the  g-factor  to  two 
decimal  places  (most  other  g-factors  for  Co  are  given  to  three  places)  allows  speculation 
that  this  band  may  also  be  misassigned.  Moore  also  has  stated  that  much  uncertainty 
may  exist  in  the  published  g-factors  since  many  factors,  including  determination  of  the 
magnetic  field  used  in  Zeeman  experiments,  resolving  power  of  the  spectrograph  and 
experimental  interpretation  of  the  observed  effects,  may  influence  the  determination  of 
the  g  value. 
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150 
The  broadness  of  the  observed  negative  band  indicates  that  another  negative  MCD 

band  may  lie  in  close  proximity  to  the  negative  band  from  the  4G11/2  state.  If  the  1°  band 

is  reassigned  to  J  =  11/2  the  observed  spectra  may  be  simulated  rather  well.    This 

simulation,  the  original  simulation  and  the  observed  spectrum  are  shown  in  Fig.  6.18.   In 

light  of  the  possible  inaccuracies  of  the  observed  g-factor  an  assignment  to  a  J  =  11/2 

level  of  some  unidentified   G  state  in  not  unwarranted.   The  expected  g-value  for  such 

a  state  is  1.27. 

With  the  two  modifications  to  the  previous  band  assignments     the  observed  MCD 

is  simulated  quite  well.  These  observations  illustrate  that  MCD  is  much  more  than  a  tool 

for  confirmation  of  assignments  but  a  technique  that  provides  a  rigorous  test  of  state 

assignments  based  on  absorption  as  well  as  Zeeman  spectra. 
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(a) 


(b) 
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Figure  6.18      Experimental  and  simulated  spectra  of  Co  atoms  in  a  Xe  matrix.    The 
simulations  have  been  shifted  to  align  with  the  experimental  spectrum. 

(a)  Simulation  using  gas  phase  data  published  in  Ref.  20. 

(b)  Experimentally  observed  spectrum. 

(c)  Simulation  with  the  unassigned  J  =  7/2  sublevel  changed  to  a  J  =  11/2 
sublevel. 


CHAPTER  VII 
FUTURE  STUDIES 


The  observation  that  depositions  of  laser  vaporized  Ti  atoms  condensed  with  a 
pulsed  flow  of  At  gas  produced  clean,  single  site  matrices  shows  great  promise  for  the 
preparation  of  matrices  using  this  setup.  Early  attempts  with  this  technique  produced 
numerous  sites  when  a  continuous  flow  of  gas  was  employed  making  spectral 
interpretation  difficult  or  even  impossible  in  some  cases.  The  atomic  studies  of  the 
remaining  first  row  transition  metals,  V  and  Sc,  should  be  completed  using  this  source. 
This  will  give  a  solid  base  of  MCD  spectral  data  for  the  atoms  which  may  then  be  used 
in  the  identification  of  spectral  features  due  to  diatomic  (or  larger)  species. 

The  pulsed  valve  system  was  originally  constructed  to  produce  diatomics  and  is 
patterned  after  the  design  of  sources  used  to  produce  clusters  for  mass  spectrometric 
studies.69  The  source  was  built  within  a  small  brass  block  similar  to  the  one  used  for 
continuous  flow  depositions  but  with  some  modifications.  A  diagram  of  the  block  is 
shown  in  Fig.  7.1.  A  pulse  valve  (General  Valve)  is  mounted  in  a  cutout  in  the  back  of 
the  block  and  the  connection  is  vacuum  sealed  with  an  O-ring  set  in  the  flange  of  the 
valve.  The  gas  pulse  travels  through  a  1  mm  I.D.  x  1  cm  channel  from  where  it  expands 
into  another  channel  (3  mm  I.D.  x  1  cm).  Careful  timing  of  the  laser  pulse  allows  the 
vaporized  metal  to  mix  with  the  expanding  gas  pulse  and  the  high  density  allows 
collisions  which  initiate  cluster  formation.      With  sufficient   backing  pressures   the 
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gas/metal  mixture  should  undergo  supersonic  expansion  on  exiting  the  second  channel. 

This  block  attaches  to  the  Displex  in  the  same  manner  as  shown  in  Fig.  3.2. 

To  insure  a  short  intense  gas  pulse  the  valve  is  backed  with  1-10  atm.  of  matrix 
gas  and  is  over-driven  by  a  400  uS,  200  V  pulse.  The  laser  trigger  and  valve  pulse 
trigger  are  variable  by  ±500  (xS  with  respect  to  their  center  positions  using  the  present 
controller.  This  allows  the  laser  to  fire  a  maximum  of  1  mS  before  (or  after)  the  pulse 
valve  trigger.  These  adjustments  are  critically  important  because  the  pulse  valve  does  not 
open  concurrently  with  its  trigger  pulse  due  to  an  inductive  charge  time  for  the 
electromagnetic  valve  actuator.  The  trigger  pulses  from  the  controller  are  initiated  from 
the  spectrometer  control  computer. 

In  the  preliminary  studies  pressures  were  kept  below  3  atm.  and  repetition  rates 
below  5  Hz.  Using  these  conditions  and  various  delay  settings  a  wide  range  of  results 
were  obtained.  Matrices  produced  ranged  from  those  that  showed  no  metal  atom  spectra 
to  those  that  contained  so  much  metal  that  all  light  from  the  monochromator  was 
absorbed.  Most  commonly  the  matrices  were  found  to  contain  only  metal  atoms  with  no 
evidence  for  dimers  or  other  cluster  species. 

For  this  technique  to  serve  as  a  reliable  source  of  diatomic  molecules  a  wide  range 
of  operating  parameters  must  be  characterized,  including  backing  pressure,  laser  power, 
clustering  chamber  length  and  timing  interactions  between  laser  firing  and  pulse  valve 
opening.  The  most  efficient  way  to  characterize  the  source  would  be  to  attach  it  to  the 
inlet  of  a  mass  spectrometer  and  determine  the  types  of  species  produced  as  these 
parameters  are  varied.   The  inclusion  of  a  small  quadrupole  mass  spectrometer  into  the 
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deposition  system,  located  directly  opposite  the  source  would  allow  optimization  of  the 

source  for  the  species  of  interest  before  deposition.    The  Displex  could  be  mounted  on 

hydraulic  lifters  to  allow  it  to  be  raised  out  of  the  path  for  source  optimization  and 

lowered  into  the  source  output  stream  for  deposition.  Differential  pumping  capability  may 

also  become  necessary  to  prevent  overload  of  the  mass  spectrometer  by  the  large  input 

of  gas  from  the  pulsed  valve.   This  modification  is  stongly  suggested. 

The  production  of  diatomics,  both  hetero-  and  homo-nuclear,  with  this  source 

should  be  possible  and  trapping  of  sufficient  quantities  should  allow  detailed  absorption 

and  MCD  studies  of  their  electronic  spectra.  The  ability  of  MCD  to  extract  excited  state 

information  and  the  dependence  of  observed  band  type  on  the  ground  state  may  allow  the 

identification  of  electronic  ground  states  for  those  species  that  have  not  previously  be 

characterized  and  should  provide  needed  information  for  the  assessment  of  theoretical 

models  for  these  species. 
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Figure  7.1        Cutaway  top  view  of  pulsed  valve,  laser  vaporization  cluster  source. 


APPENDIX 
COMPUTER  PROGRAMS 


This  appendix  contains  listings  of  several  computer  programs  used  throughout  the 
course  of  this  work.  In  several  instances  short  programs  were  written  to  plot  the  digitized 
data  to  the  screen  or  chart  recorder  and  to  convert  the  data  structure  used  for  storing 
spectral  data  to  disk  to  a  format  usable  by  other  programs.  The  spectrometer  control 
program  and  the  moments  program  will  be  listed.  The  programs  were  written  and 
compiled  for  stand-alone  operation  using  Microsoft  QuickBasic  4.5.  Mr.  Marc  Eyring 
and  the  author  have  written  and  modified  the  programs  as  needed  for  proper  functioning 
and  the  versions  used  during  this  work  are  given. 


SP.BAS    —   spectrometer  control  program 

dim 

a%(6001),b%(600 1  ),c%(500),dio%(4),fod$(  1 5),f od  1  $(1 5,20),fod2$(200),f  $(  1 5),rec(  1 00),desc$(  1 5),dat$(  1 5) 

dio%(0)=0:dio%(l)=0:dio%(2)=0:dio%(3)=0 

on  error  goto  errorhandler 

'DASH16  Initialization 

md%=0 

dio%(0)=&H300 

dio%(l)=2 

dio%(2)=3 

flag%=0 

call  DASH16(md%,dio%(0),flag%) 

if  flag%<>0  then  print  "Initialization  Error  #  ";Qag%:stop 

'Set  multiplexer  scan  limits 
md%=l 
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dio%(0)=0 

dio%(l)=l 

call  DASH16(md%,dio%(0),flag%) 

if  flag%<>0  then  print  "Error  in  setting  scan  limits.  Error  #  ";flag%:stop 

'Set  programmable  timer  or  use  external  trigger 
md%=17 

setup: 

dio%(0)=0 

gosub  dio 

print  "Input  initial  and  final  wavelengths  over  scan  range  -  " 

input  x,y 

wl=x 

gosub  pem 

print  "  ":print  "Number  of  cycles  for  signal  averaging  -  "; 

input  e:e=2*e 

print  "  "rprint  "Chart  recorder  spread:  5  A/mm(3),2.5  A/mm(6),1.25  Aymm(12),0.625  A/mm(24)" 

print  "Default  is  5  A/mm(CR)" 

90  a$=inkey$ 

if  a$ochr$(13)  and  a$<>"3"  and  a$<>"6"  and  a$<>"12"  and  a$<>"24"  goto  90 

if  a$=chr$(13)  then 

d=3 

else  d=val(a$) 
end  if 
print  "  ":print  "  ":print  "Data  will  be  collected  every  1  Angstrom" 
print  "  ":print  "Ready  to  begin  scan  (y)" 
print  "Reinitialize  scan  parameters  (n)  -  "; 
100  a$=inkey$ 

if  a$o"y"  and  a$o"Y"  and  a$o"n"  and  a$o"N"  then  goto  100 
print  a$ 

if  a$="n"  or  a$="N"  goto  setup 
if  a$="y"  or  a$="Y"  then  gosub  scan 
210  i$=inkey$ 

if  i$o"x"  and  i$o"X"  goto  210 
screen  2 
els 

print  "Do  you  wish  to  save  the  data  ?  "; 
110a$=inkey$ 

if  a$o"y"  and  a$o"Y"  and  a$o"n"  and  a$o"N"  goto  110 
print  a$ 

if  a$="Y"  or  a$="y"  then  gosub  datasave 

print  "  ":print  "Do  you  wish  to  slew  the  monochromator  back  to" 
print  "its  starting  point  ?  -  "; 
120  a$=inkey$ 

if  a$o"y"  and  a$o"Y"  and  a$o"n"  and  a$o"N"  then  goto  120 
print  a$ 

if  a$="y"  or  a$="Y"  then  gosub  slew 
print:print:print  "Type  'X'  to  return  to  DOS" 
print 
125  a$=inkey$ 
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if  a$o"X"  and  a$o"x"  goto  125 
end 

scan: 
screen  1 
color  0,1 
els 

n4=x-y:span=x-y 
line  (10,10)-(310,180)„b 
line  (10,95)-(310,95) 
locate  7,l:print "+" 
locate  18,l:print  "-" 
locate  l,16:print  "ABS  &  MCD" 
res%=int(span/300) 

if  res%<span/300  then  res%=int(span/300)+l 
spcg=100/res% 
for  i=0  to  int(span/100) 
line  (10+spcg*i,180)-(10+spcg*i,175) 
next 

locate  24,l:print  x: 
locate  24,35:print  y; 

for  il=0  to  n4 

if  (il/50)=int(il/50)  then  gosub  event 

dio%(0)=e 

dio%(2)=0 

dio%(3)=0 

md%=4 

dio%(l)=VARPTR(c%(0)) 

call  DASH16(md%,dio%(0),flag%) 

if  flag%<>0  then  print  "Error  during  Mode  4  A/D  operation.  Error  #  ";flag%:stop 

if  il=n4  goto  130 

for  i=l  to  4 

dio%(0)=4 

gosub  dio 

dio%(0)=5 

gosub  dio 

gosub  chart 

next  i 

130wl=wl-l 

gosub  pem 

gosub  datared 

if  il/res%=int(il/res%)  then 

ampmcd=(a%(il)/4095)*  170 

ampabs=(b%(il)/4095)*170 

pset  (10+il/res%,180-int(ampmcd)),l 

pset  (10+il/res%,180-int(ampabs)),2 
end  if 
next  il 

locate  25,l:print  "Enter  'X'  to  continue"; 
return 
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dio: 

md%=13 

call  DASH16(md%,dio%(0),fiag%) 

if  flag%<>0  then  print  "Error  in  Mode  13   digital  I/O.  Error  #  ";flag%:stop 

for  k=l  to  2 

next 

return 

chart: 

for  j=l  to  d 

dio%(0)=l 

gosub  dio 
dio%(0)=5 
gosub  dio 
next  j 
return 

event: 

dio%(0)=13 
gosub  dio 
for  jj=l  to  300 
next  jj 
dio%(0)=5 
gosub  dio 
return 

pern: 

md%=15 

v=2*(0.105e-3*(wl/2)*4095)+0.5 

dio%(l)=int(v) 

dio%(0)=0 

call  DASH16(md%,dio%(0)41ag%) 

if  flag%<>0  then  print  "Error  Mode  15  D/A  output.  Error  #  ";flag%:stop 

return 

datared: 

aa#=0:bb#=0 

for  k=0  to  e-2  step  2 

aa#=aa#+c%(k) 

bb#=bb#+c%(k+l) 

next  k 

aa#=aa#/(e/2) 

a%(il)=int(aa#) 

bb#=bb#/(e/2) 

b%(il)=int(bb#) 

return 

slew: 

for  i=l  to  (4*n4) 

dio%(0)=2 

gosub  dio 
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dio%(0)=3 

gosub  dio 

forjk=l  to  50 

next 

next 

return 

filedesc: 

data    "DIRECTORY    PATHNAME:","SAMPLE:","MATRIX:","DEPOS.    TEMP.","MATRIX    TEMP. 

(THERMOCOTJPLE):","MAGNETIC  FIELD  STRENGTH","FEEDBACK  VOLTAGE:" 

data    "SENSITIVITY:", "TIME   CONSTANT:","*   OF    SIGNAL   AVERAGED    CYCLES: "."INITIAL 

WAVELENGTH:" 

data  "FINAL  WAVELENGTH:","RESOLUTION:")"DATE:" 

flagl=0:lino=0 

print:print:print  "Type  'initiate'  to  create  new  configuration  file.  Enter  'mod'  to  access  file  for  modification 

410  input  a$ 

if  a$o"MOD"  and  a$o"mod"  and  a$o"INITIATE"  and  a$o"initiate"  goto  410 
if  a$="MOD"  or  a$="mod"  then  flagl=2 
print: print: print  "Enter  data  information ":print 
reen: 

for  i=l  to  14 
if  flagl=0  then 
read  f$(i) 
print  f$(i); 
elseif  flagl=l  then 
i=lino 
print  f$(i) 
elseif  flag  1=2  then 
open  "filedesc.dat"  as  #2  len=2 
field#2,  2  as  seg$ 
get#2,l 

numrec=cvi(seg$) 
dat$="":k=-l 
for  1=2  to  numrec 
get#2,l 

fod2$(l)=seg$ 
if  left$(fod2$(I),l)ochr$(34)  then 

dat$=dat$+fod2$(l) 
else  k=k+l 
fod$(k)=dat$ 
dat$="" 
end  if 
next 
i=14 

goto  155 
end  if 
150  input  fod$(i) 

if  fod$(i)=""  goto  150 
if  i=l  then  mkdir  fod$(i) 
if  flagl=l  then  i=14 
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155  next 
restore 

cls:print:print 
for  i=l  to  14 

read  f$(i) 

print  using  "##########" ;i; 

print  " "; 

print  using  "\  \";f$(i); 

print  fod$(i) 
next 
restore 
print:print:print  "Enter  'y'  if  data  is  correct,  'n'  if  you  wish  to  re-enter  -  " 

160  a$=inkey$ 

if  a$o"Y"  and  a$o"y"  and  a$o"N"  and  a$o"n"  goto  160 
if  a$="N"  or  a$="n"  then 
flagl  =  l 

161  printrinput  "Number  of  data  line  to  be  re-entered  -  ";lino 
if  liiux  1  or  lino>14  then 

print  "Data  line  out  of  range  1-14  " 
goto  161 
end  if 
goto  reen 
else 

open  "filedesc.dat"  as  #1  len=2 
field#l,  2  as  seg$ 
gosub  stringstruc 
close#l 
return 
end  if 


datasave: 

print 

print  "ABS(A),  MCD(M),  or  BOTH(B)  ?  "; 

140  d$=inkey$ 

if  d$o"A"  and  d$o"a"  and  d$o"M"  and  d$o"m"  and  d$o"B"  and  d$o"b"  then  goto  140 

print  d$ 

gosub  filedesc 

fod$(l)=left$(fod$(l),8) 

if  d$="A"  or  d$="a"  then  gosub  absorp 

if  d$="M"  or  d$="m"  then  gosub  mcd 

if  d$="B"  or  d$="b"  then  gosub  both 

return 


absorp: 

fod$(0)="ABSORPTION  FILE" 

cls:print:print:print  "ABS  FILE":print:print 

print:print:input  "Filename  (ABS)  -  ";absfile$ 

chdir  fod$(l) 

open  absfileS  as  #1  len=2 

field#l,  2  as  seg$ 
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gosub  stringstruc 
for  i=0  to  n4 

wJ=x-i 

lset  seg$=mki$(bf7c(i)) 

put#l,wl 
next 
close#l 

chdir  "\qbasic\spec" 
return 

mcd: 

cls:print:print:print  "MCD/MLD  FILE":print:print:print:print 

135  input  "Type  of  file  (MCD  or  MLD)  -  ";mdl$ 

if  mdl$o"MCD"  and  radl$o"mcd"  and  mdl$o"MLD"  and  mdl$o"rnld"  goto  135 

if  (radl$="MCD")  or  (mdl$="mcd")  then 

fod$(0)="MCD  FILE" 

mdl$="MCD" 

elseif  (mdl$="MLD")  or  (mdl$="mld")  then 

fod$(0)="MLD  FILE" 

mdl$="MLD" 

end  if 
print:print:print  "Filename  (";mdl$;")  -  "; 
input  mcdfileS 
chdir  fod$(l) 

open  mcdfile$  as  #1  len=2 
field#l,  2  as  seg$ 
gosub  stringstruc 
for  i=0  to  n4 
wl=x-i 

lset  seg$=mki$(a<7c(i)) 
put  #l,wl 
next 
close#l 

chdir  "\qbasic\spec" 
return 

both: 

gosub  absorp 

gosub  mcd 

cls:print:print:print  "ABS  Filename  -  ";fod$(l);"\";absfile$ 

print:print  mdl$;"  Filename  -  ";fod$(l);"\";mcdfile$ 

return 

stringstruc: 
numrec=2 
for  i=0  to  14 

rec(i)=len(fod$(i))/2 

if  int(rec(i))=rec(i)  then 
rec(i)=int(rec(i)) 

else  rec(i)=int(rec(i))+l 

end  if 


163 


rel=l 

for  j=l  to  rec(i) 
fodl$(i,]>mid$(fod$(i),rel,2) 
rel=rel+2 

lset  seg$=fodl$(i,j) 
put#l,numrec 
numrec=numxec+ 1 
next 

b$=chi$(34) 
lset  seg$=b$ 
put#l,numrec 
numrec=numrec+ 1 
next 

numrec=numrec-l 
lset  seg$=mki$(numrec) 
put#l,l 
return 

errorhandler: 

'error  75  means  directory  fod$(l)  already  exists 

if  err=75  then 

cls:print:print 

print  "Directory  already  exists  -  type  <CR>  to  continue" 

170  a$=inkey$ 

if  inkey$ochr$(13)  goto  170 

els 

resume  next 

end  if 

on  error  goto  0 


ALLMOMS.BAS  —  program  to  calculate  absorption  and  MCD  moments  from  digitized 

data 


DIM  b%(6001),  fod$(14),  fodl$(500),  file$(100),  pp(100) 
num  =  0 


This  program  has  been  modified  to  store  Oth.lst  and  2nd  moments  for  ABS 
'and  0th,  1st,  2nd  and  3rd  moments  for  MCD  to  be  imported  into  a 
'spreadsheet  for  data  manipulation.    DCR  4/91 


PRINT  :  PRINT  :  PRINT  "This  program  calculates  moments  for  the  specified  file.":  PRINT  :  PRINT 
GOTO  decision 

prog: 
CLS 
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IF  num  <>  0  THEN 

absfileS  =  path$  +  file$(jk) 

ELSE  INPUT  "Filename  -  ";  absfileS 
END  IF 

OPEN  absfileS  FOR  RANDOM  AS  #1  LEN  =  2 
FIELD  #1,  2  AS  seg$ 
GET  #1,  1 
numrec  =  CVI(seg$) 
dat$  =  "":  k  =  -1 
FOR  i  =  2  TO  numrec 
GET#1,  i 
fodl$(i)  =  seg$ 
IF  LEFT$(fodl$(i),  1)  <>  CHR$(34)  THEN 

dat$  =  dat$  +  fodl$(i) 
ELSE 

k  =  k+  1 

fod$(k)  =  dat$ 

PRINT  k;  fod$(k) 

dat$  =  "" 
END  IF 
NEXT 

x  =  VAL(fod$(ll)) 
y  =  VAL(fod$(12)) 
type$  =  LEFT$(fod$(0),  3) 
PRINT  "x=";  x,  "y=";  y,  "type=";  type$ 
IF  num  <>  0  THEN  RETURN 

decision: 

num  =  0 

PRINT  :  PRINT  :  PRINT  "Type  'm'  to  calculate  moments,  'q'  to  quit  -  "; 

10  a$  =  INKEYS 

IF  a$  <>  "m"  AND  a$  <>  "M"  AND  a$  <>  "Q"  AND  a$  <>  "q"  GOTO  10 

IF  a$  =  "M"  OR  a$  =  "m"  THEN 

GOTO  moments 
ELSEIF  a$  =  "Q"  OR  a$  =  "q"  THEN 

GOTO  finish 
END  IF 

moments: 

CLOSE  #1 

CLS 

PRINT  :  PRINT  "Do  you  wish  to  save  the  moments  to  disk?" 

PRINT  :  PRINT  "The  program  will  automatically  save  0th,  1st  and  2nd  moments  for  absorption  or" 

PRINT  "0th,lst,2nd  and  3rd  moments  for  MCD  if  you  answer  Yes(Y)." 

INPUT  f$ 

IF  f$  =  "Y"  OR  f$  =  "y"  THEN 

PRINT  :  PRINT  :  PRINT  "Enter  filename  for  saving  moments  (this  must  be  different  from  data  filename)" 
PRINT  :  PRINT  "Also  use  a  different  name  for  ABS  and  MCD  moments  calculations." 
PRINT  :  PRINT  "NOTE:  Filenames  must  end  with  '.DAT  to  be  used  with  Quattro." 
INPUT  name$ 
END  IF 
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CLS 

PRINT  :  PRINT  :  INPUT  "Enter  number  of  files  on  which  calculations  are  to  be  performed  -  ";  num 

PRINT  :  PRINT  :  INPUT  "Enter  type  of  files  to  be  processed:  ABS  or  MCD";  typ$ 

IF  typ$  =  "MCD"  OR  typ$  =  "mcd"  THEN 

PRINT  :  PRINT  "Do  you  wish  to  enter  the  average  energy  from  absorption  moments  calculations  for  each 

file  "; 

INPUT  d$ 

END  IF 

PRINT  :  PRINT  :  INPUT  "Enter  common  pathname  for  files  -  ";  path$ 

IF  d$  =  "Y"  OR  d$  =  "y"  THEN 

PRINT  :  PRINT  "Enter  'filename.average  energy'  at  prompt  ":  PRINT 

ELSE 

PRINT  :  PRINT  "Enter  filename(s)  at  prompt  ":  PRINT 
END  IF 

FOR  ik  =  1  TO  num 
IF  d$  =  "y"  OR  d$  =  "Y"  THEN 

INPUT  file$(ik),  pp(ik) 

ELSE 

INPUT  file$(ik) 
END  IF 
NEXT 
CLS 

PRINT  :  PRINT  :  PRINT  "Enter  number  of  bands  on  which  calculations  are  to  be  performed  -  "; 
INPUT  nb 

PRINT  :  PRINT  :  PRINT  "Input  matrix  temp,  if  different  from  indicated  temp,  (default  <CR>  is  indicated 
temp.)  -  "; 
INPUT  b$ 

IF  b$  <>  ""  THEN  fod$(5)  =  b$ 
CLS 
k=  1 

PRINT  :  PRINT  :  PRINT  "Data  for  band  #";  k 

PRINT  :  PRINT  "Do  you  want  the  moments  output  to  the  printer  ";  :  INPUT  z$ 

IF  z$  =  "y"  OR  z$  =  "Y"  THEN 

PRINT  :  PRINT  :  PRINT  "For  standard  or  NLQ  printer  mode  enter  <CR>.  In  compact  mode  enter  <CR> 
or  T" 
PRINT  "for  printing  on  left  half  of  paper,  enter  V  for  printing  on  right  half  -"; 
INPUT  s$ 

IF  s$  =  ""  THEN  s$  =  "1" 
CLS 

END  IF 

PRINT  :  PRINT  "For  wavelength  limits  of  band,  ensure  that  there  are  at  least  10  data  points  to  either  side 
of  the  limits  for  baseline  calculation." 

PRINT  :  PRINT  "Enter  long  wavelength  limit  of  band  -  ";  :  INPUT  xl 

PRINT  :  PRINT  "Enter  short  wavelength  limit  of  band  -  ";  :  INPUT  xs 

PRINT  :  PRINT  "Enter  calibration  factor  -  ";  :  INPUT  cbrt 

IF  d$  <>  "Y"  AND  d$  <>  "y"  THEN 
PRINT  :  PRINT  "Enter  approximate  band  center  -  ";  :  INPUT  p 

END  IF 

FOR  jk  =  1  TO  num 

GOSUB  prog 

sp  =  xl  -  xs 
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FOR  i  =  xl  +  5  TO  xs  -  5  STEP  -1 
GET#1,  i 

b%(xl  +  5  -  i)  =  CVI(seg$) 
NEXT 
zl  =  0 

FOR  j  =  0  TO  10 
zl  =  zl  +  b%(j) 
NEXT 
zl  =  zl/  11 
zs  =  0 

FOR  j  =  sp  TO  sp  +  10 
zs  =  zs  +  b%(j) 
NEXT 
zs  =  zs  /  11 
flag  =  0 

IF  d$  =  "Y"  OR  d$  =  "y"  THEN 
cr  =  pp(jk):  p2  =  pp(jk) 
ELSE 

cr  =  p:  p2  =  p 
END  IF 

The  basics  of  this  routine  came  from  program  'moments'  written  before  I  came 
'here  for  use  on  the  old  MCD  apparatus,  by  M.  Kregkwski. 
'Although  on  first  inspection  it  may  appear  that  the  routine  is  not  correct, 
'close  inspection  shows  that  it  is.  Variables  and  certain  proportional 
'factors  are  present  but  are  in  locations  and/or  forms  different  than  one 
'might  expect.  MWE  24  Feb  88 

1000  p2  =  2  /  p2 
a  =  (zs  -  zl)  /  (xl  -  xs) 
xl  =  1  /  xl 
z3  =  b%(5)  -  zl 
zl  =  z3  *  xl 

s  =  0:  area  =  0:  t  =  0:  u  =  0:  v  =  0 
j=6 
FOR  i  =  xl  -  1  TO  xs  STEP  -1 

x2=  1  /i 

z4  =  b%(j)  -  a  *  (xl  -  i)  -  zl 

z2  =  z4  *  i 

r  =  (x2  -  xl)  *  (zl  +  z2) 

q  =  (x2  -  xl)  *  (z3  +  z4) 

area  =  area  +  q 

s  =  s  +  r 

t  =  t  +  r  *  (x2  +  xl  -  p2) 

u  =  u  +  r  *  (x2  +  xl  -  p2)  "  2 

v  =  v  +  r  *  (x2  +  xl  -  p2)  "  3 

xl  =  x2 

zl  =z2 

z3  =  z4 

j=j  +  l 

NEXT 

ma  =  .5  *  area  *  cbrt 
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mO  =  .5  *  s  *  cbrt 

IF  type$  =  "mcd"  OR  type$  =  "MCD"  GOTO  1010 
PRINT  "flag  =  ";  flag 
pe  =  2  *  ma  /  mO 

IF  ABS(pe  -  p2)  <  IE- 10  GOTO  1010 
p2  =  mO  /  ma 
flag  =  flag  +  1 
IF  flag  =  5  THEN 

PRINT  :  PRINT  :  PRINT  :  PRINT  :  PRINT  :  PRINT  :  PRINT  "Unable  to  process  file  -  ";  absfileS 
PRINT  :  PRINT  :  PRINT  :  PRINT  :  PRINT  :  PRINT  :  PRINT 
IF  z$  =  "Y"  OR  z$  =  "y"  THEN 

WIDTH  LPRINT  160 

LPRINT  :  LPRINT  :  LPRINT  :  LPRINT  :  LPRINT  :  LPRINT 

IF  s$  =  "L"  OR  s$  =  "1"  THEN 

LPRINT  "Unable  to  process  file  -  ";  absfileS 

ELSEIF  s$  =  "R"  OR  s$  =  "r"  THEN 

LPRINT  SPC(65);  "Unable  to  process  file  -  ";  absfileS 

END  IF 

LPRINT  :  LPRINT  :  LPRINT  :  LPRINT  :  LPRINT  :  LPRINT  :  LPRINT 
END  IF 
GOTO  1100 
ELSE 

GOTO  1000 
END  IF 

1010  ml  =  2.5E+07  *  t  *  cbrt 
m2  =  1.25E+15  *  u  *  cbrt 
m3  =  6.25E+22  *  v  *  cbrt 
IF  f$  =  "Y"  OR  f$  =  "y"  THEN 
OPEN  nameS  FOR  APPEND  AS  #2 
IF  typ$  =  "ABS"  OR  typ$  =  "abs"  THEN 

WRITE  #2,  mO,  ml,  m2 
ELSEIF  typ$  =  "MCD"  OR  typ$  =  "mcd"  THEN 

WRITE  #2,  mO,  ml,  m2,  m3 
END  IF 
CLOSE  #2 
END  IF 
CLS 

PRINT  :  PRINT  "  File  ";  absfileS,  typeS 

PRINT  :  PRINT  "Matrix  temperature:  ";  fod$(5) 
PRINT  "Area  =  ";  ma  *  1E+08 
PRINT  "0th  moment  =  ";  mO,  "long  wave  limit  =  ";  xl 
PRINT  "1st  moment  =  ";  ml,  "short  wave  limit  =  ";  xs 
PRINT  "2nd  moment  =  ";  m2,  "center  of  band  =  ";  cr 
PRINT  "3rd  moment  =  ";  m3,  "zl  =  ";  zl;  "zs  =  ";  zs 
IF  typeS  =  "ABS"  THEN 

PRINT  :  PRINT  "average  energy  at  ";  mO  /  ma;  "  Angstroms" 
ELSE  PRINT 
END  IF 

PRINT  "calibration  factor  =  ";  cbrt 
IF  z$  =  "y"  OR  z$  =  "Y"  THEN 

WIDTH  LPRINT  160 
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IF  s$  =  "1"  OR  s$  =  "L"  THEN  GOSUB  prtrl 

IF  s$  =  "r"  OR  s$  =  "R"  THEN  GOSUB  prtr2 
END  IF 

1100  CLOSE  #1 
NEXT 
GOTO  decision 

prtrl: 

LPRINT  "       File  ";  absfile$,  type$ 

LPRINT  :  LPRINT  "Matrix  temperature:  ";  fod$(5) 

LPRINT  "Area  =  ";  ma  *  1E+08 

LPRINT  "0th  moment  =  ";  mO,  "long  wave  limit  =  ";  xl 

LPRINT  "1st  moment  =  ";  ml,  "short  wave  limit  =  ";  xs 

LPRINT  "2nd  moment  =  ";  m2,  "center  of  band  =  ";  cr 

LPRINT  "3rd  moment  =  ";  m3,  "zl  =  ";  zl;  "zs  =  ";  zs 

IF  type$  =  "ABS"  THEN 

LPRINT  :  LPRINT  "Average  energy  at  ";  mO  /  ma:  "  Angstroms" 

ELSE  LPRINT  :  LPRINT 

END  IF 

LPRINT  "calibration  factor  =  ";  cbrt 

LPRINT  :  LPRINT  :  LPRINT 
RETURN 

prtr2: 

LPRINT  SPQ65);  "       File  ";  absfile$,  type$ 
LPRINT  :  LPRINT  SPC(65);  "Matrix  temperature:  ";  fod$(5) 


"Area  =  ";  ma  *  1E+08 

"0th  moment  =  ";  mO,  "long  wave  limit  =  ";  xl 
"1st  moment  =  ";  ml,  "short  wave  limit  =  ";  xs 
"2nd  moment  =  ";  m2,  "center  of  band  =  ";  cr 
"3rd  moment  =  ";  m3,  "zl  =  ";  zl;  "zs  =  ";  zs 


LPRINT  SPQ65); 

LPRINT  SPQ65); 

LPRINT  SPQ65); 

LPRINT  SPQ65); 

LPRINT  SPQ65); 

IF  type$  =  "ABS"  THEN 

LPRINT  :  LPRINT  SPC(65);  "Average  energy  at  ";  mO  /  ma;  "  Angstroms" 

ELSE  LPRINT  :  LPRINT 

END  IF 

LPRINT  SPQ65);  "calibration  factor  =  ";  cbrt 

LPRINT  :  LPRINT  :  LPRINT 

RETURN 


finish: 
CLOSE  #1 
END 
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Another  program  written  by  the  author  but  not  used  for  the  work  reported  here 

was  designed  to  control  the  ramp  of  the  lab's  superconducting  magnet  and  take  data  at 

various  fields  during  the  ramp  for  saturation  studies. 

SATURATE. BAS  —  Superconducting  magnet  ramp  control 


DECLARE  SUB  datasave  (t%) 

DECLARE  SUB  chart  0 

DECLARE  SUB  persist  0 

DECLARE  SUB  nopersist  0 

DECLARE  SUB  instruct!  () 

DECLARE  SUB  fastramp  (hifield%,  fast$) 

DECLARE  SUB  fieldramp  (hifield%,  dataspace%,  directions) 

DECLARE  SUB  datacollect  (t%) 

DECLARE  FUNCTION  pem%  (wl%) 

DECLARE  SUB  DASH16  (md%,  BYVAL  dummy%,  flag%) 

DIM  dio%(4),  c%(500),  mcd%(500),  magfield%(500) 

COMMON  SHARED  dio%0,  c%0,  t  AS  INTEGER 


'this  program  is  written  in  QuickBasic  4.5. 

This  program  is  used  to  control  the  superconducting  magnet  in  constant 
'current  mode.    It  will  allow  ramping  up  to  and  down  from  30000  Gauss. 
'It  will  also  allow  the  user  to  rapidly  ramp  the  current  down  after  the 
'magnet  is  in  persistent  mode  or  up  to  take  the  magnet  out  of  persistent 
'mode. 

DEFINT  A-Z 


dio%(0)  =  0:  dio%(l)  =  0:  dio%(2)  =  0:  dio%(3)  =  0 
dp  =  0 

'DASH  16  Initialization 

md%  =  0 

dio%(0)  =  &H300 

dio%(l)  =  2 

dio%(2)  =  3 

flag%  =  0 

CALL  DASH16(md%,  VARPTR(dio%(0)),  flag%) 

IF  flag%  <>  0  THEN  PRINT  "Initialization  Error  #  ";  flag%:  STOP 

'Set  multiplexer  scan  limits 
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md%  =  1 

dio%(0)  =  0 

dio%(l)  =  1 

CALL  DASH16(md%,  VARPTR(dio%(0)),  flag%) 

IF  flag%  <>  0  THEN  PRINT  "Error  in  setting  scan  limits.  Error  #  ";  flag%:  STOP 

'Set  both  D/A  channels  to  zero 

md%  =  16 

dio%(l)  =  0 

dio%(0)  =  0 

CALL  DASH16(md%,  VARPTR(dio%(0)),  flag%) 

IF  flag%  <>  0  THEN  PRINT  "Error  Mode  16  D/A  output.  Error  #  ";  flag%:  STOP 

CALL  instruct  1 

CLS 

PRINT  :  PRINT 

INPUT  "  Choose  whether  to  ramp  UP  from  no  field  or  down(DN)  from  high  field    ",  direct! 
5  SELECT  CASE  direct! 
CASE  IS  =  "up" 

direction!  =  direct! 

INPUT  "Enter  the  wavelength  maximum  for  the  band  in  Angstroms.  ",  wl% 
pemvolt  =  pem(wl%) 
PRINT 
10  INPUT  "Enter  maximum  field  value  in  Gauss  (up  to  30000).  ",  hifield% 

IF  hifield%  >  30000  THEN 

PRINT  "  Maximum  allowed  field  is  30000  gauss." 
GOTO  10 
END  IF 
PRINT 

INPUT  "Enter  spacing  for  data  collection  in  Gauss.  ",  dataspace% 
PRINT  :  PRINT  :  PRINT  "Type  'V  when  ready  to  begin" 
DO 

choice!  =  INKEY! 
LOOP  WHILE  choice!  <>  "y" 
CALL  fieldramp(hifield%,  dataspace%,  direction!) 
CLS 

PRINT  :  PRINT  "Do  you  wish  to  save  the  data  to  disk  (y  or  n)?  " 
70  a!  =  INKEY! 

IF  a!  <>  "Y"  AND  a!  <>  "y"  AND  a!  <>  "N"  AND  a!  <>  "n"  GOTO  70 
IF  a!  =  "Y"  OR  a!  =  "y"  THEN 

CALL  datasave(t%) 
END  IF 
CLS 

PRINT  :  PRINT  "Do  you  wish  to  place  magnet  in  persistant  mode?" 
PRINT  :  PRINT  "Type  'Y'  for  persistant  mode  or  N'  to  ramp  the  magnet  down." 
20  a!  =  INKEY! 

IF  a!  <>  "Y"  AND  a!  <>  "y"  AND  a!  <>  "N"  AND  a!  <>  "n"  GOTO  20 
IF  a!  =  "y"  OR  a!  =  "Y"  THEN 
CALL  persist 
CLS 
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PRINT  :  PRINT 
fastS  =  "fastdown" 
CALL  fastramp(hifield%,  fastS) 
END  IF 

IF  a$  =  "n"  OR  a$  =  "N"  THEN 
directions  =  "down" 

PRINT  :  PRINT  "Do  you  wish  to  collect  data  on  the  ramp  down  (Y  or  N)?" 
50  b$  =  INKEY$ 

IF  b$  <>  "Y"  AND  b$  <>  "y"  AND  b$  <>  "N"  AND  b$  <>  "n"  GOTO  50 
IF  b$  =  "Y"  OR  b$  =  "y"  THEN 

PRINT  :  PRINT  "Do  you  wish  to  change  data  spacing?  (Y  or  N)   " 
60  d$  =  INKEYS 

IF  d$  <>  "Y"  AND  d$  <>  "y"  AND  d$  <>  "N"  AND  d$  <>  "n"  GOTO  60 

IF  d$  =  "Y"  OR  d$  =  "y"  THEN 
PRINT  :  INPUT  "Enter  new  data  spacing  in  gauss.    ".  dataspace% 

END  IF 
END  IF 

IF  b$  =  "N"  OR  b$  =  "n"  THEN 
dataspace%  =  0 
END  IF 

PRINT  "Type  'Y'  when  ready  to  ramp  down." 
DO 

choices  =  INKEYS 
LOOP  WHILE  choices  <>  "y" 
CALL  fieldramp(hifield%,  dataspace%,  directions) 
CLS 

PRINT  :  PRINT  "Do  you  wish  to  save  the  data  to  disk  (y  or  n)?  " 
80  a$  =  INKEYS 

IF  a$  <>  "Y"  AND  a$  <>  "y"  AND  a$  <>  "N"  AND  a$  <>  "n"  GOTO  80 
IF  a$  =  "Y"  OR  a$  =  "y"  THEN 

CALL  datasave(t%) 
END  IF 

END  IF 

END 
CASE  IS  =  "dn" 
40  INPUT  "Enter  the  present  field  in  gauss  (max=30000).  ";  hifield% 

IF  hifield%  >  30000  GOTO  40 

fastS  =  "fastup" 

PRINT  :  PRINT  "Type  'y'  when  ready  ramp  up  current  on  magnet  leads." 

DO 

e$  =  INKEYS 

LOOP  WHILE  e$  <>  "y" 

CALL  fastramp(hifield%,  fastS) 

directions  =  "down" 

CALL  nopersist 

CLS 

PRINT  :  PRINT  "Do  you  wish  to  collect  data  on  the  ramp  down  (Y  or  N)?  " 
30  b$  =  INKEYS 
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IF  b$  <>  "Y"  AND  b$  <>  "y"  AND  b$  <>  "N"  AND  b$  <>  "n"  GOTO  30 
IF  b$  =  "N"  OR  b$  =  "n"  THEN 
dataspace^  =  0 
END  IF 

IF  b$  =  "Y"  OR  b$  =  "y"  THEN 

PRINT  :  PRINT 

INPUT  "Enter  the  wavelength  maximum  for  the  band  in  Angstroms.  ",  wl% 

peravolt  =  pem(wl%) 

PRINT 

PRINT 

INPUT  "Enter  spacing  for  data  collection  in  Gauss.  ",  dataspace% 

END  IF 

PRINT  :  PRINT  :  PRINT  "Type  'Y1  when  ready  to  begin" 

DO 

choice$  =  INKEYS 
LOOP  WHILE  choices  <>  "y" 

CALL  fieldramp(hifield%,  dataspace%,  directions) 
CLS 

PRINT  :  PRINT  "Do  you  wish  to  save  the  data  to  disk  (y  or  n)?  " 
90  a$  =  INKEYS 

IF  a$  <>  "Y"  AND  a$  <>  "y"  AND  a$  <>  "N"  AND  a$  <>  "n"  GOTO  90 
IF  a$  =  "Y"  OR  a$  =  "y"  THEN 

CALL  datasave(t%) 
END  IF 
END 
END  SELECT 


SUB  chart 
SHARED  dio%0 
PRINT  "called  chart" 
md%  =  13 
FOR  j  =  1  TO  200 

dio%(0)  =  1 

CALL  DASH16(md%,  VARPTR(dio%(0)),  flag%) 

IF  flag%  <>  0  THEN  PRINT  "Error  in  Mode  13  digital  I/O.   Error*";  flag%:  STOP 

FOR  k  =  1  TO  2 

NEXT 

PRINT  dio%(0) 

dio%(0)  =  5 

CALL  DASH16(md%,  VARPTR(dio%(0)),  flag%) 

IF  flag%  <>  0  THEN  PRINT  "Error  in  Mode  13  digital  I/O.   Error#";  flag%:  STOP 

FOR  k  =  1  TO  2 

NEXT 

PRINT  dio%(0) 
NEXTj 
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END  SUB 

SUB  datacoilect  (t%) 

SHARED  dio%0,  c%0,  mcd%0 

md%  =  4 

dio%(0)  =  10 

dio%(2)  =  0 

dio%(3)  =  0 

dio%(l)  =  VARPTR(c%(0)) 

CALL  DASH16(md%,  VARPTR(dio%(0)),  flag%) 

IF  flag%  <>  0  THEN  PRINT  "Enor  during  Mode  4  AD  operation.  Error  #  ";  flag%:  STOP 

PRINT  t 

aa#  =  0:  bb#  =  0 

FOR  k  =  0  TO  18  STEP  2 

aa#  =  aa#  +  c%(k) 

bb#  =  bb#  +  c%(k  +  1) 

NEXTk 

aa#  =  aa#  /  10 

mcd%(t)  =  INT(aa#) 
PRINT  mcd%(t) 
END  SUB 

SUB  datasave  (t%) 

SHARED  magfield%0,  mcd%0 

'subroutine  to  write  sequential  data  file  to  disk 

PRINT  :  INPUT  "Enter  filename  for  data  file  ";  name$ 

PRINT  :  INPUT  "Enter  temperature  of  matrix.  ";  temp$ 

PRINT  :  INPUT  "Enter  maximum  field  ";  max$ 

PRINT  :  INPUT  "Enter  sensitivity  and  gain  ";  gainS 

OPEN  name$  FOR  OUTPUT  AS  #1 

WRITE  #1,  temp$,  max$,  gainS 

CLOSE  #1 

OPEN  name$  FOR  APPEND  AS  #1 

FOR  j  =  1  TO  (t  -  1) 

PRINT  j,  t,  magfield%(j).  mcd%(j) 

WRITE  #1,  magfield%(j),  mcd%(j) 

NEXTj 

CLOSE  #1 

END  SUB 

SUB  fastramp  (hifield%,  fastS) 
hidata%  =  INT(hifield%  *  (69.53  /  40000)  *  58.89) 
IF  hidata%  /  2  <>  INT(hidata%  /  2)  THEN 
hidata%  =  hidata%  +  1 
END  IF 
SELECT  CASE  fastS 
CASE  IS  =  "fastup" 

PRINT  "  Fast  ramp  up  proceeding." 
datapt%  =  0 
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DO  WHILE  datapt%  <=  hidata% 

md%  =  15 

dio%(l)  =  datapt% 

dio%(0)  =  1 

CALL  DASH16(md%,  VARPTR(dio%(0)),  flag%) 

IF  flag%  <>  0  THEN  PRINT  "Error  Mode  15  D/A  output.  Error  #  ";  flag%:  STOP 

SLEEP  1 

IF  datapt%  =  hidata%  THEN 
EXIT  DO 

END  IF 

datapt%  =  datapt%  +  10 

IF  datapt%  >  hidata%  THEN  datapt%  =  hidata% 

PRINT  datapt%,  hidata% 
LOOP 
CASE  IS  =  "fastdown" 

PRINT  "Fast  ramp  down  proceeding." 
datapt%  =  hidata% 

DO  WHILE  datapt%  >=  0 

md%  =  15 

dio%(l)  =  datapt% 

dio%(0)  =  1 

CALL  DASH16(md%,  VARPTR(dio%(0)),  flag%) 

IF  Dag%  <>  0  THEN  PRINT  "Error  Mode  15  D/A  output.  Error  #  ";  flag%:  STOP 

SLEEP  1 

IF  datapt%  =  0  THEN 
EXIT  DO 

END  IF 

datapt%  =  datapt%  -  10 

IF  datapt%  <  0  THEN  datapt%  =  0 
LOOP 
END  SELECT 

END  SUB 


SUB  fieldramp  (hifield%,  dataspace%,  directions) 
SHARED  magfield%0,  dio%0,  mcd%0>  c%0 
hidata%  =  INT(hifield%  *  (69.53  /  40000)  *  58.89) 

IF  hidata%  /  2  <>  INT(hidata%  /  2)  THEN 

hidata%  =  hidata%  +  1 

END  IF 
spacedata%  =  INT(dataspace%  *  (69.53  /  40000)  *  58.89) 

IF  spacedata%  /  2  <>  INT(spacedata%  /  2)  THEN 

spacedata%  =  spacedata%  +  1 

END  IF 
SELECT  CASE  directions 


CASE  IS  =  "up" 
datapt%  =  0 
atfield%  =  0 
t  =  1 
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DO  WHILE  datapt%  <=  hidata% 
md%  =  15 
dio%(l)  =  datapt% 
dio%(0)  =  1 

CALL  DASH16(md%,  VARPTR(dio%(0)),  flag%) 

IF  flag%  <>  0  THEN  PRINT  "Error  Mode  15  D/A  output.  Error  #  ";  flag%:  STOP 
SLEEP  1 
IF  dataspace%  <>  0  THEN 

IF  (datapt%  /  spacedata%)  =  INT(datapt%  /  spacedata%)  THEN 
SLEEP  4 

CALL  datacollect(t%) 
magfield%(t)  =  atfield% 
atfield%  =  atfield%  +  dataspace% 
t  =  t  +  1 
CALL  chart 
END  IF 
END  IF 
IF  datapt%  =  hidata%  THEN 

EXIT  DO 
END  IF 

PRINT  datapt%,  spacedata%,  t 
datapt%  =  datapt%  +  2 

IF  datapt%  >  hidata%  THEN  datapt%  =  hidata% 
LOOP 

CASE  IS  =  "down" 
datapt%  =  hidata% 
atfield%  =  hifield% 
t  =  1 

DO  WHILE  datapt%  >=  0 
md%  =  15 
dio%(l)  =  datapt% 
dio%(0)  =  1 

CALL  DASH16(md%,  VARPTR(dio%(0)),  flag%) 

IF  flag%  <>  0  THEN  PRINT  "Error  Mode  15  D/A  output.  Error  #  ";  flag%:  STOP 
SLEEP  1 

IF  dataspace%  <>  0  THEN 
IF  (datapt%  /  spacedata%)  =  INT(datapt%  /  spacedata%)  THEN 
SLEEP  4 

CALL  datacollect(t%) 
magfield%(t)  =  atfield% 
atfield%  =  atfield%  -  dataspace% 
t  =  t  +  1 
CALL  chart 
END  IF 
END  IF 
IF  datapt%  =  0  THEN 

EXIT  DO 
END  IF 

PRINT  datapt%,  spacedata%,  t 
datapt%  =  datapt%  -  2 
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IF  datapt%  <  0  THEN  datapt%  =  0 
LOOP 
END  SELECT 


END  SUB 


SUB  instruct  1 

CLS 

PRINT  :  PRINT  :  PRINT 

PRINT  "Instructions  for  connecting  magnet  for  computer  control." 

PRINT  :  PRINT 

PRINT  "1.  With  'MAINS'  off  and  'switch  heater'  disconnected  set  front  switch  " 

PRINT  "     to  constant  current  mode.  Turn  'POSITIVE'  knob  to  4  and  'NEGATIVE'" 

PRINT  "      knob  to  1." 

PRINT 

PRINT  "2.  Connect  BNC  cable  from  d/a  port  #2  to  BNC  on  back  of  power  supply.  " 

PRINT  "     Set  switch  to  external." 

PRINT 

PRINT  "3.  Turn  'MAINS'  on.   No  reading  should  appear  on  front  panel  current  " 

PRINT  "     meter.  (If  reading  is  not  zero  then  d/a  #2  is  not  set  for  zero  volts" 

PRINT  "     Turn  large  helipot  knob  fully  clockwise.  Again  no  reading  should  " 

PRINT  "     appear  on  the  meter." 

PRINT 

PRINT  "4.  Connect  switch  heater  and  energize  for  30  seconds.   You  are  now  " 

PRINT  "     ready  to  ramp  the  magnet." 

PRINT 

PRINT  "Enter  Y  to  procede" 

DO 

c$  =  INKEY$ 

LOOP  UNTIL  c$  =  "x" 

END  SUB 

SUB  nopersist 

CLS 

'Instruction  screen  for  taking  magnet  out  of  persistant  mode. 

PRINT  :  PRINT  :  PRINT 

PRINT  "Setup  for  taking  magnet  out  of  persistant  mode." 

PRINT  :  PRINT 

PRINT  "1.  Check  magnet  current  meter.  Current  should  now  register  corresponding" 

PRINT  "     to  the  current  magnetic  field." 

PRINT  :  PRINT 

PRINT  "2.  With  switch  heater  turned  'off  reconnect  lead  to  terminals." 

PRINT  :  PRINT 

PRINT  "3.  Energize  switch  heater  for  30  seconds." 

PRINT  :  PRINT 

PRINT  "Enter  'x'  to  procede" 

DO 

c$  =  INKEY$ 

LOOP  UNTIL  c$  =  "x" 
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END  SUB 

FUNCTION  pem  (wl) 

md%  =  15 

v  =  2  *  (.000105  *  (wl  /  2)  *  4095)  +  .5 

dio%(l)  =  INT(v) 

dio%(0)  =  0 

CALL  DASH16(md%,  VARPTR(dio%(0)),  flag%) 

IF  flag%  <>  0  THEN  PRINT  "Error  Mode  15  D/A  output.  Error  #  ";  flag%:  STOP 

END  FUNCTION 

SUB  persist 

CLS 

'Instructions  for  placing  superconducting  magnet  in  persistant  mode. 

PRINT  :  PRINT  :  PRINT 

PRINT  "Instructions  for  placing  magnet  in  persistant  mode." 

PRINT  :  PRINT 

PRINT  "1.  Turn  'switch  heater'  off  and  wait  30  seconds." 

PRINT 

PRINT  "2.  Disconnect  one  switch  heater  lead." 

PRINT 

PRINT  "Enter  V  to  procede" 

DO 

c$  =  INKEY$ 

LOOP  UNTIL  c$  =  "x" 

END  SUB 
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